Skip to content

Commit

Permalink
FDS Verification: add obst_mass_flux series
Browse files Browse the repository at this point in the history
  • Loading branch information
rmcdermo committed Jul 25, 2018
1 parent b5ea642 commit 9f08a06
Show file tree
Hide file tree
Showing 5 changed files with 171 additions and 0 deletions.
48 changes: 48 additions & 0 deletions Verification/Miscellaneous/obst_box_mass_flux.fds
@@ -0,0 +1,48 @@
&HEAD CHID='obst_box_mass_flux', TITLE='Test of MULT/OBST box exact mass flux'/

&MESH IJK=50,50,50, XB=-0.125,0.125,-0.125,0.125,-0.125,0.125/

&TIME T_END=0.02, RESTRICT_TIME_STEP=F/

&SPEC ID='HELIUM'/

&SURF ID='x1', MASS_FLUX(1)=1, SPEC_ID(1)='HELIUM', COLOR='GREEN', RAMP_MF(1)='he ramp'/ 0.01 kg/s
&SURF ID='x2', MASS_FLUX(1)=2, SPEC_ID(1)='HELIUM', COLOR='BLUE', RAMP_MF(1)='he ramp'/ 0.02 kg/s
&SURF ID='y1', MASS_FLUX(1)=3, SPEC_ID(1)='HELIUM', COLOR='PINK', RAMP_MF(1)='he ramp'/ 0.03 kg/s
&SURF ID='y2', MASS_FLUX(1)=1, SPEC_ID(1)='HELIUM', COLOR='YELLOW', RAMP_MF(1)='he ramp'/ 0.01 kg/s
&SURF ID='z1', MASS_FLUX(1)=2, SPEC_ID(1)='HELIUM', COLOR='ORANGE', RAMP_MF(1)='he ramp'/ 0.02 kg/s
&SURF ID='z2', MASS_FLUX(1)=3, SPEC_ID(1)='HELIUM', COLOR='RED', RAMP_MF(1)='he ramp'/ 0.03 kg/s

&RAMP ID='he ramp', T=0, F=0/
&RAMP ID='he ramp', T=0.01, F=1/
&RAMP ID='he ramp', T=0.02, F=1//

&MULT ID='cube array', DX=0.005,DY=0.005,DZ=0.005, I_UPPER=41,J_UPPER=41,K_UPPER=41/

&OBST XB=-0.105,-0.100,-0.105,-0.100,-0.105,-0.100
MULT_ID='cube array'
SURF_ID6='x1','x2','y1','y2','z1','z2'
SHAPE='BOX'
LENGTH=0.1
WIDTH =0.1
HEIGHT=0.1
XYZ=0,0,0 /

&VENT MB='XMIN', SURF_ID='OPEN'/
&VENT MB='XMAX', SURF_ID='OPEN'/
&VENT MB='YMIN', SURF_ID='OPEN'/
&VENT MB='YMAX', SURF_ID='OPEN'/
&VENT MB='ZMIN', SURF_ID='OPEN'/
&VENT MB='ZMAX', SURF_ID='OPEN'/

&DEVC XB=-0.125,0.125,-0.125,0.125,-0.125,0.125, QUANTITY='MASS FLUX', SPEC_ID='HELIUM', STATISTICS='SURFACE INTEGRAL', SURF_ID='x1', ID='MDOT HE X1'/
&DEVC XB=-0.125,0.125,-0.125,0.125,-0.125,0.125, QUANTITY='MASS FLUX', SPEC_ID='HELIUM', STATISTICS='SURFACE INTEGRAL', SURF_ID='x2', ID='MDOT HE X2'/
&DEVC XB=-0.125,0.125,-0.125,0.125,-0.125,0.125, QUANTITY='MASS FLUX', SPEC_ID='HELIUM', STATISTICS='SURFACE INTEGRAL', SURF_ID='y1', ID='MDOT HE Y1'/
&DEVC XB=-0.125,0.125,-0.125,0.125,-0.125,0.125, QUANTITY='MASS FLUX', SPEC_ID='HELIUM', STATISTICS='SURFACE INTEGRAL', SURF_ID='y2', ID='MDOT HE Y2'/
&DEVC XB=-0.125,0.125,-0.125,0.125,-0.125,0.125, QUANTITY='MASS FLUX', SPEC_ID='HELIUM', STATISTICS='SURFACE INTEGRAL', SURF_ID='z1', ID='MDOT HE Z1'/
&DEVC XB=-0.125,0.125,-0.125,0.125,-0.125,0.125, QUANTITY='MASS FLUX', SPEC_ID='HELIUM', STATISTICS='SURFACE INTEGRAL', SURF_ID='z2', ID='MDOT HE Z2'/

