Problemi nella sostituizione dei duct con un tratto 3D

[Francesco Roga]

Buongiorno, ho fatto girare questo caso, sostituendo due duct al tratto finale di un tunnel. Ho notato che a lungo andare si presenta un'instabilità numero e il comportamento del modello con i due duct non rispecchia quello del modello completamente 3D, infatti l'andamento dei fumi nei due modelli è completamente diverso. Alla luce di questo, vedendo il caso proposto sotto, ho dedotto che non è possibile sostituire un tratto monodimensionale ad uno tridimensionale mediante i duct. Avete riscontrato lo stesso problema ?

&HEAD CHID='prova_a_caso_2_1'/

&TIME T_END=500.0/

&DUMP RENDER_FILE='prova_a_caso_2_1.ge1', COLUMN_DUMP_LIMIT=.TRUE., DT_RESTART=300.0, DT_SL3D=0.25/

&MISC GVEC=0.0,-9.81,0.0/

&MESH ID='Mesh01', IJK=160,16,16, XB=-50.0,50.0,0.0,10.0,0.0,10.0/

&MESH ID='Mesh05', IJK=10,32,32, XB=46.875,50.0,0.0,10.0,0.0,10.0/

&REAC ID='POLYURETHANE',

      FYI='NFPA Babrauskas',

      FUEL='REAC_FUEL',

      C=6.3,

      H=7.1,

      O=2.1,

      N=1.0,

      AUTO_IGNITION_TEMPERATURE=0.0,

      SOOT_YIELD=1.0E-3,

      RADIATIVE_FRACTION=0.35/

&SURF ID='Surface01',

      COLOR='RED',

      HRRPUA=654.0,

      TMP_FRONT=300.0/

&VENT ID='Vent', SURF_ID='Surface01', XB=-6.0,2.0,0.0,0.0,5.0,8.0/ 

&VENT ID='Vent01', SURF_ID='OPEN', XB=-50.0,-50.0,0.0,10.0,0.0,10.0/ 

&VENT ID='Vent04', SURF_ID='HVAC', XB=50.0,50.0,5.0,10.0,0.0,10.0/ 

&VENT ID='Vent13', SURF_ID='HVAC', XB=50.0,50.0,0.0,5.0,0.0,10.0/ 

&HVAC ID='Node', TYPE_ID='NODE', DUCT_ID='Duct', VENT_ID='Vent04'/

&HVAC ID='Node01', TYPE_ID='NODE', DUCT_ID='Duct', AMBIENT=.TRUE., XYZ=100.0,7.5,5.0/

&HVAC ID='Duct', TYPE_ID='DUCT', AREA=50.0, PERIMETER=30.0, NODE_ID='Node','Node01', LENGTH=50.0/

&HVAC ID='Node02', TYPE_ID='NODE', DUCT_ID='Duct01', AMBIENT=.TRUE., XYZ=100.0,2.5,5.0/

&HVAC ID='Node03', TYPE_ID='NODE', DUCT_ID='Duct01', VENT_ID='Vent13'/

&HVAC ID='Duct01', TYPE_ID='DUCT', AREA=50.0, PERIMETER=30.0, NODE_ID='Node03','Node02', LENGTH=50.0/

&SLCF QUANTITY='TEMPERATURE', VECTOR=.TRUE., PBZ=5.0/

&SLCF QUANTITY='U-VELOCITY', VECTOR=.TRUE., PBZ=5.0/

&SLCF QUANTITY='PRESSURE', VECTOR=.TRUE., PBZ=5.0/

&SLCF QUANTITY='VISIBILITY', VECTOR=.TRUE., PBZ=5.0/

&SLCF QUANTITY='TEMPERATURE', VECTOR=.TRUE., PBX=50.0/

&SLCF QUANTITY='PRESSURE', VECTOR=.TRUE., PBX=50.0/

&SLCF QUANTITY='U-VELOCITY', VECTOR=.TRUE., PBX=50.0/

&TAIL /