fsimicropillar.geo 1.2 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758
  1. Mesh.ScalingFactor = 1e-6;
  2. tcfine = 0.9;
  3. Mesh.CharacteristicLengthMin = 0.5;
  4. Mesh.CharacteristicLengthMax = 4;
  5. SetFactory("OpenCASCADE");
  6. Rectangle(1) = {0, 0, 0, 350, 120, 0};
  7. Point(5) = {100, 0, 0, 100};
  8. Point(6) = {100, 60, 0, tcfine};
  9. Point(7) = {102, 60, 0, tcfine};
  10. Point(8) = {102, 62, 0, tcfine};
  11. Point(9) = {104, 60, 0, tcfine};
  12. Point(10) = {104, 0, 0, 100};
  13. Line(5) = {5, 6};
  14. Line(6) = {9, 10};
  15. Line(9) = {5, 10};
  16. Circle(7) = {6, 7, 8};
  17. Circle(8) = {8, 7, 9};
  18. Line Loop(2) = {1, 2, 3, 4};
  19. Line Loop(3) = {5, 7, 8, 6, -9};
  20. Plane Surface(2) = {3};
  21. Point(11) = {150, 120, 0, 100};
  22. Point(12) = {150, 60, 0, tcfine};
  23. Point(13) = {152, 60, 0, tcfine};
  24. Point(14) = {152, 58, 0, tcfine};
  25. Point(15) = {154, 60, 0, tcfine};
  26. Point(16) = {154, 120, 0, 100};
  27. Line(10) = {11, 12};
  28. Line(11) = {11, 16};
  29. Line(12) = {16, 15};
  30. Circle(13) = {12, 13, 14};
  31. Circle(14) = {14, 13, 15};
  32. Line Loop(4) = {10, 13, 14, -12, -11};
  33. Plane Surface(3) = {4};
  34. Coherence;
  35. fluid = 1; pillar = 2; inlet = 3; outlet = 4; sides = 5; clamp = 6;
  36. Physical Surface(fluid) = {4};
  37. Physical Surface(pillar) = {2,3};
  38. Physical Line(inlet) = {1};
  39. Physical Line(outlet) = {8};
  40. Physical Line(sides) = {2,7,9,14};
  41. Physical Line(clamp) = {15,16};