dvm force calc

gistfile1.c

    vi = w.vortices[step];

    fx = -vi.gamma * (vi.y - vi.yI)/w.dt;
    fy = -vi.gamma * (vi.xI - vi.x)/w.dt;

    /* average with previous shed vortex */
    if (step >= 1) {
      vi = w.vortices[step - 1];
      fx = (fx + (-vi.gamma * (vi.y - vi.yI)/w.dt))/2.;
      fy = (fy + (-vi.gamma * (vi.xI - vi.x)/w.dt))/2.;
    }

    for (i=0; i<=step;i++) {
      vi = w.vortices[i];
      vj = w.displaced_vortices[i];

      u = (vj.xI - vi.xI)/w.dt;
      v = (vj.yI - vi.yI)/w.dt;
      
      fy += vi.gamma * u;
      fx += vi.gamma * v;
    }

    w.vortices[step].fx = fx;
    w.vortices[step].fy = fy;