---
name: "mkgeo_obstacle"
section: "1rheolef"
platform: "linux"
package: "rheolef"
version: "7.2-4"
release: "Debian 13 (Trixie)"
platforms: ["linux"]
url: "https://new.carta.tech/mkgeo_obstacle"
---

# mkgeo\_obstacle(1rheolef)

mesh of the obstacle-in-channel (rheolef-7.2)

## NAME

mkgeo\_obstacle - mesh of the obstacle-in-channel (rheolef-7.2)

## SYNOPSIS

mkgeo\_obstacle \[options\] \[nx \[ny\]\]

## EXAMPLE

The following command build a triangle-based 2d unstructured mesh of the \[-L,L\]x\[0,c\] with a unit circle hole at origin.

```
    mkgeo_obstacle > obstacle.geo
    geo obstacle.geo
```

## DESCRIPTION

This command is convenient for building a mesh for the obstacle-in-channel geometry, as it is a very classical benchmark in complex fluid flow problems. It calls `bamg` unstructured mesh generator with anisotropy feature. The mesh files goes on name.geo where *name* is the basename for the output (see option `-name` below). The three auxiliary files required for automatic mesh generation with Rheolef combined with `bamg` are also provided:

```
        name.bamg
        name.bamgcad
        name.dmn
```

## THE GEOMETRY

`-c` *float* \
 `-L` *float*

 These options control the geometry parameters. Default is c=2 and L=2\*c.

`-cartesian` \
 `-zr`

 These options control the geometry coordinate system. Default is Cartesian.

`-quarter` \
 `-half`

 When the `quarter` option is activated, the upstream-downstream symmetry is assumed and only one quarter of the geometry is considered. Otherwise, with the `half` option, the upstream-downstream geometry is generated Default is the `quarter` option.

`-N` *int*

 The number of elements around the half of the circular obstacle. When the `quarter` option is activated, only a quarter \
\
  of the circular obstacle is considered. Default is `N=100`.

## THE DISCRETIZATION

`-hmin` *float*

 Controls the edge length on the obstacle. Default is `hmin`=0.1. Changing the density by the previous options applies the factor 1/n to `hmin` also.

## BOUNDARY DOMAINS

The boundary sides are represented by domains: `obstacle`, `axis`, `wall`, `upstream` and `downstream`.

## OTHERS OPTIONS

- `-name` *string*

   Set the basename for the output files. By default, the basename is `obstacle`.

`-[no]split`

 Split each triangle in three subtriangles by inserting the barycenter as an additional node. This is useful for using the Scott and Vogelius incompressible mixed finite element P2-P1d for velocity-pressure approximation. Default is no split.

`-[no]clean`

 Clear temporary files (this is the default).

`-[no]verbose`

 In verbose mode, print to stderr all commands and logs.

## IMPLEMENTATION

This documentation has been generated from file main/sbin/mkgeo\_obstacle.sh

## AUTHOR

Pierre  Saramito  \<Pierre.Saramito@imag.fr\>

## COPYRIGHT

Copyright   (C)  2000-2018  Pierre  Saramito  \<Pierre.Saramito@imag.fr\> GPLv3+: GNU GPL version 3 or later  \<http://gnu.org/licenses/gpl.html\>. This  is  free  software:  you  are free to change and redistribute it. There is NO WARRANTY, to the extent permitted by law.