Please Delete
mesh->Ustar[i][j][k] = mesh->U[i][j][k]
- dtdx*(0.25*((mesh->U[i+1][j][k]+mesh->U[i][j][k])*(mesh->U[i+1][j][k]+mesh->U[i][j][k]))
- 0.25*((mesh->U[i][j][k]+mesh->U[i-1][j][k])*(mesh->U[i][j][k]+mesh->U[i-1][j][k])))
- dtdy*(0.25*(mesh->U[i][j][k]+mesh->U[i][j+1][k])*(mesh->V[i][j][k]+mesh->V[i+1][j][k])
- 0.25*(mesh->U[i][j][k]+mesh->U[i][j-1][k])*(mesh->V[i][j-1][k]+mesh->V[i+1][j-1][k]))
- dtdz*(0.25*(mesh->U[i][j][k]+mesh->U[i][j][k+1])*(mesh->W[i][j][k]+mesh->W[i+1][j][k])
- 0.25*(mesh->U[i][j][k]+mesh->U[i][j][k-1])*(mesh->W[i][j][k-1]+mesh->W[i+1][j][k-1]))
+ dtdxx*(mesh->U[i-1][j][k]-2*mesh->U[i][j][k]+mesh->U[i+1][j][k])
+ dtdyy*(mesh->U[i][j-1][k]-2*mesh->U[i][j][k]+mesh->U[i][j+1][k])
+ dtdzz*(mesh->U[i][j][k-1]-2*mesh->U[i][j][k]+mesh->U[i][j][k+1]);
mesh->Vstar[i][j][k] = mesh->V[i][j][k]
- dtdy*(0.25*((mesh->V[i][j+1][k]+mesh->V[i][j][k])*(mesh->V[i][j+1][k]+mesh->V[i][j][k]))
-0.25*((mesh->V[i][j][k]+mesh->V[i][j-1][k])*(mesh->V[i][j][k]+mesh->V[i][j-1][k])))
- dtdx*(0.25*(mesh->V[i][j][k]+mesh->V[i][j+1][k])*(mesh->U[i][j][k]+mesh->U[i][j+1][k])
-0.25*(mesh->V[i][j][k]+mesh->V[i-1][j][k])*(mesh->U[i-1][j][k]+mesh->U[i-1][j-1][k]))
- dtdz*(0.25*(mesh->V[i][j][k]+mesh->V[i][j+1][k])*(mesh->W[i][j][k]+mesh->W[i][j][k+1])
-0.25*(mesh->V[i][j][k]+mesh->V[i][j][k-1])*(mesh->W[i][j][k-1]+mesh->W[i-1][j][k-1]))
+ dtdxx*(mesh->V[i-1][j][k]-2*mesh->V[i][j][k]+mesh->V[i+1][j][k])
+ dtdyy*(mesh->V[i][j-1][k]-2*mesh->V[i][j][k]+mesh->V[i][j+1][k])
+ dtdzz*(mesh->V[i][j][k-1]-2*mesh->V[i][j][k]+mesh->V[i][j][k+1]);
Rennlist Member

mesh->Ustar[i][j][k] = mesh->U[i][j][k]
- dtdx*(0.25*((mesh->U[i+1][j][k]+mesh->U[i][j][k])*(mesh->U[i+1][j][k]+mesh->U[i][j][k]))
- 0.25*((mesh->U[i][j][k]+mesh->U[i-1][j][k])*(mesh->U[i][j][k]+mesh->U[i-1][j][k])))
- dtdy*(0.25*(mesh->U[i][j][k]+mesh->U[i][j+1][k])*(mesh->V[i][j][k]+mesh->V[i+1][j][k])
- 0.25*(mesh->U[i][j][k]+mesh->U[i][j-1][k])*(mesh->V[i][j-1][k]+mesh->V[i+1][j-1][k]))
- dtdz*(0.25*(mesh->U[i][j][k]+mesh->U[i][j][k+1])*(mesh->W[i][j][k]+mesh->W[i+1][j][k])
- 0.25*(mesh->U[i][j][k]+mesh->U[i][j][k-1])*(mesh->W[i][j][k-1]+mesh->W[i+1][j][k-1]))
+ dtdxx*(mesh->U[i-1][j][k]-2*mesh->U[i][j][k]+mesh->U[i+1][j][k])
+ dtdyy*(mesh->U[i][j-1][k]-2*mesh->U[i][j][k]+mesh->U[i][j+1][k])
+ dtdzz*(mesh->U[i][j][k-1]-2*mesh->U[i][j][k]+mesh->U[i][j][k+1]);

JD
mesh->Vstar[i][j][k] = mesh->V[i][j][k]
- dtdy*(0.25*((mesh->V[i][j+1][k]+mesh->V[i][j][k])*(mesh->V[i][j+1][k]+mesh->V[i][j][k]))
-0.25*((mesh->V[i][j][k]+mesh->V[i][j-1][k])*(mesh->V[i][j][k]+mesh->V[i][j-1][k])))
- dtdx*(0.25*(mesh->V[i][j][k]+mesh->V[i][j+1][k])*(mesh->U[i][j][k]+mesh->U[i][j+1][k])
-0.25*(mesh->V[i][j][k]+mesh->V[i-1][j][k])*(mesh->U[i-1][j][k]+mesh->U[i-1][j-1][k]))
- dtdz*(0.25*(mesh->V[i][j][k]+mesh->V[i][j+1][k])*(mesh->W[i][j][k]+mesh->W[i][j][k+1])
-0.25*(mesh->V[i][j][k]+mesh->V[i][j][k-1])*(mesh->W[i][j][k-1]+mesh->W[i-1][j][k-1]))
+ dtdxx*(mesh->V[i-1][j][k]-2*mesh->V[i][j][k]+mesh->V[i+1][j][k])
+ dtdyy*(mesh->V[i][j-1][k]-2*mesh->V[i][j][k]+mesh->V[i][j+1][k])
+ dtdzz*(mesh->V[i][j][k-1]-2*mesh->V[i][j][k]+mesh->V[i][j][k+1]);
The Best Porsche Posts for Porsche Enthusiasts
BTW - when I posted this........
I was just kidding, OK???

You guys still scare me.... Seriously...
JD


