Warning: Undefined array key 1 in /home/np29546/public_html/elmerice/wiki/inc/auth.php on line 78

Warning: Cannot modify header information - headers already sent by (output started at /home/np29546/public_html/elmerice/wiki/inc/auth.php:78) in /home/np29546/public_html/elmerice/wiki/inc/auth.php on line 431

Warning: Cannot modify header information - headers already sent by (output started at /home/np29546/public_html/elmerice/wiki/inc/auth.php:78) in /home/np29546/public_html/elmerice/wiki/inc/actions.php on line 38
mesh [Elmer/Ice Wiki]

Warning: Undefined array key -1 in /home/np29546/public_html/elmerice/wiki/inc/html.php on line 1458

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
mesh [2013/04/17 23:36]
gag
mesh [2020/05/02 12:55] (current)
fgillet [2D mesh from a closed contour]
Line 1: Line 1:
-====== Meshing tools ====== +====== Mesh Generation ====== 
-Due to the small aspect ratio of ice-sheets, but also glaciers, creating a mesh for a glacier geometry or an ice-sheet geometry cannot generally be done using classical meshers. Our experience in modelling different glaciers, ice-caps and ice-sheets, is that a vertically structured mesh is very often the best solution. Such a mesh can be constructed by first meshing the footprint of the glacier contour and then extruding this footprint vertically. +Due to the small aspect ratio of ice-sheets, but also glaciers, creating a mesh for a glacier geometry or an ice-sheet geometry cannot generally be done using classical meshers. Our experience in modelling different glaciers, ice-caps and ice-sheets, is that a vertically structured mesh is very often the best solution. Such a mesh can be constructed by first meshing the footprint of the glacier contour and then extruding this footprint vertically. The extrusion can be done either before running the simulation or using the internal mesh extrusion feature in Elmer, as presented below.
  
 +
 +==== Meshing Tools ====
 This Section presents some tools that have been developed and are used to construct a mesh for a glacier type geometry. This Section presents some tools that have been developed and are used to construct a mesh for a glacier type geometry.
  
Line 10: Line 12:
 |[[mesh:extrudemesh|ExtrudeMesh]] <sup>1)</sup>    | 3D      | 2D footprint  | (x,y,z) points<sup>2)</sup>         |''1/r^m'' inverse distance approximation | |[[mesh:extrudemesh|ExtrudeMesh]] <sup>1)</sup>    | 3D      | 2D footprint  | (x,y,z) points<sup>2)</sup>         |''1/r^m'' inverse distance approximation |
 |[[mesh:mshdem|MshGlacierDEM]]   | 3D      | 3D footprint x 1m   |DEMs         |bilinear | |[[mesh:mshdem|MshGlacierDEM]]   | 3D      | 3D footprint x 1m   |DEMs         |bilinear |
-|[[mesh:mshsynthetic|MshGlacierSynthetic]]   | 3D      | 3D footprint x 1m   |functions         |N/A |+|[[mesh:mshsynthetic|MshGlacierSynthetic]]   2D & 3D      | 2D or 3D footprint x 1m   |functions         |N/A | 
  
 <sup>1)</sup> possibility to include boundary layer at top or bottom <sup>1)</sup> possibility to include boundary layer at top or bottom
 <sup>2)</sup> randomly distributed supporting points <sup>2)</sup> randomly distributed supporting points
  
-Elmer also allows to perform internal extrusion during the simulation. This feature is documented [[mesh:structuredmesh|here]]. The [[courses:courses|material course]] presented in April 2013 at Edmonton takes into account this new feature+Other usefull tools for mesh generation are [[mesh:gmsh|GMSH]], a three-dimensional finite element mesh generator and the fully automatic adaptive isotropic surface remeshing tool [[mesh:yams|YAMS ]]. 
 +==== 2D mesh from a closed contour ==== 
 + 
 +Meshing tools to create the mesh of a 2D domain from a closed contour are available  
 + [[https://github.com/ElmerCSC/elmerfem/tree/elmerice/elmerice/Meshers/GIS|here]]. 
 + 
 +The contour can be given as a shapefile. The elmer  mesh (elements and boundary elements) can be converted 
 +to shapefiles to be visualised in a GIS software (e.g. QGIS). 
 + 
 +==== Structured Mesh Facilities ==== 
 + 
 +Elmer proposes special solvers to perform efficient mesh and variables manipulations when the mesh is structured. Elmer also allows to perform internal extrusion during the simulation to create such structured mesh from a 2D-footprintThese features for structured meshes are documented [[mesh:structuredmesh|here]]. The [[courses:courses|material course]] presented in April 2013 at Edmonton takes into account some of these features. More information can also be found in the [[http://www.csc.fi/english/pages/elmer/documentation|Elmer documentation]]
    
 +
 +==== Mesh Adaptation ====
 +
 +Many problems require to have a high mesh resolution in some particular locations (e.g. close to the grounding line, to capture few-kms-width ice-streams in large-scale ice-sheet simulations etc..., see e.g. Gillet-Chaulet et al., 2012). Having a fine mesh resolution everywhere in the model domain is often not computationnally affordable. Adaptive mesh refinement (AMR) is a method where the accuracy of the solution is controlled by spatially adapting the mesh size. The main difficulty is to find efficient and reliable estimators of the numerical error to control the mesh size.
 +
 +Here, we use an error estimate based on the interpolation error, proposed by Frey and Alauzet (2005). Routines to compute the anisotropic metric defining the mesh size have been implemented as Elmer/Ice Solvers. The mesh adaptation is performed using the freely available library MMG ([[http://www.mmgtools.org]], Dapogny et al., 2014).  
 +
 +Our implementation is actually restricted to the adaptation of plane-view 2D meshes comprised of linear 3-nodes triangular elements. This procedure was previously performed using external codes and the remeshing software Yams.  The mesh adaptation features are documented [[mesh:meshadaptation|here]].
mesh.1366241807.txt.gz · Last modified: 2013/04/17 23:36 by gag
CC Attribution-Share Alike 4.0 International
www.chimeric.de Valid CSS Driven by DokuWiki do yourself a favour and use a real browser - get firefox!! Recent changes RSS feed Valid XHTML 1.0