This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
mesh [2017/07/18 09:00] fgillet [Structured Mesh Facilities] |
mesh [2020/05/02 12:55] (current) fgillet [2D mesh from a closed contour] |
||
---|---|---|---|
Line 13: | Line 13: | ||
|[[mesh: | |[[mesh: | ||
|[[mesh: | |[[mesh: | ||
+ | |||
< | < | ||
Line 18: | Line 19: | ||
Other usefull tools for mesh generation are [[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 ==== | ==== Structured Mesh Facilities ==== | ||
Line 24: | Line 33: | ||
==== Mesh Adaptation ==== | ==== 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: |