Re: Changing Geometry

[ivan]

Oh, wow, that was really helpfull! Thank you very much for that, Robin!rnrnSo, basically I was able o insert some more obstacles ( for now Im only interested in the round ones, but I really appreciate your extra tips for the square ones as well ) but it seems like I’m missing something because I’m getting some ERROR messages regarding the boundary conditions. Let me get everything here:rnrnOn // Parameters for the simulation setup// I added:rnconst T secondcenterCylinderX = 0.2;rnconst T secondcenterCylinderY = 0.2+L/2;rnconst T secondradiusCylinder = 0.035;rnrnOn /// Stores geometry information in form of material numbers/// I added:rn Vector<T,2> secondcenter(secondcenterCylinderX,secondcenterCylinderY);rn IndicatorCircle2D<T> secondcircle(secondcenter, secondradiusCylinder);rnrnThen, on /// Set material number for cylinder/// I added:rnsuperGeometry.rename(1,5,secondcircle);rnrnAnd finally, on /// Material=5 –>bouzidi, which seems to be part of /// Setting of the boundary conditions, I added:rn Vector<T,2> secondcenter(secondcenterCylinderX,secondcenterCylinderY);rn IndicatorCircle2D<T> secondcircle(secondcenter, secondradiusCylinder);rnrn and after sLattice.defineDynamics(superGeometry, 5, &instances::getNoDynamics<T,DESCRIPTOR>()); I addedrn offBc.addZeroVelocityBoundary(superGeometry, 5, secondcircle);rnrnThen when running it I get a lot of those:rnrn[BoundaryConditionInstantiator2D] ERROR: no boundary found at (17,19) ~ (1.65,0.295), in direction 3rn[BoundaryConditionInstantiator2D] ERROR: no boundary found at (17,18) ~ (1.65,0.29), in direction 1rn[BoundaryConditionInstantiator2D] ERROR: no boundary found at (17,19) ~ (1.65,0.295), in direction 2rn[BoundaryConditionInstantiator2D] ERROR: no boundary found at (17,20) ~ (1.65,0.3), in direction 3rn[BoundaryConditionInstantiator2D] ERROR: no boundary found at (17,19) ~ (1.65,0.295), in direction 1rn[BoundaryConditionInstantiator2D] ERROR: no boundary found at (17,20) ~ (1.65,0.3), in direction 2rnrnBut still, it runs. What could it be?rnrn[ And by the way, I made versions with 2, 5 and 30 obstacles, and I get more repetitions of those ERROR messages the more obstacles I have ]rnrnAtt.rn