NAME

       mkgeo_obstacle - mesh of the obstacle-in-channel (rheolef-7.1)

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.