CFD Software From ACRi -- CFDStudio, ANSWER, PORFLOW, TIDAL and RADM

************************************************************************************************
TITLe CASE 1.01 - Transient Pressure Decay to Zero
************************************************************************
/
GRID NODEs 5 BY 5 BY 5
/
DENSity = 1
FLUId thermal CONDuctivity = 1.
ROCK DENSity 1, POROsity = 1.
HYDRAULIC PROPS SS = 1, KX = 1., KY =1., KZ=1.
THERMAL PROPS CP = 1, KT = 1., DISL = 0., DIST = 0.
TRANSPORT KD = 0, DF = 1., DISL = 0., DIST = 0.
/
LOCAte (2,2,2) to (4,4,4) FIELD only
SET P = 1 for SELEcted subdomain
SET T = 1 for SELEcted subdomain
SET C = 1 for SELEcted subdomain
/
DIAGNOSTIC NODE AT 3,3,3 PRINT EVERY 10 STEPS
DEBUG GEOMETRY OFF
/
SOLVE UNTIL TIME = 0.20 IN STEPS OF 0.01
SELEct window from 1,1,3 to 5,5,3
OUTPut U,V,W,P,T,C for SELEcted window
SAVE OFF
/
END
/
/

-Back To Top

************************************************************************
TITLe CASE 1.02 - Uniform Velocity with Pressure + Normal Grad B.C.
************************************************************************
/
GRID NODEs 5 BY 5
/
SET Linear Function: P = 0. -1 X
/
/ Specify boundary conditions at the north and south boundaries
/
BOUNDARY CONDITIONS FOR P: at Y- FLUX = 0.
BOUNDARY CONDITIONS FOR P: at Y+ FLUX = 0.
/
HYDRAULIC PROPERTIES SS = 1., KX=1, KY=1
DIAGNOSTIC NODE AT 3,3 PRINT EVERY 10 STEPS
DEBUG GEOMETRY OFF
SOLVE UNTIL TIME = 1.5 IN STEPS OF 0.03
/
OUTPut variables U and P
SAVE OFF
END
/
/

-Back To Top

************************************************************************
TITLe CASE 1.03 - Linear U + Linear T with Heat Source
************************************************************************
/
GRID NODEs 5 by 5 BY 5
/
*************>>>> BOUNDARY CONDITIONS
/
BOUNdary X+ P as LINEar function: 0. -1. * Y in SELEcted subregion
BOUNdary Y+ P is LINEar function: 0. -1. * X in SELEcted subregion
BOUNdary conditions for P at Z- FLUX = 0.
BOUNdary conditions for P at Z+ FLUX = 0.
BOUNdary conditions for T at X+ value= 1.
BOUNdary conditions for T at Y- GRAD = 0.
BOUNdary conditions for T at Y+ GRAD = 0.
BOUNdary conditions for T at Z- FLUX = 0.
BOUNdary conditions for T at Z+ FLUX = 0.
/
FOR ZONE 1
DENSity RHO=1
FLUId thermal CONDuctivity = 1.
HYDRAULIC PROPS SS=1., KX=1., KY=1., KZ=1.
ROCK DENSity 1, POROsity 1
THERMAL PROPERTIES: CP=1., KT=1 ALFAL=0. ALFAT=0.
/
SOURCE FOR T: STRENGTH = 0.1
/
DIAGNOSTIC NODE AT 3,3,3 PRINT EVERY 20 STEPS
DEBUG GEOMETRY OFF
SELEct from (1,1,3) to (5,5,3)
OUTPut U,V,W,P,T for SELEcted window
/
SOLVE UNTIL TIME = 2.0 IN STEPS OF 0.01
SAVE OFF
END
/
/

-Back To Top

************************************************************************
TITLe CASE 1.04 - Rectangular Rod with Fixed End Temperatures
************************************************************************
/
GRID NODEs 5 BY 4
/
COORdinate X values at NODEs are: 0.4, 0.5, 0.6, 0.8, 1.0
/
BOUNDARY CONDITIONS FOR T at X- VAL = 0.4
BOUNDARY CONDITIONS FOR T at X+ VAL = 1.0
BOUNDARY CONDITIONS FOR T at Y- GRAD = 0.
BOUNDARY CONDITIONS FOR T at Y+ GRAD = 0.
/
DENSity = 1
FLUId thermal CONDuctivity = 1.
THERMAL PROPS CP = 1, KT = 1
/
DIAGNOSTIC NODE AT 3,3 PRINT EVERY 20 STEPS
DEBUG GEOMETRY OFF
SOLVE UNTIL TIME = 0.4 IN STEPS OF 0.005
SAVE OFF
OUTPut T
END
/
/

