Hi everyone, I am trying to create my mesh with gsmh (a channel with a sphere). I am relying on the .ceo file from this link (https://github.com/jorgensd/dolfinx_ipcs/tree/a356280a0aa49d439187a52ce13b1a4bdef0eda2) but when making the modifications to create my sphere, I get an error. Can someone guide me on this?
This is my code:
SetFactory("OpenCASCADE");
lc = 0.01;
Box(1) = {0, 0, 0, 14.0, 4.0, 5.0};
Sphere(2) = {2.5,2.0,2.5, 1.0};
BooleanDifference(3) = { Volume{1}; Delete; }{ Volume{2}; Delete; };
Mesh.ElementOrder = 1;
Field[1] = Distance;
Field[1].FacesList = {7};
Field[2] = Threshold;
Field[2].IField = 1;
Field[2].LcMin = lc;
Field[2].LcMax = 15*lc;
Field[2].DistMin = 0.025;
Field[2].DistMax = 1;
Field[7] = Min;
Field[7].FieldsList = {2};
Background Field = 7;
Physical Volume("Fluid", 5) = {3};
Physical Surface("Inlet", 1) = {1};
Physical Surface("Outlet", 2) = {6};
Physical Surface("Walls", 3) = {5, 3, 2, 4};
Physical Surface("Obstacle", 4) = {7};