![二维椭圆型流动传热程序_第1页](http://file4.renrendoc.com/view14/M07/28/2F/wKhkGWa8a1mAVgQ2AACGt6-Pg50098.jpg)
![二维椭圆型流动传热程序_第2页](http://file4.renrendoc.com/view14/M07/28/2F/wKhkGWa8a1mAVgQ2AACGt6-Pg500982.jpg)
![二维椭圆型流动传热程序_第3页](http://file4.renrendoc.com/view14/M07/28/2F/wKhkGWa8a1mAVgQ2AACGt6-Pg500983.jpg)
![二维椭圆型流动传热程序_第4页](http://file4.renrendoc.com/view14/M07/28/2F/wKhkGWa8a1mAVgQ2AACGt6-Pg500984.jpg)
![二维椭圆型流动传热程序_第5页](http://file4.renrendoc.com/view14/M07/28/2F/wKhkGWa8a1mAVgQ2AACGt6-Pg500985.jpg)
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
PAGEPAGE17二维椭圆型流动传热通用程序变量表及算例说明(本材料仅供教学参考)西安交通大学CFD&NHT/EHT研究中心陶文铨教授西安目录………………2一、FORTRAN变量表…………3二、关于程序的主要说明………………………6二维椭圆型流动和传热问题通用计算机算法方面的特点………………6各程序的主要功能…………7三种坐标系统………………8网格系统与节点命名方法…………………9十一个例题的已知条件与求解内容…………12例题1直角坐标中二维稳态无内热源的导热……………………12例题2空心圆柱内的稳态热传导…………………12例题3正方形管道内的充分发展对流换热………13例题4内壁上有直肋的环形通道内充分发展对流换热…………13例题5给定流场条件下温度场的计算……………14例题6二维突扩通道中的流动与换热……………14例题7方形通道内的复杂充分流动………………15例题8旋转圆盘上的冲击流动……………………15例题9轴对称燃烧内的瞬间燃烧过程……………16例题10带中心射流的通道内的紊流换热………16例题11有质量源的流动问题……………………17一、FORTRAN变量表ListofFORTRANVariablesACOFQuantitycalculatedbysubroutineDIFLOWtogivethecombinedconvectionanddiffusioneffect.AIM(I,J)Thecoefficient.AIP(I,J)Thecoefficient.AJM(I,J)Thecoefficient.AJP(I,J)Thecoefficient.AP(I,J)tThecoefficient;alsoinGAMSOR.APTTheunsteadyterm.AREALocalvariable,usuallytheareaofaC.V.face.ARHOLocalvariable,(area)().ARX(J)TheareaofthemainC.V.facenormaltothexdirection.ARXJ(J)ThepartofARX(J)thatoverlapsontheC.V.forV(I,J).ARXJP(J)ThepartofARX(J)thatoverlapsontheC.V.forV(I,J+1).BLBLCBLMBLPCoefficientsusedintheblockcorrection.CON(I,J)Theconstanttermbinthediscriminationequation;alsostandsforinGAMSOR.DENOMTemporarystorage.DIFFDiffusionconductanceD.DTThetimestep.DU(I,J)InfluencingU(I,J).DV(I,J)InfluencingV(I,J).F(I,J,NF)Various.FLTemporarystorageleadingtoFLOW.FLMTemporarystorageleadingtoFLOW.FLOWFLPMassflowratethroughaC.V.face.TemporarystorageleadingtoFLOW.FV(J)FVP(J)InterpolationfactorswhichgivethemassflowAtamaingridpoint,I,JasFV(J)(I,J)+FVP(I,J)(I,J+1)FX(I)Interpolationfactors,whichgivetheinterface.FXM(I)DensityRHOM(atthelocationofU(I,J))asFX(I)RHO(I,J)+FXM(I)RHO(I-1,J).FY(J)Interpolationfactors,whichgivetheinterface.FYM(J)DensityRHOM(atthelocationofV(I,J))asFY(J)RHO(I,J)+FYM(J)RHO(I,J-1).GAM(I,J)Thediffusioncoefficient.IIndexdenotingthepositioninx.IBEGIENDTemporaryvaluesusedinPRINT.IFSTThefirstvalueofIforwhichtheprint-outisarranged;usedinPRINT.IITemporaryindex.IPREFThevalueofIforthegridpoint,whichisusedasareferenceforpressure.ISTThefirstinternal-pointvalueofI.ISTFIST-1;usedinSOLVE.ITERAcounterforiterations.IT1IT2TemporaryvaluesusedinSOLVE.JIndexdenotingthepositioniny.JFLTemporaryindexusedinPRINT.JFSTSimilartoILST.JJTemporaryindexJLSTSimilartoILSTJPREFSimilartoIPREF.JSTThefirstinternal-pointvalueofJ.JSTFJST-1;usedinSOLVE.JT1JT2TemporaryvaluesusedinSOLVE.LASTThemaximumnumberofiterationsallowedbytheuser.LBLK(NF)When.TRUE.TheblockcorrectionforF(I,J,NF)isused.LPRINT(NF)When.TRUE.,F(I,J,NF)isprinted.LSOLVE(NF)When.TRUE.,WesolveforF(I,J,NF).LSTOPWhen.TRUE.Computationstops.L1ThevalueofIforthelastgridlocationinthexdirection.L2(L1-1).L3(L1-2).MODEIndexforthecoordinatesystem;=1for,=2for,=3for.M1ThevalueofIforthelastgridlocationintheydirection.M2(M1-1).M3(M1-2).NTemporarystorageforNF.NFIndexdenotingaparticular.NFMAXThelargestvalueofNFforwhichstorageisassigned.NGAMNFMAX+3;GAM(I,J)canbeconsideredasF(I,J,NGAM).NPNFMAX+1;P(I,J)canbeconsideredasF(I,J,NP).NRHONFMAX+2;RHO(I,J)canbeconsideredasF(I,J,NRHO).NTIMES(NF)ThenumberofrepetitionsofthesweepsinSOLVEforthevariableF(I,J,NF).P(I,J)Thepressurep.PC(I,J)Thepressurecorrection.PREFThepressureatthereferencepoint.PT(I)orPT(J)QT(I)orQT(J)TransformedcoefficientsintheTDMA.R(J)TheradiusrforamaingridpointI,J.REL1.0-RELAX(NF).RELAX(NF)RelaxationfactorforF(I,J,NF).RHO(I,J)Thedensity.RHOCONThevalueofforaconstant-densityproblem.RMN(J)ThevalueofradiusrforthelocationtowhichV(I,J)refers.SMAXThelargestabsolutevalueofthe“masssource”usedintheequation.SSUMThealgebraicsumofallthe“masssources”intheequation.SX(J)ScalefactorforthexdirectionatthemaingridlocationsY(J).SXMN(J)ScalefactorforthexdirectionatinterfacelocationsYV(J).TEMPTemporarystorage.TIMETimetforunsteadyproblems.TITLE(NF)AlphabetictitleforF(I,J,NF).U(I,J)Thex-directionvelocityu.V(I,J)They-directionvelocityv.VOLVolumeoftheC.V.X(I)Thevaluesofxatgridpoints.XCV(I)Thex-directionwidthsofmainC.V..XCVI(I)ThepartofXCV(I)thatoverlapsontheC.V.forU(I,J).XCVIP(I)ThepartofXCV(I)thatoverlapsontheC.V.forU(I+1,J).XCVS(I)Thex-directionwidthofthestaggeredC.V.forU(I,J).XDIF(I)ThedifferenceX(I)-X(I-1).XLThex-directionlengthofthecalculationdomain.XU(I)ThelocationsoftheC.V.faces;i.e.thelocationofU(I,J).Y(J)Thevaluesofyatgridpoints.YCV(J)They-directionwidthsofmainC.V.YCVR(J)TheareaforamainC.V.YCVRS(J)TheareafortheC.V.forv(I,J).YCVS(J)They-directionwidthofthestaggeredC.V.forV(I,J).YDIF(J)ThedifferenceY(J)-Y(J-1).YLThey-directionlengthofthecalculationdomain.YV(J)ThelocationsoftheC.V.faces;i.e.thelocationofV(I,J).二关于程序的主要说明二维椭圆型流动和传热问题通用计算机程序算法方面的特点采用原始变量法,即以速度U、V及压力P作为直接求解的变量守恒型的差分格式,离散方程系对守恒型的控制方程通过对控制容积作积分而得出的,无论网格疏密程度如何,均满足在计算区域内守恒的条件;采用区域离散化方法B,即先定控制体界面、再定节点位置采用交叉网格,速度U、V与其他变量分别存储于三套网格系统中;不同的项在空间离散化过程中去不同的型线假设,源项采用局部线性化方法;扩散——对流项采用乘方格式(但很容易转化为中心差分、迎风差分或混合格式);街面上的扩散系数采用调和平均法,而密度与流速则用线性插值;不稳态问题采用全隐格式,以保证在任何时间步长下均可获得具有物理意义的解;边界条件采用附加源项法处理;耦合的流速与压力采用SIMPLE算法来求解;迭代式的求解方法,对非线性问题,整个求解过程具有迭代性质;对于代数方程也采用迭代法求解;采用交替方向先迭代法求解代数方程并补以块修正技术以促进收敛。各方程块的主要功能主程序MAIN:规定整个计算过程的流程,决定是否停止计算。子程序SETUP1:设置与网格系统有关的,在计算过程中保持不变的几何参数及设置U、V、P’、P、ρ、CON(SC)、AP(SP)的初值。子程序SETUP2:1.建立每一轮迭代中变量U、V及P’和其他变量离散方程的系数aE、aW、aN、aS、aP及b项。所有变量公用一套数组逐个求解,逐个确定各自目标的系数,以节省内存。在建立系数过程中调用DIFLOW子程序;2调用SOLVE子程序求解代数方程;3迭代次数指标及时间步长增值;子程序DIFLOW:计算乘方格式中的系数A(|p|);子程序SOLVE:采用交替方向线迭代法并辅以块修正技术求解代数方程;子程序GRID:设置为建立网格系统所必须的量,包括设置x,y方向求解区域的宽度XL及YL设置x,y方向上各自的节点数L1及M1设置x,y方向上控制体积的界面位置,XU(I),I=2,L1;YV(J),J=2,M1。如为均分网格可确定子程序UGRID,如为非均分网格需由用户在GRID中把XU(I),YV(J)一一设置好。规定坐标系,MODE=1,2,3分别相应于直角、圆柱轴对称及极坐标当MODE不等于1时,规定径向起始点半径R(I)子程序START:设置初值,包括(1)、对不稳态问题给出初始条件(2)、对稳态问题给出迭代求解的假定值(3)、已知的边界值也可在此块中,在设置初值或假定值时一并送入,但对边界条件随时间而异的不稳态不能在此块中赋值,因在整个计算过程中此会仅执行一次。子程序DENSE:规定密度场,对于密度为常数的问题,此块中可不设任何语句,但必须保留ENTRY及RETUIRN语句子程序BOUND:规定边界条件,包括设置各变量的边界条件对各特殊问题所需的量,如Nu,fRe等等可在此块中计算,但也可在OUTPUT中设置子程序OUTPUT:打印输出,包括每做一轮迭代,输出一行信息,以观察收敛情形;调用PRINT子程序,实现二维物理量场的输出;用户所需输出的其他特殊变量子程序GAMSOR:设置扩散系数及源项,包括(1)设置控制方程中规定的源项Sc,Sp(2)对第二,三类边界条件设置与边界相邻的控制体中的附加源项,Sc,ad及sp,ad内节点及边界节点规定扩散系数,扩散采用附加源项法时令边界扩散系数为零。12、子程序UGRID:为均分网格设置界面位置13、子程序PRINT:输出二维物理场(3)三种坐标系统x,iyjx,iyj轴对称圆柱坐标,计算对θ=1弧度的区域进行;MODE=2;Y(J)可以从任何起点算起,但R(J)须从对称轴开始;R(1)维计算区的最小半径。两个区域XL,YL相等,但R(1)不同R(1)=0R(1)=0YLXLYLXLR(1)极坐标系,Z方向厚度为1;MODE=3;计算区域的θ角须小于2π;Y(J)可从任何点起算,但R(J)比从中心点开始;R(1)为计算区域边界的最小半径。yyθrMODEX(I)Y(J)R(J)SX(J)1xy1.01.02xyr1.03θyrr(4)网格系统与节点命名方法V(i,j)U(i,j)V(i,j)U(i,j)P(i,j)P(L1,M1)U(2,M2)U(2,2)V(2,M1)V(2,2)U(L1,M2)V(L1,2)(a)总体网格系统V(L2,M1)U(L1,2)U(L1,2)U(I,J)U(I,J)I,JI-1,JV(I,J)I,JI,J-1(b)速度编号的方法XDIF(3)XDIF(3)XU(3)X(I)XDIF(2)U(3,J)U(2,J)L2L3I321XU(2)X(1)X(2)X(3)XU(L1)X(L1)XU(I)XDIF(I)XCVS(I)XCVS(3)XCV(I)XCVIP(I)XCVI(I)XCVS(L2)(c)x方向的主要几何参数U(L1,J)U(L1,J)L1L1ARXJ(J)ARXJP(J)ARXJ(J)ARXJP(J)ARX(J)ycv(J-1)1y(1)yv(2)2y(2)yv(3)3y(3)M1y(M1)yv(M1)M2M3V(I,J)I,JYDIF(J)ycvs(j)(e)y方向的主要几何参数U(2,J)XCVS(3)(2,J)U(3,J)XDIF(3)(d)速度U(3,J)的控制容积(3,J)(3,J)P(1,2)P(1,1)P(1,2)P(1,1)P(2,2)P(2,1)(f)压力参考点(1,1)上压力的确定P(1,1)=P(1,2)+P(2,1)-P(2,2)变量内点下标起始值变量ISTJSTФ,P’,P22U32V23三、十一个例题的已知条件与求解内容例题1直角坐标中二维稳态无内热源的导热已知:已知:求解区域如图1所示,四个边界上的温度由下式决定;T=x+y+xy求:该区域中的温度分布。yyx21图1例题2空心圆柱内的稳态热传导已知:求解区域如图2所示。左边界为给定温度,右边界为对流,上边界为绝热,下边界为给定热流。整个计算区域内都有源项。求:该区域中的温度分布。计算方法:用附加源项法处理第2,3类边界条件图2图2x211S=100-0.5TK=1.0K=0.2(1+T/100)T=100(1+y)绝热yq=50例题3正方形管道内的充分发展对流换热已知:常物性流体在均匀壁温的正方形截面管道内作充分发展的对流换热(如图3)。求:截面上的速度分布、温度分布、fRe及Nu。计算方法:充分发展对流换热问题的处理特点图3图311TwTwyxz例题4内壁上有直肋的环形通道内充分发展对流换热已知:常物性流体在带直肋的环形通道内作充分发展层流对流换热,内表面为均匀壁温(周向),轴向则呈线性变化,外表面绝热。R1=1,R2=2,а=15o。求:截面上速度分布、温度分布及fRe、Nu。计算方法:耦合问题的一种处理方法。ααR1R2图4例题5给定流场条件下温度场的计算已知:流体流过一个直角的两表面,速度场为u=Ax,V=-AY,a=10,Tin=500,Tw=100求计算区域中温度的分布。图5图5ToutTin11Tw绝热xy例题6二维突扩通道中的流动与换热已知:由两平行平
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 数学身高比较北师大版
- 三年级上册人教版数学真题
- 四年级苏教版语文练习培养良好的学习习惯
- 人类食物与生态关系教学北师大版
- 新北师大版直角坐标系解读
- 初中数学苏教版要点解析
- 部编版人教版一年级下册科学教案动物的食性与消化
- 创新人教版字词教学策略
- 汉语音韵学练习与解析
- 掌握人教版第十三章内能复习攻略
- SYB创业培训课件-10步全
- 鞍山金和矿业有限公司台上山铁矿(铁矿)矿山地质环境保护与土地复垦方案评审意见
- 现代厨房管理教案
- 评价中层干部优缺点三篇
- 铜陵顶讯科技有限公司高导铜基电镀线项目环境影响报告书
- 阿尔茨海默症药物治疗 阿尔茨海默症药物治疗(老年安全用药课件)
- 劳动创造美好生活第一章
- 2023年网络安全知识竞赛题库及答案(共290题)
- 新媒体运营PPT完整全套教学课件
- xx矿山救护队各种记录表格
- 变压器水喷雾系统试喷方案
评论
0/150
提交评论