tunnel

[h7899]

I have a question, can't a tunnel with vertical shafts be constructed using grid boundaries? It's like not being able to add extensions to both ends of the tunnel through grid boundaries. I hope you can answer this for me. Thank you!

[mcgratta]

I do not understand the question. Can you attach a simple drawing?

[h7899]

Hello, it's like this code. I added extensions at both ends of the tunnel, but it doesn't seem to meet the usage criteria of TUNNEL-PRECONDITIONER=T. What should I do.
&HEAD CHID='1', TITLE='2'/

&MESH XB=50.6,63.1,0.0,4.5,0.0,2.3, IJK=125,46,46, MULT_ID='mesh1' /
&MULT ID='mesh1', DX=12.5, I_UPPER=15 /1474560
&MESH ID='1' XB=250.6 263 0 4.5 0 2.3 IJK=125 46 46/
&MESH ID='2' XB=38.2 50.6 0 4.5 0 2.3 IJK=125 46 46/
&MESH ID='3' XB=250.6 263 0 4.5 2.3 4.6 IJK=125 46 46/
&MESH ID='4' XB=38.2 50.6 0 4.5 2.3 4.6 IJK=125 46 46/

&MISC TMPA=20 /

&TIME T_END=2000. /
&DUMP NFRAMES=2000, DT_HRR=1., DT_DEVC=1., /

&PRES TUNNEL_PRECONDITIONER=T, MAX_PRESSURE_ITERATIONS=100 /
&CLIP MINIMUM_DENSITY=0.01/

&VENT PBX= 38.2, SURF_ID='OPEN' /
&VENT PBX=263, SURF_ID='OPEN' /
&VENT PBZ=4.6, SURF_ID='OPEN' /
&VENT XB=50.6 250.6 0 0 0 2.3, SURF_ID='CONCRETE'/
&VENT XB=50.6 250.6 4.5 4.5 0 2.3, SURF_ID='CONCRETE' /
&VENT XB=38.2 50.6 4.5 4.5 0 4.6,SURF_ID='OPEN'/
&VENT XB=38.2 50.6 0 0 0 4.6,SURF_ID='OPEN'/
&VENT XB=250.6 263 4.5 4.5 0 4.6,SURF_ID='OPEN'/
&VENT XB=250.6 263 0 0 0 4.6,SURF_ID='OPEN'/

&MATL ID ='CONCRETE',
SPECIFIC_HEAT=1.22,
CONDUCTIVITY_RAMP = 'wood_k'
DENSITY=360/
&RAMP ID = 'wood_k', T = 25, F = 0.06/
&RAMP ID = 'wood_k', T = 200, F = 0.06/
&RAMP ID = 'wood_k', T = 400, F = 0.09/
&RAMP ID = 'wood_k', T = 600, F = 0.12/
&RAMP ID = 'wood_k', T = 800, F = 0.15/

&SURF ID ='CONCRETE',
BACKING='INSULATED',
MATL_ID(1,1)='CONCRETE',
MATL_MASS_FRACTION(1,1)=1.0,
THICKNESS(1)=0.05/

&SPEC ID ='PVC', FORMULA = 'C2H3Cl' /

&SPEC ID ='AIR', BACKGROUND=.TRUE.,
SPEC_ID(1)='OXYGEN',
SPEC_ID(2)='NITROGEN',
VOLUME_FRACTION(1)=1.545,
VOLUME_FRACTION(2)=5.809/
##生成物
&SPEC ID = 'HYDROGEN CHLORIDE', LUMPED_COMPONENT_ONLY = T /
&SPEC ID = 'SOOT', FORMULA='C1', LUMPED_COMPONENT_ONLY = F /
&SPEC ID ='CARBON DIOXIDE', LUMPED_COMPONENT_ONLY=.TRUE./
&SPEC ID ='CARBON MONOXIDE', LUMPED_COMPONENT_ONLY=.TRUE./
&SPEC ID ='WATER VAPOR', LUMPED_COMPONENT_ONLY=.TRUE./
&SPEC ID ='NITROGEN', LUMPED_COMPONENT_ONLY=.TRUE./
&SPEC ID ='OXYGEN', LUMPED_COMPONENT_ONLY=.TRUE./

