NHT 第1章(续).pptVIP

下载本文档

174
0
约 60页
2017-01-26 发布于河南
举报
版权申诉

NHT 第1章(续).ppt

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

十百占消佰端厕绣敖锻没蕴抢死脓逞祁捶销掳窗伊龄移鞭圃操哪冈百人砾NHT 第1章(续)NHT 第1章(续) DO 103 J=2,M2 !边界系数设置 aW FLOW=ARX(J)*U(2,J)*RHO(1,J) DIFF=ARX(J)*GAM(1,J)/(XCV(2)*SX(J)) CALL DIFLOW !计算系数A(|P|) AIM(3,J)=ACOF+AMAX1(0.,FLOW) !系数aW ? DO 103 I=3,L2 IF(I.EQ.L2) GO TO 104 FL=U(I,J)*(FX(I)*RHO(I,J)+FXM(I)*RHO(I-1,J)) FLP=U(I+1,J)*(FX(I+1)*RHO(I+1,J)+FXM(I+1)*RHO(I,J)) FLOW=ARX(J)*0.5*(FL+FLP) !通过界面的流量 DIFF=ARX(J)*GAM(I,J)/(XCV(I)*SX(J)) GO TO 105 ?104 FLOW=ARX(J)*U(L1,J)*RHO(L1,J) !边界流量 DIFF=ARX(J)*GAM(L1,J)/(XCV(L2)*SX(J)) !边界扩导DW ?105 CALL DIFLOW ！A(|P|) AIM(I+1,J)=ACOF+AMAX1(0.,FLOW) ! 系数 aW AIP(I,J)=AIM(I+1,J)-FLOW ! 系数关系浅涕认缀蔫坑脯豁沫题嚼午另侈楷阵宏秋匪诚坯挚垛攘幸和掀戏咙诫览悠NHT 第1章(续)NHT 第1章(续) IF(J.EQ.M2) GO TO 106 FL=XCVI(I)*V(I,J+1)*(FY(J+1)*RHO(I,J+1)+FYM(J+1)*RHO(I,J)) FLM=XCVIP(I-1)*V(I-1,J+1)*(FY(J+1)*RHO(I-1,J+1)+FYM(J+1)* 1 RHO(I-1,J)) ? GM=GAM(I,J)*GAM(I,J+1)/(YCV(J)*GAM(I,J+1)+YCV(J+1)*GAM(I,J)+ 1 1.0E-30)*XCVI(I) GMM=GAM(I-1,J)*GAM(I-1,J+1)/(YCV(J)*GAM(I-1,J+1)+YCV(J+1)* 1 GAM(I-1,J)+1.E-30)*XCVIP(I-1) DIFF=RMN(J+1)*2.*(GM+GMM) GO TO 107 106 FL=XCVI(I)*V(I,M1)*RHO(I,M1) FLM=XCVIP(I-1)*V(I-1,M1)*RHO(I-1,M1) DIFF=R(M1)*(XCVI(I)*GAM(I,M1)+XCVIP(I-1)*GAM(I-1,M1))/YDIF(M1) 107 FLOW=RMN(J+1)*(FL+FLM) CALL DIFLOW ！A(|P|) AJM(I,J+1)=ACOF+AMAX1(0.,FLOW) !系数 aS＋1 AJP(I,J)=AJM(I,J+1)-FLOW !系数关系蜒荷盘搔捕度臼擎稿因宜丑桌砍醇泊敝幂帘欣映涝鼓印抠氢芒拍琅卤采砧NHT 第1章(续)NHT 第1章(续) VOL=YCVR(J)*XCVS(I) !速度控制容积体积 APT=(RHO(I,J)*XCVI(I)+RHO(I-1,J)*XCVIP(I-1)) 1/(XCVS(I)*DT) ! 非稳态项ρ/△t；DT：时间步长△t； AP(I,J)=AP(I,J)-APT ！等号后的AP（I,J) 为SP CON(I,J)=CON(I,J)+APT*U(I,J) AP(I,J)=(-AP(I,J)*VOL+AIP(I,J)+AIM(I,J)+AJP(I,J)+AJM(I,J)) 1/RELAX(NF) !计算aP CON(I,J)=CON(I,J)*VOL+REL*AP(I,J)*U(I,J) ！REL＝1－ DU(I,J)=VOL/(XDIF(I)*SX(J)) DU(I,J)=DU(I,J)/AP(I,J) ！速度修正值计算是中的de 103 CONTINUE DO 130 K=1,6 DO 130 J=JST,M2 DO 130 I=IST,L2 儿侩柱激匝诌痹哨屡醒棠汀鲍弗梯禾哼考砖笔辉消搀襟宛提壬彪