《c语言在流体力学中的应用实例》由会员分享，可在线阅读，更多相关《c语言在流体力学中的应用实例（3页珍藏版）》请在金锄头文库上搜索。
1、#include#include#define PI 3.1415926#define FOMAT1 printf(Please input the lds beginning value as ld1= ld2= :n);#define FOMAT2 scanf(ld1=%lfld2=%lfn,&ld11,&ld21);#define FOMAT3 printf(The answers are:n);#define FOMAT4 printf(ld1=%lf v1=%lfnld2=%lf v2=%lfnQ=%lf n=%dn,ld1,V1,ld2,V2,Q,n);#define FOMAT5
2、 printf(e1=%lf Re1=%lf k1=%lfne2=%lf Re2=%lf k2=%lfn,e1,fRe(fv1(ld1,ld2),d1),k1,e2,fRe(fv2(ld1,ld2),d2),k2);double Ce=0.5,L1=300.0,d1=0.6,K1=0.0015,L2=240.0,d2=0.9,K2=0.0003,V=1.0e-6,H=6.0,g=9.807;double fv1(double x1,double x2) double c1,c2,c3,v; c1=L1/d1; c2=L2/d2*pow(d1/d2,4); c3=Ce+pow(1-pow(d1/
3、d2,2),2)+pow(d1/d2,4); v=sqrt(2*g*H/(c1*x1+c2*x2+c3); return v;double fv2(double x1,double x2) double c1,c2,c3,v; c1=L1/d1*pow(d2/d1,4); c2=L2/d2; c3=Ce*pow(d2/d1,4)+pow(1-pow(d2/d1,2),2)+1; v=sqrt(2*g*H/(c1*x1+c2*x2+c3); return v;double fld(double k,double Re) double ld; ld=0.0055*(1+pow(20000*k+1.
4、0E6/Re,1/3.0); return ld;double fRe(double v,double d) double Re; Re=v*d/V; return Re;double fe(double x1,double x2) double e; e=fabs(x1-x2)/x2); return e;void main() double k1,k2; double ld11=0,ld12=0,ld21=0,ld22=0,v11=0,v12=0,v21=0,v22=0,Re1=0,Re2=0, e1=0,e2=0,ld1=0,ld2=0,V1=0,V2=0,Q=0; int n=1; k
5、1=K1/d1; k2=K2/d2; FOMAT1 FOMAT2 v11=fv1(ld11,ld21); Re1=fRe(v11,d1); ld12=fld(k1,Re1); v21=fv2(ld12,ld21); Re2=fRe(v21,d2); ld22=fld(k2,Re2); e1=fe(ld11,ld12); e2=fe(ld21,ld22); while(e1=0.001)|(e2=0.001) ld11=ld12; ld21=ld22; v11=fv1(ld11,ld21); Re1=fRe(v11,d1); ld12=fld(k1,Re1); v21=fv2(ld12,ld21); Re2=fRe(v21,d2); ld22=fld(k2,Re2); e1=fe(ld11,ld12); e2=fe(ld21,ld22); n=n+1; ld1=ld12; ld2=ld22; V1=fv1(ld1,ld2); V2=fv2(ld1,ld2); Q=V1*PI*pow(d1,2)/4.0; FOMAT3 FOMAT4 FOMAT53 / 3
