Hi, rnrnFirst of all, this is an amazing project! Congratulations to all people that some how collaborated with it. rnrnI am a new user of the OpenLB and I just finished to run one example simulation (cylinder2d). From the books that I read, using the momentum exchange approach, we can calculate the drag coefficient using the following equation: cd = (2*Fx)/(rho*D*(u^2)). Where D = diameter of the cylinder, rho = density, u = LBM velocity. rnrnThe example simulation (Re=20) gave me Drag = 5.63043, but using the equation above I can not reach a reasonable drag coefficient for Re=20. Could somebody help me?rn rnBest,rnPedro Paulo.
I assumed that the 5.63043 (drag given by the code) was the Fx in the equation. This clearly did not gave me a good (reasonable) result when I also assumed the diameter of the cylinder in number of links/nodes and u = 0.02. rn
Hi, Liliana, Pedro!rnrnI am intrested in drag force too. I simulate the fluid flow in a rectangular channel with sphere in its center. So if
drag given by the code is already the drag coefficient
, it should be equal to the 0.47 for Re = 10 000 (https://en.wikipedia.org/wiki/Drag_coefficient). But in my simulation it’s about 5.0, but it is not a constant value (fluid velocity and density are constant in this simulation). Can you explain the reason of such difference?rnAlso can you tell me, how OpenLB calculates the Reynolds number?rnrnThanks,rnAlina.
Dear Alina,rnrnsome years ago we simulated the flow around a falling ball (cf. http://optilb.org/openlb/others 3d falling ball). We got good results for Re=10,100,1000 and 5000 with D3Q19 BGK-Boltzmann model and also for Re=10000 and 100000 with D3Q19 Smago BGK-Boltzmann model.rnrnMaybe your boundary conditions are not right or your resulotion is not fine enough.rnrnBestrnMathias