mac cormack in 3-D!

austin meyer

Hello!

When doing a MacCormack predictor and corrector step in 3-D, everyone will recognize these vectors:

xflt F[5]={r[c]*u[c] ,

p[c]+r[c]*sqr(u[c]),

r[c]*u[c]* v[c] ,

r[c]*u[c]* w[c] ,

(p[c]+U[c][4])*u[c] };

xflt G[5]={r[c]*v[c] ,

r[c]*v[c]* u[c] ,

p[c]+r[c]*sqr(v[c]),

r[c]*v[c]* w[c] ,

(p[c]+U[c][4])*v[c] };

xflt H[5]={r[c]*w[c] ,

r[c]*w[c]* u[c] ,

r[c]*w[c]* v[c] ,

p[c]+r[c]*sqr(w[c]),

(p[c]+U[c][4])*w[c] };

Now, for the 2-D case, we can find our predicted or corrected "Fstar" and "Gstar" as follows:

for(eq=0;eq<5;eq++)

{

Fs[c][eq]=yj[c][j_d]*F[eq]-xj[c][j_d]*G[eq];

Gs[c][eq]=xi[c][i_d]*G[eq]-yi[c][i_d]*F[eq];

}

BUT WHAT SHOULD THOSE LINES ABOVE LOOK LIKE FOR THE 3-D CASE?

thanks!

austin