Please Delete
#154
Almost Deleted
Rennlist Member
Rennlist Member
If anyone can tell me what this code is, I'll post another Storm Troopers cameo
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->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]);
#155
Shameful Thread Killer
Rennlist Member
Rennlist Member
#157
Banned
If anyone can tell me what this code is, I'll post another Storm Troopers cameo
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->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
#158
Administrator - "Tyson"
Lifetime Rennlist
Member
Lifetime Rennlist
Member
I'll give you a hint.....
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]);
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]);