Merge pull request #14428 from mcgratta/master

FDS Source: Issue #14221. Reduce need to specify BACK props

Author: mcgratta

Source/wall.f90

@@ -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)
@@ -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
@@ -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
@@ -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.
 
@@ -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
 
@@ -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
@@ -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
@@ -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

Verification/Heat_Transfer/ht1d_pile.fds

@@ -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' /

Verification/Heat_Transfer/ht3d_ibeam.fds

@@ -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/

Verification/Heat_Transfer/ht3d_pile.fds

@@ -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' /

Verification/Heat_Transfer/ht3d_sphere_24.fds

@@ -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 /

Verification/Heat_Transfer/ht3d_sphere_48.fds

@@ -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 /

Verification/Heat_Transfer/ht3d_sphere_96.fds

@@ -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 /

Verification/Pyrolysis/part_baking_soda_420K.fds

@@ -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))

Verification/Pyrolysis/part_baking_soda_450K.fds

@@ -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))

Verification/Pyrolysis/part_baking_soda_500K.fds

@@ -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))