//Code for 3D heat conduction using flux based approach //Created by G.B. Arivazhagan. Thanks to Prof. Atul Sharma who taught me CFDHT. #include #include #include #include #include using namespace std; double ***x,***y,***z,***xx,***xy,***xz,***Qxold,***Qyold,***Qzold; double ***qx,***qy,***qz,***T,***q,***qold,***Told,***unstead; void printit(int); int Imax=66; //Max grid points int Jmax=66; int Kmax=258; int imax=Imax-1; //To denote the flux indices int jmax=Jmax-1; int kmax=Kmax-1; int i,j,k; int main() { int ii,ij,ik,iter;//have used two indices; imax,jmax for flux and Imax and Jmax for Temperature for convenience double dx,dy,dz,dt,Lx,Ly,Lz; //grid properties Lx=1.0; Ly=1.0; Lz=1.0; iter=0; //To keep track of iteration imax=Imax-1; //To denote the flux indices jmax=Jmax-1; kmax=Kmax-1; double rho=7750; //Physical parameters double cp=500; double kk=16.2; double qvolx=0; //Volumetric heat generation in x and y direction double qvoly=0; double qvolz=0; clock_t t1,t2; t1=clock(); // double unstead[kmax-1][imax-1][jmax-1]; //Unsteadiness used for convergence double alpha; //Courant number double unsteadi=1; x=new double**[Kmax]; y=new double**[Kmax]; z=new double**[Kmax]; xx=new double**[kmax]; xy=new double**[kmax]; xz=new double**[kmax]; Qxold=new double**[kmax]; Qyold=new double**[kmax]; Qzold=new double**[kmax]; qold=new double**[kmax]; q=new double**[kmax]; unstead=new double**[kmax-1]; for(k=0;k

3D colour code