-Back To Top

************************************************************************
TITLe CASE 1.05 - Cylindrical Ring with Fixed End T and C
************************************************************************
/
GRID NODEs 3 BY 5
COORdinate R values at NODEs are: 0.4, 0.5, 0.6, 0.8, 1.0
$
BOUNDARY CONDITIONS FOR T Y- VAL=.0837
BOUNDARY CONDITIONS FOR T Y+ VAL=1.0
$
BOUNDARY CONDITIONS FOR C Y- VAL=.4
BOUNDARY CONDITIONS FOR C Y+ VAL=1.0
$
DENSity 1
FLUId thermal CONDuctivity = 1.
ROCK DENSity = 1, POROsity =1
THERMAL PROPS CP=1, KT=1
TRANSPORT KD = 0., DFM = 1.
/
DIAGNOSTIC TIME T C NODE AT 2,3 PRINT EVERY 20 STEPS
DEBUG GEOMETRY OFF
OUTPut variables T and C
/
SOLVE UNTIL TIME = 0.5 IN STEPS OF .005
SAVE OFF
END
/
/

-Back To Top

************************************************************************
TITLe CASE 1.06 - End Temperature Fixed Uniform Velocity Field
************************************************************************
/
GRID NODEs 5,3
$
BOUNDARY CONDITIONS FOR T X- VALUE = 1.
BOUNDARY CONDITIONS FOR T X+ VALUE = 0.
$
SET U 1.
DENSity 1
FLUId thermal CONDuctivity = 1.
ROCK DENSity 1, POROsity = 1
THERMAL PROPS CP=1, KT=1
DIAGNOSTIC NODE AT 3,2 PRINT EVERY 10 STEPS
DEBUG GEOMETRY OFF
SOLVE UNTIL TIME =.4 IN STEPS OF .01
SAVE OFF
OUTPut variable T
END
/

-Back To Top

************************************************************************
TITLe CASE 1.07 - Comparison of ADI Results with SOR (Next Set)
************************************************************************
/
GRID NODEs 5 BY 5 BY 5
COORdinate X RANGe 4
COORdinate Y RANGe 4
COORdinate Z RANGe 4
/X DIRECTION
BOUNDARY FOR P BOUNDARY at X-: VALUE = 10
BOUNDARY FOR P BOUNDARY at X+: VALUE = 0
/Y DIRECTION
BOUNDARY FOR P BOUNDARY at Y-: VALUE = 10
BOUNDARY FOR P BOUNDARY at Y+: VALUE = 10
/Z DIRECTION
BOUNDARY FOR P BOUNDARY at Z-: VALUE = 10
BOUNDARY FOR P BOUNDARY at Z+: VALUE = 0
/
MATRIX FOR P BY ADI
DENSity = 1
ROCK DENSity 1, POROSITY = 1.
HYDRAULIC PROPS SS = 1, KX = 1., KY =1., KZ=1.
/
SELEct window (1,1,2) to (100,100,2)
OUTPut U,P in XY format for SELEcted window
DIAGNOSTIC NODE AT 3,3,3 PRINT EVERY 10 STEPS
DEBUG GEOMETRY OFF
/
SOLVE FOR TIME = 1.0 TIME STEP 0.1
SAVE OFF
END
/
/

-Back To Top

************************************************************************
TITLe CASE 1.08 - Comparison of SOR Results with ADI (Previous Set)
************************************************************************
/
GRID NODEs 5 BY 5 BY 5
COORdinate X RANGe 4
COORdinate Y RANGe 4
COORdinate Z RANGe 4
/X DIRECTION
BOUNDARY FOR P BOUNDARY at X-: VALUE =10
BOUNDARY FOR P BOUNDARY at X+: VALUE =0
/Y DIRECTION
BOUNDARY FOR P BOUNDARY at Y-: VALUE =10
BOUNDARY FOR P BOUNDARY at Y+: VALUE =10
/Z DIRECTION
BOUNDARY FOR P BOUNDARY at Z-: VALUE =10
BOUNDARY FOR P BOUNDARY at Z+: VALUE =0
/
MATRIX FOR P BY SOR
DENSity = 1
ROCK DENSity 1, POROSITY = 1.
HYDRAULIC PROPS SS = 1, KX = 1., KY =1., KZ=1.
/
SELEct window (1,1,2) to (100,100,2)
OUTPut U,P in XY format for SELEcted window
DIAGNOSTIC NODE AT 3,3,3 PRINT EVERY 10 STEPS
DEBUG GEOMETRY OFF
/
SOLVE FOR TIME = 1.0 TIME STEP 0.1
SAVE OFF
END
/
/

