一维对流扩散问题TDMA算法.docVIP

April 24, 2005 一维对流扩散问题TDMA算法 Filed under: 乱七八糟 一维对流扩散问题TDMA算法，支持中心差分，迎风，混合，指数，乘方格式，还可以输出精确解 用法直接编译运行，不带参数运行时会自动输出使用方法: upwind –left=(左端点处的值) –right=(右端点处的值) –pdelta=(网格Pe数)–scheme=(格式，可选center/upwind/mixed/exponent/power/precise) –output=(输出文件名) 如：upwind –left=100 –right=200 –pdelta=5 –scheme=upwind –output=Pe5Upwind.dat参数顺序无关紧要。 注意现在的代码输出相对值(phi-phi0)/(phi_L-phi0)，想输出绝对值phi的自己改代码吧:)另存为upwind.c编译即可 #include stdio.h #include math.h int main(int argc, char **argv){ double max(double number1,double number2); void parseParameter(int *argc, char **argv, double *PHI_0, double *PHI_L, double *P_Delta, char *scheme, char *outputfilename); const int GRIDAMOUNT=20;//网格数目为20，由于两端点已知，实际未知数个数为19个。 double PHI_0=100; double PHI_L=200; double P_Delta=1; //网格Pe数，取为1,5,10等。 char scheme[10]=upwind; int PRECISE=0; char outputfilename[100]=output.dat; double a_E[GRIDAMOUNT-1]; double a_W[GRIDAMOUNT-1]; double a_P[GRIDAMOUNT-1]; double A_P_Delta[GRIDAMOUNT-1]; double P_TDMA[GRIDAMOUNT-1]; double Q_TDMA[GRIDAMOUNT-1]; double PHI[GRIDAMOUNT-1]; int i; FILE *dataoutput; parseParameter(argc, argv, PHI_0, PHI_L, P_Delta, scheme, outputfilename); if(!strcmp(scheme,upwind)){for (i=0;iGRIDAMOUNT-1;i++){A_P_Delta[i]=1;} } else if(!strcmp(scheme,center)) {for (i=0;iGRIDAMOUNT-1;i++){A_P_Delta[i]=1-0.5*fabs(P_Delta);} } else if(!strcmp(scheme,mixed)) {for (i=0;iGRIDAMOUNT-1;i++){A_P_Delta[i]=max(0,1-0.5*fabs(P_Delta));} } else if(!strcmp(scheme,exponent)) {for (i=0;iGRIDAMOUNT-1;i++){A_P_Delta[i]=fabs(P_Delta)/(exp(fabs(P_Delta))-1);} } else if(!strcmp(scheme,power)) {for (i=0;iGRIDAMOUNT-1;i++){A_P_Delta[i]=max(0,pow((1-0.1*fabs(P_Delta)),5));} } else if(!strcmp(scheme,precise)) {PRECISE=1;} else{fprintf(stderr,\nUnknown scheme %s!\nAccepted schemes are: upwind|center|mixed|exponent|power|precise\n,scheme);exit(2); } //A_P_Delta[i]=1-0.5*fabs(P_Delta); //A_P_Delta[i]=1; //A_P_Delta[i]=max(0,1-0.5*fabs(P_Delta));