FSI on a rigid 3D flap

[sfraniatte]

Dear community,

I am interested in performing an FSI simulation involving a rigid 3D flap (available as an STL file) immersed in a turbulent flow. The flap is mounted on a rotational axis and subjected to a resisting torque.

While reviewing the examples provided in the fsi folder, I noticed that there does not seem to be a direct functionality for adding such a resisting torque.

My question is therefore: is it feasible to perform this simulation using the existing capabilities of OpenLB, or would additional development work be required?

Thanks in advance !

Best regards,
Sylvain

[Adrian]

Currently the idea is that OpenLB provides the framework for efficiently embedding structures into porosities and obtaining their boundary integrals while the user plugs in their own structure model. For a rotating flap with resisting torque I would consider the flap as a point mass with orientation and inertia. You can start with the rigid valve 2d case and add such a torque there to test (~ 1 line of code)

[sfraniatte]

Thank you for this quick answer ! I noticed that only the taylorCouette3d example works properly on GPU though. The rigidValve2d example does not compile and the deformableParticleInShearFlow3d runs but all the fields remain equal to zero. It is not a big deal though.

Best regards,
Sylvain

[Adrian]

Rigid valve definitely works on GPU on my system – can you post the error?

Deformable particle doesn’t work in the public version on GPUs but will be ported in the next release.

In terms of FSI there is also the GPU enabled centrifugal pump case.

[sfraniatte]

Interesting. I have donloaded the last release on GitLab to retry if it works and the compilation still does not work with this error :

../../../src/fsi/platform/cpu/integral.h(77): error: no instance of function template “lambda [](olb::meta::id<FIELD>, const olb::FieldD<T, DESCRIPTOR, FIELD> &)->auto” matches the argument list
argument types are: (olb::meta::id<olb::fields::fsi::ELEMENT_TORQUE>, const olb::concepts::placeholder::Parameters2D::value_t)
object type is: lambda [](olb::meta::id<FIELD>, const olb::FieldD<T, DESCRIPTOR, FIELD> &)->auto
((set(meta::id<FIELDS>{}, fields.template get<FIELDS>())), …);
^

I have CUDA 12.9 with an Ampere generation GPU. The others cases work properly and it runs properly with gcc and mpic++.

I did not know for the centrifugal pump case. Thank you !

[Author]

Posts