This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Last revision Both sides next revision | ||
mesh:meshadaptation [2017/07/21 11:56] tzwinger [Mesh Adaptation (2D)] |
mesh:meshadaptation [2018/03/06 16:08] fgillet [Mesh Adaptation (2D)] |
||
---|---|---|---|
Line 1: | Line 1: | ||
====== Mesh Adaptation (2D) ====== | ====== Mesh Adaptation (2D) ====== | ||
This page describes how to use the tools dedicated to mesh adaptation using the Mmg remeshing library (http:// | This page describes how to use the tools dedicated to mesh adaptation using the Mmg remeshing library (http:// | ||
+ | |||
+ | **< | ||
+ | ** | ||
+ | **Since elmerice rev. f6873f2(6/ | ||
Using the remeshing solver require to install Mmg and provide the direction to the header files and library using the CMake arguments | Using the remeshing solver require to install Mmg and provide the direction to the header files and library using the CMake arguments | ||
Line 48: | Line 52: | ||
Computing second derivatives in linear elements in not straightforward. Following Buscaglia and Dari (1997), this is done in two steps: | Computing second derivatives in linear elements in not straightforward. Following Buscaglia and Dari (1997), this is done in two steps: | ||
* First, we compute the nodal gradient < | * First, we compute the nodal gradient < | ||
- | * Second, we compute the hessian matrix //**H**// by solving the diffusive equation < | + | * Second, we compute the hessian matrix //**H**// by solving the diffusive equation < |
==== Metric intersection ==== | ==== Metric intersection ==== | ||
Line 68: | Line 72: | ||
The [[solvers: | The [[solvers: | ||
+ | |||
+ | ==== Transient mesh adaptation ==== | ||
+ | |||
+ | During a transient simulation the progression of a physical phenomena, //e.g.// the location of the grounding line, in the computational domain is usually not known in advance. \\ | ||
+ | If the mesh has been adapted using informations available at time //t//, the mesh may not be optimal anymore at time //t+n.dt//. Keeping a mesh sufficiently refined in areas where the physical phenomena progress requires to adapt the mesh frequently, introducing additionnal interpolation errors when solutions are transfered to the adapted meshes. \\ | ||
+ | The metric intersection formula given above allows to compute a metric map that contain informations at different time-steps. Iteratively solving the transient problem allows to adapt and refine the mesh in all regions where the phenomena progress. | ||
+ | |||
+ | The algorithm described in Alauzet //et al.// (J. Comp. Phys., 2007) has been implemented using a bash script that allows to iterate .sif files for the physical transient simulation and the mesh adaptation. Details are given {{ : | ||
==== Examples ==== | ==== Examples ==== | ||
- | An example for isotropic mesh adaptation can be found under [ELMER_TRUNK]/ | + | * An example for isotropic mesh adaptation can be found under |
- | Examples for anisotropic mesh adaptation can be found under [ELMER_TRUNK]/ | + | * [ELMER_TRUNK]/ |
+ | | ||
+ | * [ELMER_TRUNK]/ | ||
+ | * [ELMER_TRUNK]/ | ||
+ | * An example for transient | ||
+ | * [ELMER_TRUNK]/ |