&SLCF PBY=0, QUANTITY='VELOCITY', VECTOR=T/
&SLCF PBY=0, QUANTITY='MASS FRACTION', SPEC_ID='HELIUM', CELL_CENTERED=T/

&TAIL /
40 changes: 40 additions & 0 deletions Verification/Miscellaneous/obst_cone_mass_flux.fds
@@ -0,0 +1,40 @@
&HEAD CHID='obst_cone_mass_flux', TITLE='Test of MULT/OBST cone exact mass flux'/

&MESH IJK=50,50,50, XB=-0.125,0.125,-0.125,0.125,-0.125,0.125/

&TIME T_END=0.02, RESTRICT_TIME_STEP=F/

&SPEC ID='HELIUM'/

&SURF ID='top', MASS_FLUX(1)=0.3088059488033187, SPEC_ID(1)='HELIUM', COLOR='GREEN', RAMP_MF(1)='he ramp'/ 0.01 kg/s
&SURF ID='bottom', MASS_FLUX(1)=2.5464790894703250, SPEC_ID(1)='HELIUM', COLOR='BLACK', RAMP_MF(1)='he ramp'/ 0.02 kg/s

&RAMP ID='he ramp', T=0, F=0/
&RAMP ID='he ramp', T=0.01, F=1/
&RAMP ID='he ramp', T=0.02, F=1/

&MULT ID='cube array', DX=0.005,DY=0.005,DZ=0.005, I_UPPER=41,J_UPPER=41,K_UPPER=41/

&OBST XB=-0.105,-0.100,-0.105,-0.100,-0.105,-0.100
MULT_ID='cube array'
SURF_IDS='top','top','bottom'
SHAPE='CONE'
RADIUS=0.05
HEIGHT=0.2
XYZ=0,0,-0.1
ORIENTATION=0,0,1/

&VENT MB='XMIN', SURF_ID='OPEN'/
&VENT MB='XMAX', SURF_ID='OPEN'/
&VENT MB='YMIN', SURF_ID='OPEN'/
&VENT MB='YMAX', SURF_ID='OPEN'/
&VENT MB='ZMIN', SURF_ID='OPEN'/
&VENT MB='ZMAX', SURF_ID='OPEN'/

&DEVC XB=-0.125,0.125,-0.125,0.125,-0.125,0.125, QUANTITY='MASS FLUX', SPEC_ID='HELIUM', STATISTICS='SURFACE INTEGRAL', SURF_ID='top', ID='MDOT HE TOP'/
&DEVC XB=-0.125,0.125,-0.125,0.125,-0.125,0.125, QUANTITY='MASS FLUX', SPEC_ID='HELIUM', STATISTICS='SURFACE INTEGRAL', SURF_ID='bottom', ID='MDOT HE BOTTOM'/

&SLCF PBY=0, QUANTITY='VELOCITY', VECTOR=T/
&SLCF PBY=0, QUANTITY='MASS FRACTION', SPEC_ID='HELIUM', CELL_CENTERED=T/

&TAIL /
42 changes: 42 additions & 0 deletions Verification/Miscellaneous/obst_cylinder_mass_flux.fds
@@ -0,0 +1,42 @@
&HEAD CHID='obst_cylinder_mass_flux', TITLE='Test of MULT/OBST cylinder exact mass flux'/

&MESH IJK=50,50,50, XB=-0.125,0.125,-0.125,0.125,-0.125,0.125/

&TIME T_END=0.02, RESTRICT_TIME_STEP=F/

&SPEC ID='HELIUM'/

&SURF ID='top', MASS_FLUX(1)=1.2732395447351626, SPEC_ID(1)='HELIUM', COLOR='GRAY', RAMP_MF(1)='he ramp'/ 0.01 kg/s
&SURF ID='side', MASS_FLUX(1)=0.3183098861837907, SPEC_ID(1)='HELIUM', COLOR='ORANGE', RAMP_MF(1)='he ramp'/ 0.02 kg/s
&SURF ID='bottom', MASS_FLUX(1)=3.8197186342054880, SPEC_ID(1)='HELIUM', COLOR='BLACK', RAMP_MF(1)='he ramp'/ 0.03 kg/s

&RAMP ID='he ramp', T=0, F=0/
&RAMP ID='he ramp', T=0.01, F=1/
&RAMP ID='he ramp', T=0.02, F=1/

