资源描述：
《相变传热与流体流动数值分析作业》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、相变传热与流体流动数值分析作业6相变传热与流体流动数值分析作业6学院（系）：能源与动力学院专业：能源与环境工程学生姓名：王佳琪学号：21110060指导教师：完成日期：大连理工大学DalianUniversityofTechnologyPoiseuilleFlowusingLBMSubjectsConsidertheproblemsoftheconstantpressuregradientdrivenconstantPoiseuilleflowin2DpipelineusingD2Q9model.Analyticalsolution:Simulationc
2、onditions:1.LinewidthisH=50,lengthisL=200,inletpressureof1.01,outletpressureof1.0,fluidmotionviscosity1/6.2.Boundaryconditions:theupanddownreboundboundaryformat,ortheleftandlightboundary(noequilibriumextrapolationformatorZou/Heboundary)SolutionThefollowingisthegeneralprograms://王佳
3、琪-作业6.cpp:定义控制台应用程序的入口点。//#include"stdafx.h"#include<stdio.h>#include<math.h>#include<string.h>#include<stdlib.h>#defineQ9#defineNX100//40#defineNY20#defineP_in1.01//1.001#defineP_out1.0inte[Q][2]={{0,0},{1,0},{0,1},{-1,0},{0,-1},{1,1},{-1,1},{-1,-1},{1,-1}
4、};doublew[Q]={4.0/9,1.0/9,1.0/9,1.0/9,1.0/9,1.0/36,1.0/36,1.0/36,1.0/36};doublep[NX+1][NY+1],u[NX+1][NY+1][2],u0[NX+1][NY+1][2],fp[NX+1][NY+1][Q],Fp[NX+1][NY+1][Q],U[NX+1][NY+1];inti,j,k,ip,jp,n;//,krever[Q]doublec,dx,dy,dt,rho0,p0,LX,LY,tau_f,niu,error;//dtpvoidinit();doublefpeq(
5、intk,doublep,doubleu[2]);voidcollision();voidboundary();voidstreaming();voidmacro();voidoutput(intm);voidError();intmain(){init();for(i=0;i<=NX;i++)for(j=0;j<=NY;j++){U[i][j]=(P_in-P_out)*((LY/2)*(LY/2)-(LY/2-j*LY/NY)*(LY/2-j*LY/NY))/(2*LX*niu);}for(n=0;;n++){collision();str
6、eaming();macro();boundary();if(n%100==0){Error();printf("The%dthcomputationresult:",n);printf("Theu,vofpoint(NX/2,NY/2)is:%f,%f",u[NX/2][NY/2][0],u[NX/2][NY/2][1]);printf("Themaxrelativeerrorofuvis:%f",error);if(n>=1000){if(n%1000==0)output(n);if(error<1.0e-8)break;//跳
7、出循环}}}output(n);for(j=0;j<=NY;j++)printf("U[0][%d]=%f",j,U[0][j]);return0;}voidinit(){dx=1.0;dy=1.0;LX=dx*NX;LY=dy*NY;niu=0.025;//运动粘度系数dt=dx;c=dx/dt;tau_f=0.5+3*niu/dt;rho0=1.0;p0=(rho0*c*c)/3;printf("dt=%f,c=%f,p0=%f",dt,c,p0);for(i=0;i<=NX;i++)for(j=0;j<=NY;j++){u[
8、i][j][0]=0.0;u[i][j][1]=0.0;p[i][
