CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Main CFD Forum (https://www.cfd-online.com/Forums/main/)
-   -   Apame: 3D panel code question (https://www.cfd-online.com/Forums/main/91429-apame-3d-panel-code-question.html)

 pm3d August 10, 2011 08:21

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

 pm3d August 12, 2011 03:05

I think nobody is interested in panel codes any more! :rolleyes:

 All times are GMT -4. The time now is 02:45.