Calculation of Drag, Wall shear stress, Lift

Hello

I have a simulation around the airfoil the field files generated as .f00000 and a binder file .nek5000. I am trying to see if it is possible to extract surface of the airfoil as shown in the video below or calculate the following:

  1. Wall shear stress and the profile from the wall if possible
  2. Calculation of lift and drag forces

How to calculate aerodynamic forces with Paraview (youtube.com)

There was a previous discussion of calcuation of Cp vs x/c in the post
Manipulation of field files - ParaView Support - ParaView

Example field files : Re10000

Thank you
Pranay

I loaded the data set you posted. What I see is a slab of data with a hole in the middle that I expect is the cross section of your airfoil.

In the video you reference, the speaker extracts a block of the surface to compute the shear-wall stress on. Your data only has one block, so that won’t work. I don’t know if you can direct Nek5000 to write out this block, but you should be able to extract it in ParaView without too much trouble.

While looking at your data, I noticed that your mesh is not well connected. That is, when you have two hexahedra sharing a face, they are not sharing the same points. That is going to mess up ParaView’s neighbor calculations. You can fix this by running the Clean to Grid filter.

You can get the polygons around that hole in your data by applying the Extract Surface filter (after Clean to Grid). After you apply that filter, things will look the same because you are also extracting the outer boundary of the mesh. You can remove the outer boundary that you don’t want by applying the Extract Cells By Region filter.

In the Cells By Region filter, change the Intersect With parameter to Box. Then, adjust the planes of the box to be inside the boundary. Check the box that says Extract intersected.

Here is a statefile of what I did:
extract-airfoil-surface.pvsm (647.4 KB)

With this surface you should be able to follow the video’s instructions of computing normals and calculating the forces.

2 Likes

Thanks a lot, will look into it.
Appreciate the quick response

Another question, is it possible to extract wall shear stress around the surface or does it have to come from NEK5000 as we do it in OpenFOAM?

Maybe. I don’t know anything about computing wall shear stress myself. There are several references for trying to do this:

That said, it might be better to have Nek5000 compute the wall shear stress if you can. Several people reported problems with computing it post processing including inaccurate velocities on the surface due to boundary conditions.