Update on the development of the OpenFOAM CFD workbench (CFDOF)

[Worufu]

Aha. The issue seems to occur when the 'digit grouping symbol' in your Windows localisation settings is set to '.' Something is obviously being printed internally as "xxx.0" and this is being interpreted as a multiplication by 10. I will dig around and see if I can find the culprit, but for now you could work around it by changing your digit grouping symbol (1000s separator) in the Windows control panel.

Oh so good. It works! Yes the problem is the italian decimal separator and grouping symbol (comma instead dot and dot instead quotation mark). In the future can this problem be written in the CFDOF readme? I think I will not the only one to have this external issue . Thanks again.

[oliveroxtoby]

Aha. The issue seems to occur when the 'digit grouping symbol' in your Windows localisation settings is set to '.' Something is obviously being printed internally as "xxx.0" and this is being interpreted as a multiplication by 10. I will dig around and see if I can find the culprit, but for now you could work around it by changing your digit grouping symbol (1000s separator) in the Windows control panel.

Oh so good. It works! Yes the problem is the italian decimal separator and grouping symbol (comma instead dot and dot instead quotation mark). In the future can this problem be written in the CFDOF readme? I think I will not the only one to have this external issue . Thanks again.

Found the error - my bad! It should be working now with your normal settings, I hope.
Thanks for reporting the issue.

[thschrader]

just succeeded with a first example
laminar_flow_sphere.png
this is a very nice way to play with open-foam. thanks for this nice add-on.

Hi looo,
is that a "half-model" you are using or is the sphere fully integrated in the cube?
I get an error when i want to simulate a cylinder in a channel.
(patch without faces in boundary-dict)
Thomas
OS: Windows 10
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.17.11939 (Git)
Build type: Release
Branch: master
Hash: 6e7952ec672895900eec0c2a25807b25befba818
Python version: 2.7.8
Qt version: 4.8.7
Coin version: 4.0.0a
OCC version: 7.1.0
Locale: German/Germany (de_DE)

cylinder_in_channel.JPG (86.68 KiB) Viewed 833 times

missing_faces.JPG (58.53 KiB) Viewed 833 times

[looo]

is that a "half-model" you are using or is the sphere fully integrated in the cube?

I used a box and a sphere and removed the sphere from the box. To select the sphere you can use a clipping plane from coin:
https://forum.freecadweb.org/viewtopic. ... 10#p170965

But I have no idea what your problem is. I made a 2d-case too and it worked out of the box.

[HoWil]

Found the error - my bad! It should be working now with your normal settings, I hope.
Thanks for reporting the issue.

Hi Oliver,
can you please explain in more detail what kind of error you corrected. There is a similar problem with the material dialogue in fem-wb (see https://forum.freecadweb.org/viewtopic. ... 3e#p187725) and I am sure your infos could help here .
Thanks in advance,
HoWil

[HoWil]

is that a "half-model" you are using or is the sphere fully integrated in the cube?

I used a box and a sphere and removed the sphere from the box. To select the sphere you can use a clipping plane from coin:
https://forum.freecadweb.org/viewtopic. ... 10#p170965

But I have no idea what your problem is. I made a 2d-case too and it worked out of the box.

In such cases you could also use the "ListSelect" developed by m42kus ... see : https://github.com/drhooves/SelectionTools
EDIT: a nice screencast can be found here https://forum.freecadweb.org/viewtopic. ... or#p181397
BR HoWil

[thschrader]

is that a "half-model" you are using or is the sphere fully integrated in the cube?

I used a box and a sphere and removed the sphere from the box. To select the sphere you can use a clipping plane from coin:
https://forum.freecadweb.org/viewtopic. ... 10#p170965

But I have no idea what your problem is. I made a 2d-case too and it worked out of the box.

In such cases you could also use the "ListSelect" developed by m42kus ... see : https://github.com/drhooves/SelectionTools
EDIT: a nice screencast can be found here https://forum.freecadweb.org/viewtopic. ... or#p181397
BR HoWil

My fault. I produced an "empty" boundary condition...
I love this WB

cylinder_channel.JPG (35.44 KiB) Viewed 792 times

[thschrader]

We probably don't have the spare time or funding to implement this any time soon, so editing case dictionaries is the best bet. You would want the 'forces' function object, see the OpenFOAM user guide - https://cfd.direct/openfoam/user-guide/ ... 390006.2.2 - and the template you would need to use can be copied from <OpenFOAM directory>/etc/caseDicts/postProcessing/forces/forcesIncompressible

Of course, if anyone wants to add this functionality to the workbench, that would be great.

No need
You can calculte the resulting forces in paraview by projecting/integrating the pressure field along
the x,y,z-Axes. I have stolen the idea from this guy:
https://www.youtube.com/watch?v=J944HOj_4b0

[thschrader]

@oliveroxtoby
Here is a little info about a calculation I have done.
The pic above (cylinder in channel) is a test case. The channel is 6 m long, 2 m high and 100 mm "thick".

Theory:
Diameter cylinder is D=400 mm. Fluid is water with rho=1000 kg/m^3, inlet speed 1,5 m/s.
This gives a pressure of q=1085 N/m^2, Re=5,9*10^5 ==> drag-coeff=0,8 according to DIN4131.
The resulting force on an infinite long cylinder is F=cf*D*q=348 N/m.

openfoam-result, see spreadsheet in paraview pic:
the resulting drag-force on the 10 cm long simulation-cylinder is 35,3 N, so for a 1m part it is 353 N!
I think that is a perfect matching!

(btw: as far as I understood the input-files, openfoam caculates a "normalized" pressure field,
p_openfoam=p_real/density.
So to get the real pressure, you must multiply the calculated field with fluid-density. OK?)
regards Thomas

OS: Windows 10
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.17.11939 (Git)
Build type: Release
Branch: master
Hash: 6e7952ec672895900eec0c2a25807b25befba818
Python version: 2.7.8
Qt version: 4.8.7
Coin version: 4.0.0a
OCC version: 7.1.0
Locale: German/Germany (de_DE)

drag_lift_forces.JPG (140.05 KiB) Viewed 741 times

[oliveroxtoby]

When will this be part of the addon manager ?
It really full featured workbench already

It is now available in the addons manager (CfdOF) - thanks @yorick.