This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
mesh [2012/11/30 13:27] tzwinger |
mesh [2020/05/02 12:55] fgillet [2D mesh from a closed contour] |
||
---|---|---|---|
Line 1: | Line 1: | ||
- | ====== | + | ====== |
- | 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 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: | |[[mesh: | ||
|[[mesh: | |[[mesh: | ||
- | |[[mesh: | + | |[[mesh: |
< | < | ||
< | < | ||
+ | |||
+ | Other usefull tools for mesh generation are [[mesh: | ||
+ | ==== 2D mesh from a closed contour ==== | ||
+ | |||
+ | Meshing tools to create the mesh of a 2D domain from a closed contour are available | ||
+ | | ||
+ | |||
+ | 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-footprint. These features for structured meshes are documented [[mesh: | ||
+ | |||
+ | |||
+ | ==== 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:// | ||
+ | |||
+ | 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: |