CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   Main CFD Forum (http://www.cfd-online.com/Forums/main/)
-   -   mac cormack in 3-D! (http://www.cfd-online.com/Forums/main/6760-mac-cormack-3-d.html)

austin meyer October 29, 2003 10:34

mac cormack in 3-D!
 
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


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