Does the forces.f90 file require modifications for 3D cases ?

[airwarriorg91]

Hello,
I am trying to validate the flow around a sphere at Re=3700. I have used a coarse domain of 5D x 3D x 3D with discretization (321, 321, 120) and slip boundary conditions in y-directions whereas periodic b.c. in z-direction (I have attached the input.i3d file for your reference). inputi3d.txt

The issue is the values of CL and CD, I am getting doesn't feel right especially if I use stretching in y-direction (iibm=2 for stretching case) while keeping everything else same. Does the force.f90 file require any modifications to calculate the CL and CD for 3D cases ? From what I understood, the subroutine is calculating the CL and CD for each plane in z direction and averaging it.

Without Stretching:

With Stretching:

Thank you for the help and suggestions !

[slaizet]

Hi,
The current implementation in the master branch is designed for extruded shapes in the Z-direction (like the cylinder), so you would need to change the code to accommodate for the volume control surfaces in the Z-direction.
Cheers,
Sylvain

[airwarriorg91]

Hello !
Thank you for the response. Can you please share the source for the method of computation of the forces ? I am not able to find the reference mentioned in the forces.f90 file.

[slaizet]

Have a look at this report for the equations but you just need to do in Z-direction what is done in the X- and Y-direction.
Best,
Sylvain
HOESSLER.pdf

[airwarriorg91]

Thank you for the reference ! I will have a look at it.

[airwarriorg91]

Hello @slaizet !
I tried modified the force.f90 file based on the equations provided in the reference. This is the result I got from computing the lift and drag using the already available code (2D CV) and using the new code I added (3D CV). The simulation is for a flow around cylinder. The input file used is 'input_test.i3d'.

Thanks a lot ! I wish to contribute to xcompact3d by completing the code and making it available for others users too. So, I would need some guidance from your side.

Edit: Earlier, I had faced some issues related to parallelization which I solved. Now, the code works fine on both my local PC and HPC.