mkgeo_obstacle(1rheolef)

📦
version 7.2 · 4 versions
, in debian trixie, ubuntu noble, ubuntu questing

💖 mkgeo_obstacle — mesh of the obstacle-in-channel 💖

okay bestie, so this is the rheolef-7.2 version of this page, and it's literally about building meshes of that classic obstacle-in-channel thing. you're gonna eat this up, i promise. 🦋

✨ NAME ✨

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

🌸 SYNOPSIS 🌸

mkgeo_obstacle [options] [nx [ny]]

🎀 EXAMPLE 🎀

so like, this command right here builds a triangle-based 2d unstructured mesh of the [-L,L]x[0,c] with a unit circle hole at origin. iconic geometry tbh.

mkgeo_obstacle > obstacle.geo
geo obstacle.geo

💅 DESCRIPTION 💅

bestie. this command is the convenient way to build a mesh for the obstacle-in-channel geometry, because it's a super classical benchmark in complex fluid flow problems. it calls the bamg unstructured mesh generator with the anisotropy feature. 🔥

the mesh files go on name.geo where name is the basename for the output (peep the -name option below, she's important).

and not just that — the three auxiliary files required for automatic mesh generation with Rheolef combined with bamg are also provided:

name.bamg
name.bamgcad
name.dmn

no cap, all three. you're welcome. 💜

🌈 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. periodt. 💅

🫶 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 (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.

ok so the no-warranty part is real, bestie, don't come for me 😭

and that's the whole vibe check — this was rheolef, Version 7.2, page mkgeo_obstacle(1rheolef). go build some slay meshes, you got this 👑💖

📦 where this page comes from