stokes-darcy.geo

This file on gitlab

The stokes-darcy.geo file given below, describes the domain for the Stokes-Darcy equation example problem. It can be used to generate the respective mesh using the following command:

gmsh stokes-darcy.geo -1 -2 -o stokes-darcy.msh

// SPDX-FileCopyrightText: Copyright © DUNE-FUFEM Project contributors, see file AUTHORS.md
// SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception OR LGPL-3.0-or-later
 
// Generate msh file with $ gmsh stokes-darcy.geo -1 -2 -o stokes-darcy.msh
 
Point(1) = {0,   0,   0};
Point(2) = {5,   0,   0};
Point(3) = {5,   2,   0};
Point(4) = {4.5, 2,   0};
Point(5) = {3.5, 1,   0};
Point(6) = {1.5, 1,   0};
Point(7) = {0.5, 2,   0};
Point(8) = {0,   2,   0};
 
 
Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Line(4) = {4, 5};
Line(5) = {5, 6};
Line(6) = {6, 7};
Line(7) = {7, 8};
Line(8) = {8, 1};
Line(9) = {7, 4};
 
// Stokes subdomain
Curve Loop(1) = {4, 5, 6, 9};
Plane Surface(1) = {1};
Physical Surface(1) = {1};
 
// Darcy subdomain
Curve Loop(2) = {1, 2, 3, 4, 5, 6, 7, 8};
Plane Surface(2) = {2};
Physical Surface(2) = {2};
 
// Stokes Dirichlet boundary
Physical Curve(1) = {9};
 
// Darcy Dirichlet boundary
Physical Curve(2) = {1};
 
Mesh.MshFileVersion = 2.2;
Mesh.MeshSizeMax = 0.5;
Mesh 2;