Apame: 3D panel code question
Hi
If someone have looked at the Matlab or Fortran implementation of Apame 3D panel code (3dpanelmethod.com/), please help me with the following question. For the following part of code, WHY the free-stream velocity component in transversal direction "gp" is NOT taken into account? Is this because there is no side-slip angle? The vectors u, p, o are unit vectors in longitudinal, transversal and perpendicular direction, while qu, qp, qo is the induced velocity in longitudinal, transversal and perpendicular direction. ---------------------- % induced velocity in perpendicular direction qo = (p1(i,j)*o1(i,j) + p2(i,j)*o2(i,j) + p3(i,j)*o3(i,j))*(qp - (u1(i,j)*p1(i,j) + u2(i,j)*p2(i,j) + u3(i,j)*p3(i,j))*qu) % free-stream velocity components gu, go in longitudinal and perpendicular direction gu = u1(i,j)*Vx + u2(i,j)*Vy + u3(i,j)*Vz; go = o1(i,j)*Vx + o2(i,j)*Vy + o3(i,j)*Vz; % resultant velocity vx(i,j) = (-qu + gu)*u1(i,j) + (-qo + go)*o1(i,j); vy(i,j) = (-qu + gu)*u2(i,j) + (-qo + go)*o2(i,j); vz(i,j) = (-qu + gu)*u3(i,j) + (-qo + go)*o3(i,j); v(i,j) = realsqrt(vx(i,j)^2+vy(i,j)^2+vz(i,j)^2); ---------------------- If I have to consider the side-slip, does the following correction make sense? ---------------------- gu = u1(i,j)*Vx + u2(i,j)*Vy + u3(i,j)*Vz; gp = p1(i,j)*Vx + p2(i,j)*Vy + p3(i,j)*Vz; go = o1(i,j)*Vx + o2(i,j)*Vy + o3(i,j)*Vz; vx(i,j) = (-qu + gu)*u1(i,j) + (-qp + gp)*p1(i,j) + (-qo + go)*o1(i,j); vy(i,j) = (-qu + gu)*u2(i,j) + (-qp + gp)*p2(i,j) + (-qo + go)*o2(i,j); vz(i,j) = (-qu + gu)*u3(i,j) + (-qp + gp)*p3(i,j) + (-qo + go)*o3(i,j); ---------------------- Any help in this regard is appreciated. Regards |
I think nobody is interested in panel codes any more! :rolleyes:
|
All times are GMT -4. The time now is 17:58. |