Cuboid Geometry-Drag and Lift
OpenLB – Open Source Lattice Boltzmann Code › Forums › on OpenLB › General Topics › Cuboid Geometry-Drag and Lift
- This topic has 5 replies, 2 voices, and was last updated 1 year, 7 months ago by Adrian.
November 3, 2021 at 11:07 am #6145
I would like to ask if I could calculate drag and lift forces on my cuboid geometry. Is it possible to calculate it? I tried implementing the Drag setup as per the cylinder2d/3d setup but I am getting infinite value for my cuboid geometry. Is it only that drag and lift can be calculated for spherical bodies?
ThankyouNovember 3, 2021 at 2:31 pm #6146
No, the functors used to calculate lift and drag in e.g. cylinder(2,3)d support arbitrary shapes.
However, I am unsure what you mean by “drag of cuboid geometry”. The
CuboidGeometryclass describes the block decomposition of the simulation domain and not the obstacle geometry of which you want to calculate flow coefficients.November 3, 2021 at 2:36 pm #6148
Thank-you for your prompt reply.
When I tried implementing the drag and lift for my inlet and outlet which are in cuboid shape,
SuperLatticePhysDrag3D<T,DESCRIPTOR> drag( sLattice, superGeometry, 3, converter ); // for inlet
SuperLatticePhysDrag3D<T,DESCRIPTOR> drag2( sLattice, superGeometry, 4, converter ); // for outlet
drag( dragA, input1 );
clout << “; drag_inlet=” << dragA << “; lift_inlet=” << dragA << endl;
drag2( dragB, input2 );
clout << “; drag_outlet=” << dragB << “; lift_outlet=” << dragB << endl;
I am getting inf value for drag_inlet and drag_outlet and same for lift as well. Could you tell me the reason for this inf value?
I did not understand the meaning of the sentence you said in your comment:
The CuboidGeometry class describes the block decomposition of the simulation domain and not the obstacle geometry of which you want to calculate flow coefficients.
Thank-youNovember 8, 2021 at 10:27 am #6158
The drag functor assumes that the given material number / indicator describes the boundary of the object for which the drag and lift coefficients are to be computed. You can see this in action e.g. in our cylinder(2,3)d examples.
What I meant by my other comment is that the cuboid geometry and material numbers describe different concepts. Asking for flow coefficients of a
CuboidGeometrydoesn’t make sense in OpenLB terminology. However, rereading your initial comment: Maybe you did not mean the
CuboidGeometryclass but the flow around cuboid shapes?November 8, 2021 at 10:49 am #6161
Thank-you for the reply.
Yes,actually I mean the flow around the cuboid shapes.
As per my above comment, I have described my drag definition at the inlet and outlet. However, I am getting infinite value at the inlet and outlet.
I tried setting the drag definition even for bstep2d/3d (Cuboid Shape geometry) and found similar infinite value for drag coefficient and lift. That is the reason I asked whether Drag force is calculated using Stokes Law. If it is calculating by Stokes Law, the Radius of Curvature in the cuboid geometry would be taken to be infinite.
So, could you please help me to imply the drag and lift definition for cuboid shapes?
Thank-youNovember 8, 2021 at 11:10 am #6163
If you mean the flow around the cuboid shape that you have to provide the material number of the cuboid shapes (an this number must only describe the cuboids). I repeat that the drag functor will not work when provided in / outflow material numbers in your approach.
For these kinds of questions I also kindly refer you to our upcoming Spring School. In the practical sessions there we also offer support for getting started with your own cases.
- You must be logged in to reply to this topic.