[ivan]

Oh, thank you again for that, Robin, it really helped!rnrnBut now the challenge is: I would like to plot the pressure drop around the obstacle over time, so, I would need the pressure drops in every instant. For that, I thought of simply copy and paste the results printed on the terminal by the function [getResults], but then I realised it’s not even possible because I am not able to visualize the whole printed data due to the limited dimensions of the terminal and even if possible, would be quite exausting the copy all the data manually ( specially in some edited examples of mine, in which there are more obstacles ). Also, I would like to plot the velocity in two or three points around the obstacle, also over the time, and I simply don’t get these data from the terminal. For that, I was considering using Paraview, but for now I was only able to plot the Pressure/Velocity between two different points ( or in a line ) for a single instant, which is not what I want. Any ideas?rnrnAtt.

[ivan]

Hi, Robin,rnrnBasically I have edited the cylinder example in order to add more round obstacles to the flow and now I would like to do some analysis regarding the velocity / Reynolds number and maybe even the pressure drop. Do you think itś easier to do it with getResults function than with Paraview? And how is it done exactly?

[ivan]

Thank you for that Robin, it really helped! 🙂

[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

[ivan]

But right now I’m dealing with the cylinder2d example, in which I have to add more obstacles. Any tips of how to start it?

[ivan]

Yes, I was able to change the dimensions of the step but I wanted to move it, to make is further from the channel inlet.

[ivan]

Yes, thank you! That’s exactly what I was looking for. Still, is there any way I can also move the step in the “”bstep2d”” example and not only changing its dimensions?

[ivan]

Yup! Seems to work now! Thank you very much 🙂

[ivan]

That’s what I tried to do at first. For example, I have renamed the “”bstep2d (copy)”” folder to “”Flow&Obstacle””, then, after editing some conditions ( like the length of the step ) and trying to compile everything again with the command “”Make”” I got the following message:rnrn””depalma@zenon:~/olb-1.0r0/examples/Flow&Obstacle$ makernCreate dependencies for bstep2d.cpprn/bin/sh: 3: Obstacle/bstep2d.d: not foundrnmake: *** [/home/depalma/olb-1.0r0/examples/Flow&Obstacle/bstep2d.d] Error 127rndepalma@zenon:~/olb-1.0r0/examples/Flow&Obstacle$””rnrnIs there anything else I should change besides the name of the folder?rn

[ivan]

Basically I would like to keep the original example and at the same time have a copy, which I could edit as much as I wanted to.

[ivan]

Yes, I was able to compile the original example using the command “”Make”” and also run the bstep2d executable through the terminal and visualize the GIFs and the simulation itself with Paraview.

[ivan]

Okay, just wrote you. And thank you for that :))

[Author]

Posts