&MULT ID='cube array', DX=0.005,DY=0.005,DZ=0.005, I_UPPER=41,J_UPPER=41,K_UPPER=41/

&OBST XB=-0.105,-0.100,-0.105,-0.100,-0.105,-0.100
MULT_ID='cube array'
SURF_IDS='top','side','bottom'
SHAPE='CYLINDER'
RADIUS=0.05
HEIGHT=0.2
XYZ=0,0,-0.1
ORIENTATION=0,0,1/

&VENT MB='XMIN', SURF_ID='OPEN'/
&VENT MB='XMAX', SURF_ID='OPEN'/
&VENT MB='YMIN', SURF_ID='OPEN'/
&VENT MB='YMAX', SURF_ID='OPEN'/
&VENT MB='ZMIN', SURF_ID='OPEN'/
&VENT MB='ZMAX', SURF_ID='OPEN'/

&DEVC XB=-0.125,0.125,-0.125,0.125,-0.125,0.125, QUANTITY='MASS FLUX', SPEC_ID='HELIUM', STATISTICS='SURFACE INTEGRAL', SURF_ID='top', ID='MDOT HE TOP'/
&DEVC XB=-0.125,0.125,-0.125,0.125,-0.125,0.125, QUANTITY='MASS FLUX', SPEC_ID='HELIUM', STATISTICS='SURFACE INTEGRAL', SURF_ID='side', ID='MDOT HE SIDE'/
&DEVC XB=-0.125,0.125,-0.125,0.125,-0.125,0.125, QUANTITY='MASS FLUX', SPEC_ID='HELIUM', STATISTICS='SURFACE INTEGRAL', SURF_ID='bottom', ID='MDOT HE BOTTOM'/

&SLCF PBY=0, QUANTITY='VELOCITY', VECTOR=T/
&SLCF PBY=0, QUANTITY='MASS FRACTION', SPEC_ID='HELIUM', CELL_CENTERED=T/

&TAIL /
5 changes: 5 additions & 0 deletions Verification/Miscellaneous/obst_mass_flux.csv
@@ -0,0 +1,5 @@
s,kg/s,kg/s,kg/s
Time,MDOT_1,MDOT_2,MDOT_3
0,0,0,0
0.01,0.01,0.02,0.03
0.02,0.01,0.02,0.03
36 changes: 36 additions & 0 deletions Verification/Miscellaneous/obst_sphere_mass_flux.fds
@@ -0,0 +1,36 @@
&HEAD CHID='obst_sphere_mass_flux', TITLE='Test of MULT/OBST sphere exact mass flux'/

&MESH IJK=50,50,50, XB=-0.125,0.125,-0.125,0.125,-0.125,0.125/

&TIME T_END=0.02, RESTRICT_TIME_STEP=F/

&SPEC ID='HELIUM'/

&SURF ID='shape', MASS_FLUX(1)=0.0795774715459477, SPEC_ID(1)='HELIUM', COLOR='ORANGE', RAMP_MF(1)='he ramp'/ 0.01 kg/s

&RAMP ID='he ramp', T=0, F=0/
&RAMP ID='he ramp', T=0.01, F=1/
&RAMP ID='he ramp', T=0.02, F=1/

&MULT ID='cube array', DX=0.005,DY=0.005,DZ=0.005, I_UPPER=41,J_UPPER=41,K_UPPER=41/

&OBST XB=-0.105,-0.100,-0.105,-0.100,-0.105,-0.100
MULT_ID='cube array'
SURF_ID='shape'
SHAPE='SPHERE'
RADIUS=0.1/

&VENT MB='XMIN', SURF_ID='OPEN'/
&VENT MB='XMAX', SURF_ID='OPEN'/
&VENT MB='YMIN', SURF_ID='OPEN'/
&VENT MB='YMAX', SURF_ID='OPEN'/
&VENT MB='ZMIN', SURF_ID='OPEN'/
&VENT MB='ZMAX', SURF_ID='OPEN'/

&DEVC XB=-0.125,0.125,-0.125,0.125,-0.125,0.125,
QUANTITY='MASS FLUX', SPEC_ID='HELIUM', STATISTICS='SURFACE INTEGRAL', SURF_ID='shape', ID='MDOT HE'/

&SLCF PBY=0, QUANTITY='VELOCITY', VECTOR=T/
&SLCF PBY=0, QUANTITY='MASS FRACTION', SPEC_ID='HELIUM', CELL_CENTERED=T/

&TAIL /

0 comments on commit 9f08a06

Please sign in to comment.