Reply To: Wall Shear Stress on Bouzidi Wall
Due to recent bot attacks we have changed the sign-up process. If you want to participate in our forum, first register on this website and then send a message via our contact form.
› Forums › OpenLB › General Topics › Wall Shear Stress on Bouzidi Wall › Reply To: Wall Shear Stress on Bouzidi Wall
December 11, 2023 at 1:02 pm
#8011
TSchenkel
Participant
Thanks.
ad 1: That is what I figured. I checked that by setting up the STLReader and modified the functor to output the normals (will need those for other postprocessing). They did indeed all point towards the centre.
ad 2: How difficult would it be to set up such an analytical functor? Just starting to get my head round OpenLB (and object-oriented C++). Where would I start looking for something to use as a template for such a functor?
Re 3D WSS vector, does the following change to latticePhysWallShearStress3D.hh look correct to output the WSS vector in 3D in the global coordinate system?
l49:
SuperLatticePhysF3D<T, DESCRIPTOR>(sLattice, converter, 3), // changed this from 1 to 3 to give WSS vector
l155ff:
T WSS[3];
WSS[0] = traction[0] - tractionNormalComponent[0];
WSS[1] = traction[1] - tractionNormalComponent[1];
WSS[2] = traction[2] - tractionNormalComponent[2];
// magnitude of the wall shear stress vector
// output[0] = util::sqrt(WSS[0]*WSS[0] + WSS[1]*WSS[1] + WSS[2]*WSS[2]);
//
// Output all 3 components of WSS vector
output[0]= WSS[0];
output[1]= WSS[1];
output[2]= WSS[2];