PVC

&SPEC ID ='PRODUCT_3',
SPEC_ID(1)='HYDROGEN CHLORIDE', VOLUME_FRACTION(1)=1.0,
SPEC_ID(2)='WATER VAPOR', VOLUME_FRACTION(2)=1.0,
SPEC_ID(3)='CARBON MONOXIDE', VOLUME_FRACTION(3)=0.11,
SPEC_ID(4)='CARBON DIOXIDE', VOLUME_FRACTION(4)=0.99,
SPEC_ID(5)='NITROGEN', VOLUME_FRACTION(5)=5.809
SPEC_ID(6)='SOOT', VOLUME_FRACTION(6)=0.9 /

&REAC FUEL='PVC',
HEAT_OF_COMBUSTION=17700,
SOOT_YIELD=0.172
CO_YIELD=0.05
RADIATIVE_FRACTION=0.456
SPEC_ID_NU='PVC','AIR','PRODUCT_3',
NU=-1,-1, 1 /

&MATL ID = 'pvc'
DENSITY = 1500
EMISSIVITY = 0.9
CONDUCTIVITY = 0.026
SPECIFIC_HEAT = 6.266
HEAT_OF_COMBUSTION = 2025/

&MATL ID = 'T'
DENSITY = 8940
CONDUCTIVITY = 387
SPECIFIC_HEAT = 0.38/

&SURF ID = 'DL'
COLOR = 'PURPLE'
BURN_AWAY = .TRUE.
MATL_ID(1:2,1) = 'pvc','T'
THICKNESS(1:2) = 0.006,0.004
HRRPUA = 270
RAMP_Q = 'ramp_q'
IGNITION_TEMPERATURE = 419/
&RAMP ID = 'ramp_q', T = 0, F = 0/
&RAMP ID = 'ramp_q', T =30, F = 1/
&RAMP ID = 'ramp_q', T =270, F = 1/
&RAMP ID = 'ramp_q', T =510, F = 0.5/
&RAMP ID = 'ramp_q', T =750, F = 0/
&SURF ID='fire',HRRPUA=9881,COLOR='RED',
BURN_DURATION=158/

&OBST ID = 'ra1 ' , XB = 55 245 0.2 0.3 0.6 0.7 SURF_ID = 'DL'
MULT_ID='1'/
&OBST ID = 'ra2 ' , XB = 55 245 4.2 4.3 0.6 0.7 SURF_ID = 'DL'
MULT_ID='2'/
&MULT ID='1'
DY=0.2,J_UPPER=6
DZ=0.3,K_UPPER=4/
&MULT ID='2'
DY=-0.2,J_UPPER=6
DZ=0.3,K_UPPER=4/

&OBST ID='burn1', XB=149.8 150.2 0.7 1.1 0.3 0.4 SURF_ID6='INERT','INERT','INERT','INERT','INERT','fire'/
&SLCF PBY=2.3, QUANTITY='TEMPERATURE', CELL_CENTERED=.TRUE. /
&SLCF PBY=2.3, QUANTITY='VELOCITY', VECTOR=.TRUE. /
&SLCF PBY=2.3, QUANTITY='VOLUME FRACTION', SPEC_ID='OXYGEN', CELL_CENTERED=.TRUE. /
&SLCF PBY=2.3, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', CELL_CENTERED=.TRUE. /
&SLCF PBY=2.3, QUANTITY='PRESSURE', CELL_CENTERED=.TRUE. /

&TAIL /

[mcgratta]

I will take a look at the case.

[mcgratta]

You cannot use TUNNEL_PRECONDITIONER in this case. You can increase the pressure iterations, which you have done, or you can remove the extra meshes at the end.

[mcgratta]

The case failed for me after 100 s. I suggest that you remove the extensions at the end of the tunnel and use TUNNEL_PRECONDITIONER=T.