




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、问题描述图1为20℃的水在长度为150mm,宽度为10mm的管道中流动,流入管道速度为0.2m/s,流出管道背压为1atm,基于simple算法对整个流场进行计算,计算管长100mm处流速并与fluent计算结果对比。图1流动系统流动计算:该流动问题为二维定常无内热源不可压缩层流流动。控制方程离散在二维直角坐标系中,对流—扩散方程的通用形式为:图2网格编号针对本问题,其连续性方程和动量方程为:交错网格下动量方程的离散:上式积分得到:对流项采用一阶迎风格式,扩散项采用中心差分,MATLAB程序:%%%%%%%%%%%%%%%%%%%%%%%%dx=5e-4;dy=5e-4;den=998;dyna=1001.6e-6;a=0.01/dy+3;%x节点%b=0.15/dx+1;%y节点%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%边界条件的设置、初始值的设置%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%T=0.2*ones(a,b);U1_old=zeros(a,b);P1_old=zeros(a,b);%初值%U1_old(3:a-2,:)=T(3:a-2,:);%加速收敛,边界条件%uw_sum(1,b)=0;ue_sum(1,b)=0;de1_sum(1,b)=0;m=1;A(1:a-4,2:b-1)=0.2;%初值为1%e22=1;whilee22>0.01%整个流场迭代%ifm>4breakendm=m+1;forj=2:1:b-1;ifj>2breakendn=1;%迭代次数%e=1;%初次迭代%B=[0;0;0.0823455;0.149211;0.196652;0.224194;0.23674;0.241087;0.242138;0.242214;0.242125;0.242083;0.242125;0.242214;0.242138;0.241087;0.23674;0.224195;0.196652;0.149212;0.0823459;0;0];U1_old(3:a-2,1)=B(3:a-2,1);U1_old(3:a-2,2)=B(3:a-2,1);whilee>0.001%迭代,算两列,其他为已知%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%对(i,j)点求uxin%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%fori=3:1:a-2Fe1(i,j)=den*(U1_old(i,j)+U1_old(i,j+1))*dy/2;%对流强度%Fw1(i,j)=den*(U1_old(i,j-1)+U1_old(i,j))*dy/2;ap_u1(i,j)=max(Fe1(i,j),0)+max(-Fw1(i,j),0)+4*dyna;%离散方程对应系数%ae_u1(i,j)=dyna+max(-Fe1(i,j),0);aw_u1(i,j)=dyna+max(Fw1(i,j),0);an_u1(i,j)=dyna;as_u1(i,j)=dyna;A1_u(i,i-1)=-an_u1(i,j);A1_u(i,i)=ap_u1(i,j);A1_u(i,i+1)=-as_u1(i,j);b1_u(i,j)=aw_u1(i,j)*U1_old(i,j-1)+ae_u1(i,j)*U1_old(i,j+1)...+(P1_old(i,j)-P1_old(i,j+1))*dy/2;%(i,j+1)点的源项%endA1_u(1,1)=1;%附加点%A1_u(2,2)=1;%边界点%A1_u(a-1,a-1)=1;A1_u(a,a)=1;b1_u(1,j)=0;%附加点%b1_u(2,j)=0;%边界点%b1_u(a-1,j)=0;b1_u(a,j)=0;U1_new(:,j)=inv(A1_u)*b1_u(:,j);%得到(i,j+1)点的速度计算值%%对(i,j)点求p%uw_sum(1,j)=0;ue_sum(1,j)=0;de1_sum(1,j)=dy/ap_u1(2,j);fori=3:1:a-2uw_sum(1,j)=uw_sum(1,j)+U1_old(i,j-1);ue_sum(1,j)=ue_sum(1,j)+U1_new(i,j);de1_sum(1,j)=de1_sum(1,j)+dy/ap_u1(i,j);endP1_fix(1,j)=(uw_sum(1,j)-ue_sum(1,j))/de1_sum(1,j);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%对(i,j)点进行速度修正%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%fori=3:1:a-2U1_fix(i,j)=(dy/ap_u1(i,j))*P1_fix(1,j);%速度修正%endU1_old(3:a-2,j)=U1_new(3:a-2,j)+U1_fix(3:a-2,j);P1_old(2:a-1,j)=P1_old(2:a-1,j)+0.4*P1_fix(1,j);ee=uw_sum(1,j)-ue_sum(1,j);e=max(max(abs(ee)));%判断的是否合理?%n=n+1;N(1,j)=n;ifn>1000breakendA((a-4)*m+1:(a-4)*(m+1),2:b-1)=U1_old(3:a-2,2:b-1);%赋值,即m从1开始%e
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 社区孤独症患者关怀的个人方法计划
- 分配任务与跟进进度的秘书工作计划
- 品牌资产与市场风险的关联研究计划
- 全面提升团队绩效的年度规划计划
- 强化职能部门之间的配合计划
- 如何选择合适的品牌推广渠道计划
- 2025年物位仪项目建议书
- 2025年数字货币金融项目发展计划
- 2025年沥青试验仪器项目发展计划
- 新员工入职引导及工作流程简明教程
- 音乐学科阅读方案
- 2024-2030年中国医药设备市场发展分析及市场趋势与投资方向研究报告
- 基于新课标学习的教材解读及教学建议部编《道德与法治》二年级下册
- 《劳动教育》+劳动项目八《栽种小树》课件
- 《社区康复》课件-第二章 社区康复的内容
- 泪道狭窄与阻塞的护理
- 铣床工安全技术操作规程培训
- 宏观报告:美国房地产市场的五大特点
- 视频监控系统验收测试报告
- 医院骨科专病数据库建设需求
- 湖北省普通高中2022-2023学年高一下学期学业水平合格性考试模拟化学(二)含解析
评论
0/150
提交评论