-Back To Top

************************************************************************
TITLe CASE 1.09 - Bouyancy Affected Flow in a HLW Repository
************************************************************************
**** INEL BENCHMARK TEST CASE BT4
**** Magnuson S.W., R.G. Baca & A.J. Sondrup, August 1990.
**** Independent Verification and Benchmark Testing of
**** The PORFLO-3 Computer Code, Version 1.0
**** EGG-BG-9175; INEL, Idaho Falls, ID 83415
************************************************************************
/
USER AJ SONDRUP $ last revision 7/10/90
/
GRID NODEs 65 by 69
/
COORdinate X NODE values are:
0 50 100 150 200 240 280 320 350 374 394 410 422 432 438
442 445 447 449 451 453 456 460 464 468 472 476 480 484 488
492 496 500 504 508 512 516 520 524 528 532 536 540 544 547
549 551 553 555 558 562 568 578 590 606 626 650 680 720 760
800 850 900 950 1000
COORdinate Y NODE values are:
-1000 -950 -900 -850 -800 -760 -730 -710 -690 -675 -660 -650
-640 -630 -620 -610 -600 -590 -580 -570 -562 -554 -548 -544
-540 -536 -532 -528 -524 -520 -516 -512
-508 -504 -500 -496 -492
-488 -484 -480 -476 -472 -468 -464 -460
-456 -452 -446 -438 -430 -420 -410 -400 -390 -380 -370 -360
-350 -340 -325 -310 -290 -270 -240 -200 -150 -100 -50 0
/
*****BOUNDARY CONDITIONS*********************************************
BOUNdary T: X+ flux = 0.
BOUNdary T: X- flux = 0.
BOUNdary T: Y+ value = 30.
BOUNdary T: Y- value = 30.
/
BOUNdary P: X+ flux = 0.
BOUNdary P: X- flux = 0.
/
*****INITIAL CONDITIONS**********************************************
SET T 30.
SET H 10.
/
*****FLUID PROPERTIES************************************************
/ Remember to change fluid properties if you change reference temperature
/
GRAV 0., -9.81
REFErence T = 30. base=273.15
DENSity 995.7 POWEr law a=0.2 tcrit=374.15
VISCosity 0. EXPOnential law b=1436.0
FLUId SPECific heat = 4184.0
FLUId thermal CONDuctivity = 2.0E+7
/
/*****SOLID PROPERTIES************************************************
HYDRaulic ss=1.0e-04, kx=3.0e-02 ky=3.0e-02 kz=3.0e-02
ROCK density=2400.0, porosities 0.01, 0.01, 0.01
THERmal sh=930.0, cond=5.6e+07, alphal=50.0, alphat=2.0
/
/*****SOURCE**********************************************************
LOCAte subregion from (20,33) to (46 37)
SOURce for T = 2.E+07 per unit VOLUme in SELected zone
/
/*****OPERATIONAL CONTROL*********************************************
CONVergence with REFErence to T: local mode; value = 1.E-5
MATRIX for T 3
/
DIAGnostic node (33,35) print every 10 steps
DEBUG GEOMETRY OFF
SELEct window from (11,11) to (60,60) interval (3,3)
OUTPut for SELEcted window
/
SOLVe for 5.0 years in steps of 0.25
!!!<in> OUTPut P and T NOW for SELEcted window
/
SOLVe for 20.0 years in steps of 0.25
!!!<in> OUTPut NOW for SELEcted window
/
!!!<in> SOLVe for 25.0 years in step of 0.25
/
save off
/
END

-Back To Top

1 2 3 4 5 6 7 8 9 10

# Back to PORFLOW EXAMPLES DATABASE TOC