TALK=T;RUN( 1, 3)
TEXT(Domain-partitioning technique.
TITLE
DISPLAY
This example illustrates the Domain-partitioning technique.
This steady case simulate simple flow in xy plane. On the west
side there is an
inlet
and the north and east sides are
outlets
.
The south side is a
wall
with constant temperature.
In the west-south corner of the domain area there is a
source of a pollution concentration.
The solved area is divided into two parts.
Outlet
---------------------
I ! ! ! O
U1 n ! ! ! u
--> l ! 1st run ! 2nd run ! t
e ! ! ! l
Y ! t !C ! ! e
! --------------------- t
! Wall
!----- X
The solution in each part is made in a separate run.
The TROB$
1 transfer object on the east boundary is formed
at the end of the first run by means of an '(export' In-Form stat$
ement.
It stores the values of the outlet mass flux and the dependent
variables values on this boundary to transfer them to the second $
run.
The second run reads the information at the west boundary from the
TROB1 ob$
ject by means of an '(import' In-Form statement.
The third run simulates the flow without the domain-partitioning
technique for comparison with the results of the previous runs.
The Q1 contains PHOTON USE commands.
ENDDIS
PHOTON USE
p
phi1
gr z 1
set ve ref 0.2
ve z 1 sh
msg velocity vectors of first run
pause
ve cl;red
msg pressure contours of first run
con P1 z 1 fi;0.0001
pause
con cl;red
set con scale range on
con h1 z 1 fi;0 1;0.0001
set con scale range off
msg temperature contours of first run
pause
con cl;red
set con scale range on
con c1 z 1 fi;0 1;0.0001
set con scale range off
msg comcentration contours of first run
pause
p
phi2
gr z 1
set ve ref 0.2
ve z 1 sh
msg velocity vectors of second run
pause
ve cl;red
msg pressure contours of second run
con P1 z 1 fi;0.0001
pause
con cl;red
set con scale range on
con h1 z 1 fi;0 1;0.0001
set con scale range off
msg temperature contours of second run
pause
con cl;red
set con scale range on
con c1 z 1 fi;0 1;0.0001
set con scale range off
msg comcentration contours of second run
pause
p
phi3
gr z 1
set ve ref 0.2
ve z 1 sh
msg velocity vectors of third run
pause
ve cl;red
msg pressure contours of third run
con P1 z 1 fi;0.0001
pause
con cl;red
set con scale range on
con h1 z 1 fi;0 1;0.0001
*set con scale range off
msg temperature contours of third run
pause
con cl;red
set con scale range on
con c1 z 1 fi;0 1;0.0001
set con scale range off
msg comcentration contours of third run
enduse
First run
---------
GROUP 1. Run title and other preliminaries
TEXT(domain-partitioning technique. First run
Problem data
REAL(UINL,VINL,HINL,UMOV,HMOV,POUT,PCOF)
UINL=0.1 ! Inlet X velocity
VINL=0.0 ! Inlet Y velocity
HINL=1.0 ! Inlet enthalpy
UMOV=0.0 ! Moving wall velocity
HMOV=0.0 ! Moving wall enthalpy
POUT=0.0 ! Outlet pressure
PCOF=1.E+3 ! Outlet pressure coefficient
GROUP 3. x-direction grid specification
GRDPWR(X,20,0.2,1.0)
GROUP 4. Y-direction grid specification
GRDPWR(Y,20,0.2,1.0)
GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1,H1,C1)
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)
RHO1=1.189
ENUL=1.E-5; PRNDTL(H1)=0.7
GROUP 11. Initialization of variable or porosity fields
FIINIT(U1)=UINL
GROUP 13. Boundary conditions and special sources
! Inlet sources
PATCH(IN,West,1,1,1,NY,1,NZ,1,1)
COVAL(IN,P1,FIXFLU,RHO1*UINL)
COVAL(IN,U1,ONLYMS,UINL)
COVAL(IN,V1,ONLYMS,VINL)
COVAL(IN,H1,ONLYMS,HINL)
! Outlets sources
PATCH(OUT,EAST,NX,NX,1,NY,1,NZ,1,1)
COVAL(OUT,P1,PCOF,POUT)
PATCH(NOUT,NORTH,1,NX,NY,NY,1,NZ,1,1)
COVAL(NOUT,P1,PCOF,POUT)
COVAL(NOUT,H1,ONLYMS,SAME)
! Wall sources
PATCH(SW,SWALL,1,NX,1,1,1,NZ,1,1)
COVAL(SW,U1,1.,UMOV)
COVAL(SW,H1,1.,HMOV)
! source of a pollution concentration
PATCH(SOURCE,CELL,1,1,1,1,1,NZ,1,1)
COVAL(SOURCE,C1,FIXFLU,0.001)
! creation export transfer object
PATCH(PAT1,EAST,NX,NX,1,NY,1,NZ,1,1)
(EXPORT in TROB1 at PAT1)
GROUP 15. Termination of sweeps
LSWEEP=800
Group 18. Limits
VARMIN(C1)=0.
GROUP 22. Spot-value print-out
IYMON=NY/2; IZMON=10; TSTSWP=-1
distil=t
EX(P1)=1.364E-04
EX(U1)=9.971E-02
EX(V1)=1.457E-03
EX(H1)=9.800E-01
EX(C1)=4.900E-02
NSAVE=PHI1
LSG57=F
NXPRIN=1;NYPRIN=1
STOP
Second run
---------
GROUP 1. Run title and other preliminaries
TEXT(domain-partitioning technique. Second run
Problem data
REAL(UINL,VINL,HINL,UMOV,HMOV,POUT,PCOF)
UINL=0.1 ! Inlet X velocity
VINL=0.0 ! Inlet Y velocity
HINL=1.0 ! Inlet enthalpy
UMOV=0.0 ! Moving wall velocity
HMOV=0.0 ! Moving wall enthalpy
POUT=0.0 ! Outlet pressure
PCOF=1.E+3 ! Outlet pressure coefficient
GROUP 3. x-direction grid specification
GRDPWR(X,20,0.2,1.0)
GROUP 4. Y-direction grid specification
GRDPWR(Y,20,0.2,1.0)
GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1,H1,C1)
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)
RHO1=1.189
ENUL=1.E-5; PRNDTL(H1)=0.7
GROUP 11. Initialization of variable or porosity fields
FIINIT(U1)=UINL
GROUP 13. Boundary conditions and special sources
! Reading import transfer object
PATCH(PAT2,WEST,1,1,1,NY,1,NZ,1,1)
(IMPORT from TROB1 at PAT2)
! Outlets sources
PATCH(OUT,EAST,NX,NX,1,NY,1,NZ,1,1)
COVAL(OUT,P1,PCOF,POUT)
PATCH(NOUT,NORTH,1,NX,NY,NY,1,NZ,1,1)
COVAL(NOUT,P1,PCOF,POUT)
COVAL(NOUT,H1,ONLYMS,SAME)
! Wall sources
PATCH(SW,SWALL,1,NX,1,1,1,NZ,1,1)
COVAL(SW,U1,1.,UMOV)
COVAL(SW,H1,1.,HMOV)
GROUP 15. Termination of sweeps
LSWEEP=800
Group 18. Limits
VARMIN(C1)=0.
GROUP 22. Spot-value print-out
IYMON=NY/2; IZMON=10; TSTSWP=-1
distil=t
EX(P1)=7.119E-05
EX(U1)=9.724E-02
EX(V1)=2.698E-03
EX(H1)=9.438E-01
EX(C1)=6.385E-02
NSAVE=PHI2
LSG57=F
NXPRIN=1;NYPRIN=1
STOP
Third run
---------
GROUP 1. Run title and other preliminaries
TEXT(standard calculation. Third run
Problem data
REAL(UINL,VINL,HINL,UMOV,HMOV,POUT,PCOF)
UINL=0.1 ! Inlet X velocity
VINL=0.0 ! Inlet Y velocity
HINL=1.0 ! Inlet enthalpy
UMOV=0.0 ! Moving wall velocity
HMOV=0.0 ! Moving wall enthalpy
POUT=0.0 ! Outlet pressure
PCOF=1.E+3 ! Outlet pressure coefficient
GROUP 3. x-direction grid specification
GRDPWR(X,40,0.4,1.0)
GROUP 4. Y-direction grid specification
GRDPWR(Y,20,0.2,1.0)
GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1,H1,C1)
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)
RHO1=1.189
ENUL=1.E-5; PRNDTL(H1)=0.7
GROUP 11. Initialization of variable or porosity fields
FIINIT(U1)=UINL
GROUP 13. Boundary conditions and special sources
! Inlet sources
PATCH(IN,West,1,1,1,NY,1,NZ,1,1)
COVAL(IN,P1,FIXFLU,RHO1*UINL)
COVAL(IN,U1,ONLYMS,UINL)
COVAL(IN,V1,ONLYMS,VINL)
COVAL(IN,H1,ONLYMS,HINL)
! Outlets sources
PATCH(OUT,EAST,NX,NX,1,NY,1,NZ,1,1)
COVAL(OUT,P1,PCOF,POUT)
PATCH(NOUT,NORTH,1,NX,NY,NY,1,NZ,1,1)
COVAL(NOUT,P1,PCOF,POUT)
COVAL(NOUT,H1,ONLYMS,SAME)
! Wall sources
PATCH(SW,SWALL,1,NX,1,1,1,NZ,1,1)
COVAL(SW,U1,1.,UMOV)
COVAL(SW,H1,1.,HMOV)
! source of a pollution concentration
PATCH(SOURCE,CELL,1,1,1,1,1,NZ,1,1)
COVAL(SOURCE,C1,FIXFLU,0.001)
GROUP 15. Termination of sweeps
LSWEEP=800
Group 18. Limits
VARMIN(C1)=0.
GROUP 22. Spot-value print-out
IYMON=NY/2; IZMON=10; TSTSWP=-1
distil=t
EX(P1)=1.016E-04
EX(U1)=9.865E-02
EX(V1)=2.119E-03
EX(H1)=9.694E-01
EX(C1)=5.740E-02
NSAVE=PHI3
LSG57=F
NXPRIN=1;NYPRIN=1
STOP
STOP