Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
26 changes: 21 additions & 5 deletions Source/wall.f90
Original file line number Diff line number Diff line change
Expand Up @@ -1818,7 +1818,7 @@ SUBROUTINE SOLID_HEAT_TRANSFER(NM,T,DT_BC,PARTICLE_INDEX,WALL_INDEX,CFACE_INDEX,
INTEGER, INTENT(IN), OPTIONAL:: WALL_INDEX,PARTICLE_INDEX,CFACE_INDEX,THIN_WALL_INDEX
REAL(EB) :: RDT,DTMP,QDXKF,QDXKB,RR,RFACF,RFACB,RFACF2,RFACB2,Q_RAD_IN_B,RFLUX_UP,RFLUX_DOWN,EMISSIVITY_BACK, &
VOLSUM,KAPSUM,DXF,DXB,HTCF,HTCB,Q_RAD_OUT,Q_RAD_OUT_OLD,Q_CON_F,Q_CON_B,&
Q_LIQUID_F,Q_LIQUID_B,LAYER_DIVIDE,TMP_GAS_BACK,GEOM_FACTOR,DT_BC_SUB_OLD,&
Q_LIQUID_F,Q_LIQUID_B,LAYER_DIVIDE,TMP_BACK,TMP_GAS_BACK,GEOM_FACTOR,DT_BC_SUB_OLD,&
DEL_DOT_Q_SC,Q_DOT_G_PP,Q_DOT_G_PP_NET,Q_DOT_O2_PP,Q_DOT_O2_PP_NET,R_SURF,U_SURF,V_SURF,W_SURF,T_BC_SUB,DT_BC_SUB,&
Q_NET_F,Q_NET_B,TMP_RATIO,KODXF,KODXB,H_S,T_NODE,C_S,H_NODE,VOL,&
RADIUS,HTC_LIMIT,CP1,CP2,DENOM,THICKNESS,DT_FO,DDSUM,NODE_RDT(NWP_MAX)
Expand All @@ -1837,7 +1837,7 @@ SUBROUTINE SOLID_HEAT_TRANSFER(NM,T,DT_BC,PARTICLE_INDEX,WALL_INDEX,CFACE_INDEX,
INTEGER :: NWP_NEW,I_GRAD,IZERO,SURF_INDEX,BACKING,NWP,I,NL,N,OBST_INDEX,&
N_CELLS,ITMP,ITER,BACK_MESH,BACK_INDEX,BACK_WALL_INDEX
CHARACTER(MESSAGE_LENGTH) :: MESSAGE
LOGICAL :: ISOLATED_THIN_WALL,ISOLATED_THIN_WALL_BACK,REMESH_LAYER(MAX_LAYERS),REMESH_CHECK,&
LOGICAL :: ISOLATED_THIN_WALL,ISOLATED_THIN_WALL_BACK,REMESH_LAYER(MAX_LAYERS),REMESH_CHECK,DIRICHLET_BACK,&
CELL_ZERO(MAX_LAYERS),TMP_CHECK(MAX_LAYERS)
TYPE(WALL_TYPE), POINTER :: WC,WC_BACK
TYPE(THIN_WALL_TYPE), POINTER :: TW,TW_BACK
Expand Down Expand Up @@ -2127,6 +2127,12 @@ SUBROUTINE SOLID_HEAT_TRANSFER(NM,T,DT_BC,PARTICLE_INDEX,WALL_INDEX,CFACE_INDEX,

CASE(VOID) ! Non-insulated backing to an ambient void

IF (SF%DIRICHLET_BACK) THEN
DIRICHLET_BACK = .TRUE.
TMP_BACK = SF%TMP_BACK
ELSE
DIRICHLET_BACK = .FALSE.
ENDIF
IF (SF%TMP_GAS_BACK>0._EB) THEN
TMP_GAS_BACK = TMP_0(BC%KK) + EVALUATE_RAMP(T-T_BEGIN,SF%RAMP(TIME_TGB)%INDEX)*(SF%TMP_GAS_BACK-TMP_0(BC%KK))
ELSE
Expand All @@ -2148,6 +2154,7 @@ SUBROUTINE SOLID_HEAT_TRANSFER(NM,T,DT_BC,PARTICLE_INDEX,WALL_INDEX,CFACE_INDEX,
EMISSIVITY_BACK = 0._EB
Q_LIQUID_B = 0._EB
TMP_GAS_BACK = B1%TMP_B
DIRICHLET_BACK = .FALSE.

CASE(EXPOSED) ! The backside is exposed to gas in current or adjacent mesh.

Expand All @@ -2160,6 +2167,12 @@ SUBROUTINE SOLID_HEAT_TRANSFER(NM,T,DT_BC,PARTICLE_INDEX,WALL_INDEX,CFACE_INDEX,
Q_LIQUID_B = 0._EB
ENDIF
Q_CON_B = HTCB*DTMP
IF (SF_BACK%DIRICHLET_FRONT) THEN
DIRICHLET_BACK = .TRUE.
TMP_BACK = SF_BACK%TMP_FRONT
ELSE
DIRICHLET_BACK = .FALSE.
ENDIF

END SELECT

Expand Down Expand Up @@ -2756,6 +2769,8 @@ SUBROUTINE SOLID_HEAT_TRANSFER(NM,T,DT_BC,PARTICLE_INDEX,WALL_INDEX,CFACE_INDEX,

IF (SF%DIRICHLET_FRONT) THEN
RFACF2 = -1._EB
B1%TMP_F = TMP_0(BC%KKG) + &
EVALUATE_RAMP(T-T_BEGIN,SF%RAMP(TIME_TEMP)%INDEX,TAU=SF%RAMP(TIME_TEMP)%TAU)*(SF%TMP_FRONT-TMP_0(BC%KKG))
QDXKF = 2._EB*B1%TMP_F
ELSEIF (.NOT.RADIATION .OR. ONE_D%INTERNAL_RADIATION .OR. ISOLATED_THIN_WALL) THEN
RFACF = 0.5_EB*HTCF
Expand All @@ -2767,8 +2782,9 @@ SUBROUTINE SOLID_HEAT_TRANSFER(NM,T,DT_BC,PARTICLE_INDEX,WALL_INDEX,CFACE_INDEX,
QDXKF = (HTCF*B1%TMP_G + Q_LIQUID_F + B1%Q_RAD_IN + 3._EB*B1%EMISSIVITY*SIGMA*B1%TMP_F**4) / (KODXF+RFACF)
ENDIF

IF (SF%DIRICHLET_BACK) THEN
IF (DIRICHLET_BACK) THEN
RFACB2 = -1._EB
B1%TMP_B = TMP_0(BC%KK) + EVALUATE_RAMP(T-T_BEGIN,SF%RAMP(TIME_TB)%INDEX)*(TMP_BACK-TMP_0(BC%KK))
QDXKB = 2._EB*B1%TMP_B
ELSEIF (.NOT.RADIATION .OR. ONE_D%INTERNAL_RADIATION .OR. ISOLATED_THIN_WALL_BACK) THEN
RFACB = 0.5_EB*HTCB
Expand Down Expand Up @@ -3251,17 +3267,17 @@ SUBROUTINE PYROLYSIS(N_MATS,MATL_INDEX,SURF_INDEX,IIG,JJG,KKG,TMP_S,TMP_F,Y_O2_F
IF (PRESENT(R_DROP)) THEN
LENGTH_SCALE = 2._EB*R_DROP
ELSE

LENGTH_SCALE = SF%CONV_LENGTH
ENDIF
RE_L = RHO_FILM*VEL*LENGTH_SCALE/MU_FILM
SELECT CASE(SF%GEOMETRY)
CASE DEFAULT ; SHERWOOD = 0.037_EB*SC_FILM**ONTH*RE_L**0.8_EB
CASE(SURF_SPHERICAL) ; SHERWOOD = 2._EB + 0.6_EB*SC_FILM**ONTH*SQRT(RE_L)
END SELECT

H_MASS = MAX(H_MASS_DNS,SHERWOOD*D_FILM/LENGTH_SCALE)

ENDIF

ENDIF H_MASS_IF

ENDIF IF_DO_EVAPORATION
Expand Down
7 changes: 2 additions & 5 deletions Verification/Heat_Transfer/ht1d_pile.fds
Original file line number Diff line number Diff line change
Expand Up @@ -9,11 +9,8 @@

&RADI RADIATION=F /

&SURF ID='HOT', TMP_FRONT=500, RAMP_T='temp', TMP_BACK=20., RAMP_TMP_BACK='temp', COLOR='RED', VARIABLE_THICKNESS=T /
&SURF ID='COLD', TMP_FRONT= 20, RAMP_T='temp', TMP_BACK=500., RAMP_TMP_BACK='temp', COLOR='BLUE', VARIABLE_THICKNESS=T /

&RAMP ID='temp', T=0, F=1 /
&RAMP ID='temp', T=1, F=1 /
&SURF ID='HOT', TMP_FRONT=500, TAU_T=0., COLOR='RED', VARIABLE_THICKNESS=T /
&SURF ID='COLD', TMP_FRONT= 20, TAU_T=0., COLOR='BLUE', VARIABLE_THICKNESS=T /

&VENT MB='XMIN', SURF_ID='OPEN' /
&VENT MB='XMAX', SURF_ID='OPEN' /
Expand Down
2 changes: 0 additions & 2 deletions Verification/Heat_Transfer/ht3d_ibeam.fds
Original file line number Diff line number Diff line change
Expand Up @@ -15,15 +15,13 @@

&SURF ID='ADIABATIC STEEL', COLOR='GRAY', HEAT_TRANSFER_COEFFICIENT=0, HT3D=T, CELL_SIZE=0.01 /
&SURF ID='HOT STEEL', COLOR='RED', TMP_FRONT=800, TAU_T=0., HT3D=T, CELL_SIZE=0.01 /
&SURF ID='HOT STEEL2', COLOR='GRAY', HEAT_TRANSFER_COEFFICIENT=0, TMP_BACK=800, HT3D=T, CELL_SIZE=0.01 /

&MATL ID = 'STEEL'
SPECIFIC_HEAT = 0.60
CONDUCTIVITY = 45.
DENSITY = 7850./

&VENT XB=.2,.2,0,.2,.02,.08, SURF_ID='HOT STEEL'/
&VENT XB=-.2,-.2,0,.2,.02,.08, SURF_ID='HOT STEEL2'/

&BNDF QUANTITY='WALL TEMPERATURE', CELL_CENTERED=T/
&BNDF QUANTITY='NET HEAT FLUX', CELL_CENTERED=T/
Expand Down
7 changes: 2 additions & 5 deletions Verification/Heat_Transfer/ht3d_pile.fds
Original file line number Diff line number Diff line change
Expand Up @@ -9,13 +9,10 @@

&RADI RADIATION=F /

&SURF ID='HOT', TMP_FRONT=500, RAMP_T='temp', TMP_BACK=20., RAMP_TMP_BACK='temp', COLOR='RED', HT3D=T, CELL_SIZE=0.01 /
&SURF ID='COLD', TMP_FRONT= 20, RAMP_T='temp', TMP_BACK=500., RAMP_TMP_BACK='temp', COLOR='BLUE', HT3D=T, CELL_SIZE=0.01 /
&SURF ID='HOT', TMP_FRONT=500, TAU_T=0., COLOR='RED', HT3D=T, CELL_SIZE=0.01 /
&SURF ID='COLD', TMP_FRONT= 20, TAU_T=0., COLOR='BLUE', HT3D=T, CELL_SIZE=0.01 /
&SURF ID='SIDE', HEAT_TRANSFER_COEFFICIENT=0, COLOR='GRAY', HT3D=T, CELL_SIZE=0.01 /

&RAMP ID='temp', T=0, F=1 /
&RAMP ID='temp', T=1, F=1 /

&VENT MB='XMIN', SURF_ID='OPEN' /
&VENT MB='XMAX', SURF_ID='OPEN' /
&VENT MB='YMIN', SURF_ID='OPEN' /
Expand Down
1 change: 0 additions & 1 deletion Verification/Heat_Transfer/ht3d_sphere_24.fds
Original file line number Diff line number Diff line change
Expand Up @@ -16,7 +16,6 @@
&SURF ID='SPHERE'
HT3D=T
TMP_FRONT=20.
TMP_BACK=20.
COLOR='ORANGE' /

&MULT ID='cube array', DX=0.01, DY=0.01, DZ=0.01, I_UPPER=21, J_UPPER=21, K_UPPER=21 /
Expand Down
1 change: 0 additions & 1 deletion Verification/Heat_Transfer/ht3d_sphere_48.fds
Original file line number Diff line number Diff line change
Expand Up @@ -23,7 +23,6 @@
&SURF ID='SPHERE'
HT3D=T
TMP_FRONT=20.
TMP_BACK=20.
COLOR='ORANGE' /

&MULT ID='cube array', DX=0.005, DY=0.005, DZ=0.005, I_UPPER=41, J_UPPER=41, K_UPPER=41 /
Expand Down
1 change: 0 additions & 1 deletion Verification/Heat_Transfer/ht3d_sphere_96.fds
Original file line number Diff line number Diff line change
Expand Up @@ -17,7 +17,6 @@
&SURF ID='SPHERE'
HT3D=T
TMP_FRONT=20.
TMP_BACK=20.
COLOR='ORANGE' /

&MULT ID='cube array', DX=0.0025, DY=0.0025, DZ=0.0025, I_UPPER=81, J_UPPER=81, K_UPPER=81 /
Expand Down
12 changes: 6 additions & 6 deletions Verification/Pyrolysis/part_baking_soda_420K.fds
Original file line number Diff line number Diff line change
Expand Up @@ -70,13 +70,13 @@
&INIT ID='part_init first-order', PART_ID='agent_part first-order', XB=-0.05,0.05,-0.05,0.05,0.95,1.05, N_PARTICLES = 1 /
&INIT ID='part_init spherical' , PART_ID='agent_part spherical', XB=-0.05,0.05,-0.05,0.05,0.95,1.05, N_PARTICLES = 1 /

&DEVC XYZ=0,0,1, QUANTITY='WALL THICKNESS', INIT_ID='part_init first-order', CONVERSION_FACTOR=2.E6, UNITS='mu-m', ID='diam first-order'/
&DEVC XYZ=0,0,1, QUANTITY='INSIDE WALL TEMPERATURE', INIT_ID='part_init first-order', DEPTH=1.E-9, ID='TMP first-order'/
&DEVC XYZ=0,0,1, QUANTITY='SOLID DENSITY', MATL_ID='Sodium Bicarbonate first-order', INIT_ID='part_init first-order', DEPTH=1.E-9, ID='RHO_S first-order'/
&DEVC QUANTITY='WALL THICKNESS', INIT_ID='part_init first-order', CONVERSION_FACTOR=2.E6, UNITS='mu-m', ID='diam first-order'/
&DEVC QUANTITY='INSIDE WALL TEMPERATURE', INIT_ID='part_init first-order', DEPTH=1.E-9, ID='TMP first-order'/
&DEVC QUANTITY='SOLID DENSITY', MATL_ID='Sodium Bicarbonate first-order', INIT_ID='part_init first-order', DEPTH=1.E-9, ID='RHO_S first-order'/

&DEVC XYZ=0,0,1, QUANTITY='WALL THICKNESS', INIT_ID='part_init spherical', CONVERSION_FACTOR=2.E6, UNITS='mu-m', ID='diam spherical'/
&DEVC XYZ=0,0,1, QUANTITY='INSIDE WALL TEMPERATURE', INIT_ID='part_init spherical', DEPTH=1.E-6, ID='TMP spherical'/
&DEVC XYZ=0,0,1, QUANTITY='SOLID DENSITY', MATL_ID='Sodium Bicarbonate spherical', INIT_ID='part_init spherical', DEPTH=1.E-9, ID='RHO_S spherical'/
&DEVC QUANTITY='WALL THICKNESS', INIT_ID='part_init spherical', CONVERSION_FACTOR=2.E6, UNITS='mu-m', ID='diam spherical'/
&DEVC QUANTITY='INSIDE WALL TEMPERATURE', INIT_ID='part_init spherical', DEPTH=1.E-6, ID='TMP spherical'/
&DEVC QUANTITY='SOLID DENSITY', MATL_ID='Sodium Bicarbonate spherical', INIT_ID='part_init spherical', DEPTH=1.E-9, ID='RHO_S spherical'/

&CTRL ID='MY POWER', FUNCTION_TYPE='POWER', INPUT_ID='RHO_S spherical','CONSTANT', CONSTANT=0.33333/
&DEVC ID='SPHERE DIA', XYZ=0,0,1, QUANTITY='CONTROL VALUE', CTRL_ID='MY POWER', CONVERSION_FACTOR=0.38444/ CONVERSION_FACTOR=DIAMETER/(RHO_S(0)^(1/3))
Expand Down
12 changes: 6 additions & 6 deletions Verification/Pyrolysis/part_baking_soda_450K.fds
Original file line number Diff line number Diff line change
Expand Up @@ -70,13 +70,13 @@
&INIT ID='part_init first-order', PART_ID='agent_part first-order', XB=-0.05,0.05,-0.05,0.05,0.95,1.05, N_PARTICLES = 1 /
&INIT ID='part_init spherical' , PART_ID='agent_part spherical', XB=-0.05,0.05,-0.05,0.05,0.95,1.05, N_PARTICLES = 1 /

&DEVC XYZ=0,0,1, QUANTITY='WALL THICKNESS', INIT_ID='part_init first-order', CONVERSION_FACTOR=2.E6, UNITS='mu-m', ID='diam first-order'/
&DEVC XYZ=0,0,1, QUANTITY='INSIDE WALL TEMPERATURE', INIT_ID='part_init first-order', DEPTH=1.E-9, ID='TMP first-order'/
&DEVC XYZ=0,0,1, QUANTITY='SOLID DENSITY', MATL_ID='Sodium Bicarbonate first-order', INIT_ID='part_init first-order', DEPTH=1.E-9, ID='RHO_S first-order'/
&DEVC QUANTITY='WALL THICKNESS', INIT_ID='part_init first-order', CONVERSION_FACTOR=2.E6, UNITS='mu-m', ID='diam first-order'/
&DEVC QUANTITY='INSIDE WALL TEMPERATURE', INIT_ID='part_init first-order', DEPTH=1.E-9, ID='TMP first-order'/
&DEVC QUANTITY='SOLID DENSITY', MATL_ID='Sodium Bicarbonate first-order', INIT_ID='part_init first-order', DEPTH=1.E-9, ID='RHO_S first-order'/

&DEVC XYZ=0,0,1, QUANTITY='WALL THICKNESS', INIT_ID='part_init spherical', CONVERSION_FACTOR=2.E6, UNITS='mu-m', ID='diam spherical'/
&DEVC XYZ=0,0,1, QUANTITY='INSIDE WALL TEMPERATURE', INIT_ID='part_init spherical', DEPTH=1.E-6, ID='TMP spherical'/
&DEVC XYZ=0,0,1, QUANTITY='SOLID DENSITY', MATL_ID='Sodium Bicarbonate spherical', INIT_ID='part_init spherical', DEPTH=1.E-9, ID='RHO_S spherical'/
&DEVC QUANTITY='WALL THICKNESS', INIT_ID='part_init spherical', CONVERSION_FACTOR=2.E6, UNITS='mu-m', ID='diam spherical'/
&DEVC QUANTITY='INSIDE WALL TEMPERATURE', INIT_ID='part_init spherical', DEPTH=1.E-6, ID='TMP spherical'/
&DEVC QUANTITY='SOLID DENSITY', MATL_ID='Sodium Bicarbonate spherical', INIT_ID='part_init spherical', DEPTH=1.E-9, ID='RHO_S spherical'/

&CTRL ID='MY POWER', FUNCTION_TYPE='POWER', INPUT_ID='RHO_S spherical','CONSTANT', CONSTANT=0.33333/
&DEVC ID='SPHERE DIA', XYZ=0,0,1, QUANTITY='CONTROL VALUE', CTRL_ID='MY POWER', CONVERSION_FACTOR=0.38444/ CONVERSION_FACTOR=DIAMETER/(RHO_S(0)^(1/3))
Expand Down
12 changes: 6 additions & 6 deletions Verification/Pyrolysis/part_baking_soda_500K.fds
Original file line number Diff line number Diff line change
Expand Up @@ -70,13 +70,13 @@
&INIT ID='part_init first-order', PART_ID='agent_part first-order', XB=-0.05,0.05,-0.05,0.05,0.95,1.05, N_PARTICLES = 1 /
&INIT ID='part_init spherical' , PART_ID='agent_part spherical', XB=-0.05,0.05,-0.05,0.05,0.95,1.05, N_PARTICLES = 1 /

&DEVC XYZ=0,0,1, QUANTITY='WALL THICKNESS', INIT_ID='part_init first-order', CONVERSION_FACTOR=2.E6, UNITS='mu-m', ID='diam first-order'/
&DEVC XYZ=0,0,1, QUANTITY='INSIDE WALL TEMPERATURE', INIT_ID='part_init first-order', DEPTH=1.E-9, ID='TMP first-order'/
&DEVC XYZ=0,0,1, QUANTITY='SOLID DENSITY', MATL_ID='Sodium Bicarbonate first-order', INIT_ID='part_init first-order', DEPTH=1.E-9, ID='RHO_S first-order'/
&DEVC QUANTITY='WALL THICKNESS', INIT_ID='part_init first-order', CONVERSION_FACTOR=2.E6, UNITS='mu-m', ID='diam first-order'/
&DEVC QUANTITY='INSIDE WALL TEMPERATURE', INIT_ID='part_init first-order', DEPTH=1.E-9, ID='TMP first-order'/
&DEVC QUANTITY='SOLID DENSITY', MATL_ID='Sodium Bicarbonate first-order', INIT_ID='part_init first-order', DEPTH=1.E-9, ID='RHO_S first-order'/

&DEVC XYZ=0,0,1, QUANTITY='WALL THICKNESS', INIT_ID='part_init spherical', CONVERSION_FACTOR=2.E6, UNITS='mu-m', ID='diam spherical'/
&DEVC XYZ=0,0,1, QUANTITY='INSIDE WALL TEMPERATURE', INIT_ID='part_init spherical', DEPTH=1.E-6, ID='TMP spherical'/
&DEVC XYZ=0,0,1, QUANTITY='SOLID DENSITY', MATL_ID='Sodium Bicarbonate spherical', INIT_ID='part_init spherical', DEPTH=1.E-9, ID='RHO_S spherical'/
&DEVC QUANTITY='WALL THICKNESS', INIT_ID='part_init spherical', CONVERSION_FACTOR=2.E6, UNITS='mu-m', ID='diam spherical'/
&DEVC QUANTITY='INSIDE WALL TEMPERATURE', INIT_ID='part_init spherical', DEPTH=1.E-6, ID='TMP spherical'/
&DEVC QUANTITY='SOLID DENSITY', MATL_ID='Sodium Bicarbonate spherical', INIT_ID='part_init spherical', DEPTH=1.E-9, ID='RHO_S spherical'/

&CTRL ID='MY POWER', FUNCTION_TYPE='POWER', INPUT_ID='RHO_S spherical','CONSTANT', CONSTANT=0.33333/
&DEVC ID='SPHERE DIA', XYZ=0,0,1, QUANTITY='CONTROL VALUE', CTRL_ID='MY POWER', CONVERSION_FACTOR=0.38444/ CONVERSION_FACTOR=DIAMETER/(RHO_S(0)^(1/3))
Expand Down
2 changes: 1 addition & 1 deletion Verification/Pyrolysis/surf_mass_two_species_cart.fds
Original file line number Diff line number Diff line change
Expand Up @@ -58,7 +58,7 @@

&INIT ID = 'init_1', PART_ID='Wood_part', XB=0.,1.,0.,1.,0.,1., N_PARTICLES=1, MASS_PER_VOLUME = 0.1 /

&DEVC XB=0.,1.,0.,1.,0.,1., QUANTITY='WALL TEMPERATURE', INIT_ID='init_1', ID='Particle T', SPATIAL_STATISTIC='MAX'/
&DEVC QUANTITY='WALL TEMPERATURE', INIT_ID='init_1', ID='Particle T', SPATIAL_STATISTIC='MAX'/
&DEVC XB=0.,1.,0.,1.,0.,1., QUANTITY='DENSITY', ID='fuel gas mass', SPATIAL_STATISTIC='VOLUME INTEGRAL', SPEC_ID='FUEL' /
&DEVC XB=0.,1.,0.,1.,0.,1., QUANTITY='DENSITY', ID='water vapor mass', SPATIAL_STATISTIC='VOLUME INTEGRAL', SPEC_ID='WATER VAPOR' /
&DEVC XB=0.,1.,0.,1.,0.,1., QUANTITY='MPUV', PART_ID='Wood_part', ID='solid mass', SPATIAL_STATISTIC='VOLUME INTEGRAL' /
Expand Down
2 changes: 1 addition & 1 deletion Verification/Pyrolysis/surf_mass_two_species_cyl.fds
Original file line number Diff line number Diff line change
Expand Up @@ -56,7 +56,7 @@

&INIT ID = 'init_1', PART_ID='Wood_part', XB=0.,1.,0.,1.,0.,1., N_PARTICLES=1, MASS_PER_VOLUME = 0.1 /

&DEVC XB=0.,1.,0.,1.,0.,1., QUANTITY='WALL TEMPERATURE', INIT_ID='init_1', ID='Particle T', SPATIAL_STATISTIC='MAX'/
&DEVC QUANTITY='WALL TEMPERATURE', INIT_ID='init_1', ID='Particle T', SPATIAL_STATISTIC='MAX'/
&DEVC XB=0.,1.,0.,1.,0.,1., QUANTITY='DENSITY', ID='fuel gas mass', SPATIAL_STATISTIC='VOLUME INTEGRAL', SPEC_ID='FUEL' /
&DEVC XB=0.,1.,0.,1.,0.,1., QUANTITY='DENSITY', ID='water vapor mass', SPATIAL_STATISTIC='VOLUME INTEGRAL', SPEC_ID='WATER VAPOR' /
&DEVC XB=0.,1.,0.,1.,0.,1., QUANTITY='MPUV', PART_ID='Wood_part', ID='solid mass', SPATIAL_STATISTIC='VOLUME INTEGRAL' /
Expand Down
2 changes: 1 addition & 1 deletion Verification/Pyrolysis/surf_mass_two_species_spher.fds
Original file line number Diff line number Diff line change
Expand Up @@ -55,7 +55,7 @@

&INIT ID = 'init_1', PART_ID='Wood_part', XB=0.,1.,0.,1.,0.,1., N_PARTICLES=1, MASS_PER_VOLUME = 0.1 /

&DEVC XB=0.,1.,0.,1.,0.,1., QUANTITY='WALL TEMPERATURE', INIT_ID='init_1', ID='Particle T', SPATIAL_STATISTIC='MAX'/
&DEVC QUANTITY='WALL TEMPERATURE', INIT_ID='init_1', ID='Particle T', SPATIAL_STATISTIC='MAX'/
&DEVC XB=0.,1.,0.,1.,0.,1., QUANTITY='DENSITY', ID='fuel gas mass', SPATIAL_STATISTIC='VOLUME INTEGRAL', SPEC_ID='FUEL' /
&DEVC XB=0.,1.,0.,1.,0.,1., QUANTITY='DENSITY', ID='water vapor mass', SPATIAL_STATISTIC='VOLUME INTEGRAL', SPEC_ID='WATER VAPOR' /
&DEVC XB=0.,1.,0.,1.,0.,1., QUANTITY='MPUV', PART_ID='Wood_part', ID='solid mass', SPATIAL_STATISTIC='VOLUME INTEGRAL' /
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -76,15 +76,15 @@

Ball output
&DEVC XB=0,1,0,1,0,1, QUANTITY='PARTICLE RADIATION LOSS', ID='Q_rad ball', SPATIAL_STATISTIC='VOLUME INTEGRAL', TEMPORAL_STATISTIC='TIME INTEGRAL'/
&DEVC XB=0,1,0,1,0,1, QUANTITY='WALL ENTHALPY', INIT_ID='ball', ID='Delta e_ball', RELATIVE=T /
&DEVC QUANTITY='WALL ENTHALPY', INIT_ID='ball', ID='Delta e_ball', RELATIVE=T /

Cylinder output
&DEVC XB=2,3,0,1,0,1, QUANTITY='PARTICLE RADIATION LOSS', ID='Q_rad cyl', SPATIAL_STATISTIC='VOLUME INTEGRAL' TEMPORAL_STATISTIC='TIME INTEGRAL'/
&DEVC XB=2,3,0,1,0,1, QUANTITY='WALL ENTHALPY', INIT_ID='cylinder', ID='Delta e_cyl', RELATIVE=T /
&DEVC QUANTITY='WALL ENTHALPY', INIT_ID='cylinder', ID='Delta e_cyl', RELATIVE=T /

Plate output
&DEVC XB=4,5,0,1,0,1, QUANTITY='PARTICLE RADIATION LOSS', ID='Q_rad plate', SPATIAL_STATISTIC='VOLUME INTEGRAL' TEMPORAL_STATISTIC='TIME INTEGRAL'/
&DEVC XB=4,5,0,1,0,1, QUANTITY='WALL ENTHALPY', INIT_ID='plate', ID='Delta e_plate', RELATIVE=T /
&DEVC QUANTITY='WALL ENTHALPY', INIT_ID='plate', ID='Delta e_plate', RELATIVE=T /

Plate object output

Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -67,15 +67,15 @@

Ball output
&DEVC XB=0,1,0,1,0,1, QUANTITY='PARTICLE RADIATION LOSS', ID='Q_rad ball', SPATIAL_STATISTIC='VOLUME INTEGRAL', TEMPORAL_STATISTIC='TIME INTEGRAL'/
&DEVC XB=0,1,0,1,0,1, QUANTITY='WALL ENTHALPY', INIT_ID='ball', ID='Delta e_ball', RELATIVE=T /
&DEVC QUANTITY='WALL ENTHALPY', INIT_ID='ball', ID='Delta e_ball', RELATIVE=T /

Cylinder output
&DEVC XB=2,3,0,1,0,1, QUANTITY='PARTICLE RADIATION LOSS', ID='Q_rad cyl', SPATIAL_STATISTIC='VOLUME INTEGRAL' TEMPORAL_STATISTIC='TIME INTEGRAL'/
&DEVC XB=2,3,0,1,0,1, QUANTITY='WALL ENTHALPY', INIT_ID='cylinder', ID='Delta e_cyl', RELATIVE=T /
&DEVC QUANTITY='WALL ENTHALPY', INIT_ID='cylinder', ID='Delta e_cyl', RELATIVE=T /

Plate output
&DEVC XB=4,5,0,1,0,1, QUANTITY='PARTICLE RADIATION LOSS', ID='Q_rad plate', SPATIAL_STATISTIC='VOLUME INTEGRAL' TEMPORAL_STATISTIC='TIME INTEGRAL'/
&DEVC XB=4,5,0,1,0,1, QUANTITY='WALL ENTHALPY', INIT_ID='plate', ID='Delta e_plate', RELATIVE=T /
&DEVC QUANTITY='WALL ENTHALPY', INIT_ID='plate', ID='Delta e_plate', RELATIVE=T /

Plate object output

Expand Down