channel.geo 1.0 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344
  1. boxlen = 2.0;
  2. boxwidth = 0.8;
  3. cylinderradius = 0.07;
  4. cylinderxpos = 0.25;
  5. // Tiny misalignment:
  6. cylinderypos = boxwidth/2 + boxwidth/1000;
  7. lc = 0.02;
  8. lcfine = 0.002;
  9. Point(1) = {0,0,0, lc};
  10. Point(2) = {boxlen,0,0, lc};
  11. Point(3) = {boxlen,boxwidth,0, lc};
  12. Point(4) = {0,boxwidth,0, lc};
  13. Point(5) = {cylinderxpos,cylinderypos,0, lcfine};
  14. Point(6) = {cylinderxpos+cylinderradius,cylinderypos,0, lcfine};
  15. Point(7) = {cylinderxpos,cylinderypos+cylinderradius,0, lcfine};
  16. Point(8) = {cylinderxpos-cylinderradius,cylinderypos,0, lcfine};
  17. Point(9) = {cylinderxpos,cylinderypos-cylinderradius,0, lcfine};
  18. Line(1) = {1, 2};
  19. Line(2) = {2, 3};
  20. Line(3) = {3, 4};
  21. Line(4) = {4, 1};
  22. Circle(5) = {6, 5, 7};
  23. Circle(6) = {7, 5, 8};
  24. Circle(7) = {8, 5, 9};
  25. Circle(8) = {9, 5, 6};
  26. Line Loop(1) = {4, 1, 2, 3};
  27. Line Loop(2) = {6, 7, 8, 5};
  28. Plane Surface(1) = {1, 2};
  29. Plane Surface(2) = {2};
  30. fluid = 1; cylinder = 2; inlet = 3; outlet = 4;
  31. Physical Surface(fluid) = 1;
  32. Physical Surface(cylinder) = 2;
  33. Physical Line(inlet) = 4;
  34. Physical Line(outlet) = 2;