Both sides previous revision
Previous revision
Next revision
|
Previous revision
|
solvers:glads [2018/10/27 18:28] gag [Examples] |
solvers:glads [2019/09/06 13:17] (current) gagliar [Known bugs] |
==== Known bugs ==== | ==== Known bugs ==== |
If in parallel a moulin belongs to two partitions, the flux from the moulin is taken into account twice. This is may be not a bug, but the partitioning should avoid to have moulins at the boundary of partitions. The python tool ''makemoulin.py'' (in ''[ELMER_TRUNK]/elmerice/Meshers'') take care to have no duplicated moulins on the boundaries of partitions. | If in parallel a moulin belongs to two partitions, the flux from the moulin is taken into account twice. This is may be not a bug, but the partitioning should avoid to have moulins at the boundary of partitions. The python tool ''makemoulin.py'' (in ''[ELMER_TRUNK]/elmerice/Meshers'') take care to have no duplicated moulins on the boundaries of partitions. |
| |
| If running the solver on a 3d internally extruded mesh, one should specify ''Preserve Edges = True'' and ''Preserve Baseline = True'' in the Simulation section. |
| |
| If running on a "true" 3d mesh, the ''GlaDSchannelsSolver'' has to be executed on the 3d body (not the bottom surface body as for the two other GlaDS solvers). |
| |
| |
| |
* ''GlaDSchannelOut'' has two functions: declare that the ''Channel Area'' variable is an edge variable (''Element = "n:0 e:1"'') and create output vtu files for edge variables. | * ''GlaDSchannelOut'' has two functions: declare that the ''Channel Area'' variable is an edge variable (''Element = "n:0 e:1"'') and create output vtu files for edge variables. |
| |
Since version Version 8.3 (Rev: b213b0c8), ''GlaDSchannelOut'' works for parallel simulation (no more vtk or acscii output, only vtu). These solvers only work in transient. They can be executed either on a 2d plane view mesh defining the bedrock or on the boundary of a 3d mesh. More details about the specificity of the solvers are given below. | Since version Version 8.3 (Rev: b213b0c8), ''GlaDSchannelOut'' works for parallel simulation (no more vtk or acscii output, only vtu). These solvers only work in transient. They can be executed either on a 2d plane view mesh defining the bedrock or on the boundary of a 3d mesh. If using internal extrusion within Elmer see the [[mesh:structuredmesh|structured mesh]] page for essential keywords to preserve baseline and edges. More details about the specificity of the solvers are given below. |
| |
| |
Usage: | Usage: |
<code> | <code> |
python makempoulin.py --meshdir mesh_dir --moulin moulin_file --partition number_of_partition | python makemoulin.py --meshdir mesh_dir --moulin moulin_file --partition number_of_partition |
</code> | </code> |
where ''moulin_file'' is an ascii file which contains the (x,y) coordinates of the moulins. The same file has to be used in gmsh so that nodes with the moulin coordinates already exist. | where ''moulin_file'' is an ascii file which contains the (x,y) coordinates of the moulins. The same file has to be used in gmsh so that nodes with the moulin coordinates already exist. |
Examples using the //GlaDS// Solver can be found in ''[ELMER_TRUNK]/elmerice/Tests/''. | Examples using the //GlaDS// Solver can be found in ''[ELMER_TRUNK]/elmerice/Tests/''. |
| |
==== Reference ==== | ==== References ==== |
The description of the GlaDS model is in:\\ | The description of the GlaDS model is in:\\ |
* Werder M.A., I.J. Hewitt, C.G. Schoof and G.E. Flowers, 2013. Modeling channelized and distributed subglacial drainage in two dimensions. Journal of Geophysical Research: Earth Surface, 118(4), 2140-2158.\\ | * Werder M.A., I.J. Hewitt, C.G. Schoof and G.E. Flowers, 2013. Modeling channelized and distributed subglacial drainage in two dimensions. Journal of Geophysical Research: Earth Surface, 118(4), 2140-2158.\\ |