GROUP 1. Run title and other preliminaries TEXT(Non-uniform moving-wall velocity TITLE DISPLAY In-Form is used for specifyinhg non-uniform inlet, outlet and wall boundary conditions, for: steady laminar flow and heat transfer in a plane channel. ------> Wall velocity proportional to x //////////////////// North Wall /////////////////////// ------------------------------------------------------- parabolic -> outlet resistance mass-flux --> varies with velocity -> y value and constant ---------------------------------------------------- enthalpy //////////////////// South Wall /////////////////////// ^ y ------> Wall velocity proportional to x |-------> x-direction The inlet flow has a parabolic profile The outlet coefficient = A + B * sin[ pi *( y/yv]ast - 1/2 )] where A and B are constants, yvlast is channel width. North and south walls moved with velocity proportionally to longitudinal distance, x Their enthalpy is a piece-wise-linear function of x. South-wall boundary condition uses PATCH of SWALL type. North wall boundary condition uses PATCH of NORTH type and therefore need LAMWALL option. Solution should be symmetrical about channel centre plane. The Q1 contains PHOTON USE commands ENDDIS PHOTON USE p phi 1 8 1 msg picture is enlarged 8 times in y direction gr ou z 1 msg contours of H1 with Pr = 0.7 con h1 z 1 fi;0.001;upause 2 msg velocity vectors vec z 1 enduse REAL(VELMAX,UWALLMAX,TWALLMAX,ACON,BCON,PI) VELMAX=0.1; UWALLMAX=0.1; TWALLMAX=1.0; ACON=1.0; BCON=0.5 PI=3.14159 GROUP 3. X-direction grid specification GRDPWR(X,20,0.2,1.0) GROUP 4. Y-direction grid specification GRDPWR(Y,20,0.01,1.0) GROUP 7. Variables stored, solved & named SOLVE(P1,V1,U1,H1) GROUP 8. Terms (in differential equations) & devices ** De-activate the built-in source term TERMS(H1,N,Y,Y,Y,Y,Y) GROUP 9. Properties of the medium (or media) ENUL=1.545E-5; PRNDTL(H1)=0.7 ;RHO1=1.189 ! air at room temperature GROUP 13. Boundary conditions and special sources Inlet PATCH(INL,WEST,1,1,1,NY,1,1,1,1) ! parabolic velocity, uniform H1 (SOURCE of P1 at INL is RHO1*VELMAX*(1.-(2.*YG/YVLAST-1)^2)) (SOURCE of U1 at INL is VELMAX*(1.-(2.*YG/YVLAST-1)^2) with ONLYMS) (SOURCE of V1 at INL is 0.0 with ONLYMS) (SOURCE of H1 at INL is TWALLMAX*0.2 with ONLYMS) Walls CHAR(FORM) FORM=pwl3(xg,0,0,.4*xulast,-.5*:twallmax:,.5*xulast,$ 2.0*:twallmax:,xulast,:twallmax:) ! piece-wise linear formula for H1 PATCH(SW,SWALL,1,NX,1,1,1,NZ,1,1) (SOURCE of U1 at SW is UWALLMAX*XU/XULAST) (SOURCE of H1 at SW is :FORM:) PATCH(NW,NORTH,1,NX,NY,NY,1,NZ,1,1) (SOURCE of U1 at NW is UWALLMAX*XU/XULAST with LAMWALL) (SOURCE of H1 at NW is :FORM: with LAMWALL) Outlet PATCH(OUT,EAST,NX,NX,1,NY,1,NZ,1,1) FORM=-P1*(:ACON:+:BCON:*COS(:PI:*(YG/:YVLAST:-0.5))) (SOURCE of P1 at OUT is :FORM: with LINE) GROUP 15. Termination of sweeps LSWEEP=50 SELREF=T; RESFAC=0.1 GROUP 17. Under-relaxation devices RELAX(U1,FALSDT,0.1); RELAX(V1,FALSDT,0.1) GROUP 22. Spot-value print-out IYMON=19; IXMON=10; TSTSWP=-1 #conprom #maxabs #endpause distil=t