Sostituzione 3-D con 1-D

[Francesco Roga]

Buongiorno

Sto studiando un tunnel ideale di 40m. Lo scopo è quello di passare dal modello 3-D ad uno 1-D sostituendo un tratto del tunnel con un damper che riproduca gli effetti del tratto sostituito. Tuttavia a inizio prova (10s circa) mi riporta un'instabilità numerica. Ho notato che il problema è dovuto alle dimensioni del damper, le quali sono : area 25m^2 (sto considerando una sezione quadrata di lato 5m ), perimetro 20m, lunghezza fissata (del tratto che ha sostituito) di 10m. Nelle impostazioni ho messo una sezione non circolre.Per dimensioni minori, i test danno risultati. Sapreste spiegarmi perchè il 1-D non coincide con il 3-D, nonostante sia state rispettate le dimensioni ?

[luca.mi...@gmail.com]

Prova a mandarci il tuo file fds, se no risulta difficile aiutarti.

[Francesco Roga]

come posso inviare il file? Mette che posso inviare solo foto

[luca.mi...@gmail.com]

o fai un file winrar con dentro il file fds oppure copi direttamente qui il contenuto.

[Francesco Roga]

il problema è proprio questo, ho fatto entrambe le cose che lei ha detto, ma nel caso di file winrar mi dice che si è verificato un errore. Nel caso di copia diretta ,mi apre il file in basso, ma non nella conversazione

[luca.mi...@gmail.com]

No. intendo che potresti copiare proprio il testo del file fds qui. Non mi sembra un file complesso. non credo siano tante righe

[Francesco Roga]

2JF_1D_PROVA.fds

Generated by PyroSim - Version 2019.1.0515

1-mar-2021 13.50.49

&HEAD CHID='2JF_1D_PROVA'/

&TIME T_END=3600.0/

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

&MESH ID='Mesh01', IJK=33,10,10, XB=0.0,15.0,0.0,5.0,0.0,5.0/

&MESH ID='Mesh02', IJK=33,10,10, XB=25.0,40.0,0.0,5.0,0.0,5.0/

&REAC ID='POLYURETHANE',

      FYI='NFPA Babrauskas',

      FUEL='REAC_FUEL',

      C=6.3,

      H=7.1,

      O=2.1,

      N=1.0,

      SOOT_YIELD=0.1,

      RADIATIVE_FRACTION=0.35/

&DEVC ID='FED', QUANTITY='FED', XYZ=1.0,3.0,1.0/

&DEVC ID='FED01', QUANTITY='FED', XYZ=5.0,3.0,1.0/

&DEVC ID='FED02', QUANTITY='FED', XYZ=9.0,3.0,1.0/

&DEVC ID='FED03', QUANTITY='FED', XYZ=13.0,3.0,1.0/

&DEVC ID='FED06', QUANTITY='FED', XYZ=25.0,3.0,1.0/

&DEVC ID='FED07', QUANTITY='FED', XYZ=29.0,3.0,1.0/

&DEVC ID='FED08', QUANTITY='FED', XYZ=33.0,3.0,1.0/

&DEVC ID='FED09', QUANTITY='FED', XYZ=37.0,3.0,1.0/

&DEVC ID='FED10', QUANTITY='FED', XYZ=1.0,3.0,2.5/

&DEVC ID='FED11', QUANTITY='FED', XYZ=5.0,3.0,2.5/

&DEVC ID='FED12', QUANTITY='FED', XYZ=9.0,3.0,2.5/

&DEVC ID='FED13', QUANTITY='FED', XYZ=13.0,3.0,2.5/

&DEVC ID='FED16', QUANTITY='FED', XYZ=25.0,3.0,2.5/

&DEVC ID='FED17', QUANTITY='FED', XYZ=29.0,3.0,2.5/

&DEVC ID='FED18', QUANTITY='FED', XYZ=33.0,3.0,2.5/

&DEVC ID='FED19', QUANTITY='FED', XYZ=37.0,3.0,2.5/

&DEVC ID='FED20', QUANTITY='FED', XYZ=1.0,3.0,4.0/

&DEVC ID='FED21', QUANTITY='FED', XYZ=5.0,3.0,4.0/

&DEVC ID='FED22', QUANTITY='FED', XYZ=9.0,3.0,4.0/

&DEVC ID='FED23', QUANTITY='FED', XYZ=13.0,3.0,4.0/

&DEVC ID='FED26', QUANTITY='FED', XYZ=25.0,3.0,4.0/

&DEVC ID='FED27', QUANTITY='FED', XYZ=29.0,3.0,4.0/

&DEVC ID='FED28', QUANTITY='FED', XYZ=33.0,3.0,4.0/

&DEVC ID='FED29', QUANTITY='FED', XYZ=37.0,3.0,4.0/

&DEVC ID='temp', QUANTITY='TEMPERATURE', XYZ=1.0,3.0,1.0/

&DEVC ID='temp01', QUANTITY='TEMPERATURE', XYZ=5.0,3.0,1.0/

&DEVC ID='temp02', QUANTITY='TEMPERATURE', XYZ=9.0,3.0,1.0/

&DEVC ID='temp03', QUANTITY='TEMPERATURE', XYZ=13.0,3.0,1.0/

&DEVC ID='temp06', QUANTITY='TEMPERATURE', XYZ=25.0,3.0,1.0/

&DEVC ID='temp07', QUANTITY='TEMPERATURE', XYZ=29.0,3.0,1.0/

&DEVC ID='temp08', QUANTITY='TEMPERATURE', XYZ=33.0,3.0,1.0/

&DEVC ID='temp09', QUANTITY='TEMPERATURE', XYZ=37.0,3.0,1.0/

&DEVC ID='temp10', QUANTITY='TEMPERATURE', XYZ=1.0,3.0,2.5/

&DEVC ID='temp11', QUANTITY='TEMPERATURE', XYZ=5.0,3.0,2.5/

&DEVC ID='temp12', QUANTITY='TEMPERATURE', XYZ=9.0,3.0,2.5/

&DEVC ID='temp13', QUANTITY='TEMPERATURE', XYZ=13.0,3.0,2.5/

&DEVC ID='temp16', QUANTITY='TEMPERATURE', XYZ=25.0,3.0,2.5/

&DEVC ID='temp17', QUANTITY='TEMPERATURE', XYZ=29.0,3.0,2.5/

&DEVC ID='temp18', QUANTITY='TEMPERATURE', XYZ=33.0,3.0,2.5/

&DEVC ID='temp19', QUANTITY='TEMPERATURE', XYZ=37.0,3.0,2.5/

&DEVC ID='temp20', QUANTITY='TEMPERATURE', XYZ=1.0,3.0,4.0/

&DEVC ID='temp21', QUANTITY='TEMPERATURE', XYZ=5.0,3.0,4.0/

&DEVC ID='temp22', QUANTITY='TEMPERATURE', XYZ=9.0,3.0,4.0/

&DEVC ID='temp23', QUANTITY='TEMPERATURE', XYZ=13.0,3.0,4.0/

&DEVC ID='temp26', QUANTITY='TEMPERATURE', XYZ=25.0,3.0,4.0/

&DEVC ID='temp27', QUANTITY='TEMPERATURE', XYZ=29.0,3.0,4.0/

&DEVC ID='temp28', QUANTITY='TEMPERATURE', XYZ=33.0,3.0,4.0/

&DEVC ID='temp29', QUANTITY='TEMPERATURE', XYZ=37.0,3.0,4.0/

&DEVC ID='vel', QUANTITY='VELOCITY', XYZ=1.0,3.0,1.0/

&DEVC ID='vel01', QUANTITY='VELOCITY', XYZ=5.0,3.0,1.0/

&DEVC ID='vel02', QUANTITY='VELOCITY', XYZ=9.0,3.0,1.0/

&DEVC ID='vel03', QUANTITY='VELOCITY', XYZ=13.0,3.0,1.0/

&DEVC ID='vel06', QUANTITY='VELOCITY', XYZ=25.0,3.0,1.0/

&DEVC ID='vel07', QUANTITY='VELOCITY', XYZ=29.0,3.0,1.0/

&DEVC ID='vel08', QUANTITY='VELOCITY', XYZ=33.0,3.0,1.0/

&DEVC ID='vel09', QUANTITY='VELOCITY', XYZ=37.0,3.0,1.0/

&DEVC ID='vel10', QUANTITY='VELOCITY', XYZ=1.0,3.0,2.5/

&DEVC ID='vel11', QUANTITY='VELOCITY', XYZ=5.0,3.0,2.5/

&DEVC ID='vel12', QUANTITY='VELOCITY', XYZ=9.0,3.0,2.5/

&DEVC ID='vel13', QUANTITY='VELOCITY', XYZ=13.0,3.0,2.5/

&DEVC ID='vel16', QUANTITY='VELOCITY', XYZ=25.0,3.0,2.5/

&DEVC ID='vel17', QUANTITY='VELOCITY', XYZ=29.0,3.0,2.5/

&DEVC ID='vel18', QUANTITY='VELOCITY', XYZ=33.0,3.0,2.5/

&DEVC ID='vel19', QUANTITY='VELOCITY', XYZ=37.0,3.0,2.5/

&DEVC ID='vel20', QUANTITY='VELOCITY', XYZ=1.0,3.0,4.0/

&DEVC ID='vel21', QUANTITY='VELOCITY', XYZ=5.0,3.0,4.0/

&DEVC ID='vel22', QUANTITY='VELOCITY', XYZ=9.0,3.0,4.0/

&DEVC ID='vel23', QUANTITY='VELOCITY', XYZ=13.0,3.0,4.0/

&DEVC ID='vel26', QUANTITY='VELOCITY', XYZ=25.0,3.0,4.0/

&DEVC ID='vel27', QUANTITY='VELOCITY', XYZ=29.0,3.0,4.0/

&DEVC ID='vel28', QUANTITY='VELOCITY', XYZ=33.0,3.0,4.0/

&DEVC ID='vel29', QUANTITY='VELOCITY', XYZ=37.0,3.0,4.0/

&DEVC ID='vel output', QUANTITY='VELOCITY', XYZ=7.2,4.2,2.45/

&DEVC ID='vel output01', QUANTITY='VELOCITY', XYZ=27.2,4.2,2.45/

&SURF ID='burn',

      COLOR='RED',

      HRRPUA=100.0,

      RAMP_Q='burn_RAMP_Q',

      TMP_FRONT=300.0/

&RAMP ID='burn_RAMP_Q', T=120.0, F=0.6/

&RAMP ID='burn_RAMP_Q', T=200.0, F=1.0/

&OBST ID='Obstruction', XB=6.0,6.5,4.0,4.3,2.5,2.55, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=6.0,6.5,4.0,4.3,2.3,2.35, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=6.0,6.5,4.0,4.05,2.3,2.55, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=6.0,6.5,4.3,4.35,2.3,2.55, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=6.5,6.55,4.0,4.35,2.3,2.55, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=6.55,7.05,4.0,4.3,2.5,2.55, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=6.55,7.05,4.0,4.3,2.3,2.35, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=6.55,7.05,4.0,4.05,2.3,2.55, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=6.55,7.05,4.3,4.35,2.3,2.55, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=26.0,26.5,4.0,4.3,2.5,2.55, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=26.0,26.5,4.0,4.3,2.3,2.35, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=26.0,26.5,4.0,4.05,2.3,2.55, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=26.0,26.5,4.3,4.35,2.3,2.55, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=26.5,26.55,4.0,4.35,2.3,2.55, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=26.55,27.05,4.0,4.3,2.5,2.55, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=26.55,27.05,4.0,4.3,2.3,2.35, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=26.55,27.05,4.0,4.05,2.3,2.55, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=26.55,27.05,4.3,4.35,2.3,2.55, SURF_ID='INERT'/ 

&VENT ID='Vent', SURF_ID='OPEN', XB=0.0,0.0,4.973799E-14,5.0,-1.953993E-14,5.0/ 

&VENT ID='Vent11', SURF_ID='OPEN', XB=40.0,40.0,3.019807E-14,5.0,-3.019807E-14,5.0/ 

&VENT ID='Vent13', SURF_ID='burn', XB=6.666667,10.30303,0.0,0.0,1.5,3.5/ 

&VENT ID='Vent17', SURF_ID='HVAC', XB=6.5,6.5,4.0,4.35,2.3,2.55/ 

&VENT ID='Vent19', SURF_ID='HVAC', XB=6.55,6.55,4.0,4.35,2.3,2.55/ 

&VENT ID='Vent20', SURF_ID='HVAC', XB=26.5,26.5,4.0,4.35,2.3,2.55/ 

&VENT ID='Vent21', SURF_ID='HVAC', XB=26.55,26.55,4.0,4.35,2.3,2.55/ 

&VENT ID='Vent22', SURF_ID='HVAC', XB=15.0,15.0,4.973799E-14,5.0,-1.953993E-14,5.0/ 

&VENT ID='Vent23', SURF_ID='HVAC', XB=25.0,25.0,4.973799E-14,5.0,-1.953993E-14,5.0/ 

&HVAC ID='Node', TYPE_ID='NODE', DUCT_ID='Duct03', AMBIENT=.TRUE., XYZ=6.5,4.2,2.4/

&HVAC ID='Node04', TYPE_ID='NODE', DUCT_ID='Duct03', AMBIENT=.TRUE., XYZ=6.55,4.2,2.4/

&HVAC ID='Duct03', TYPE_ID='DUCT', DIAMETER=0.3048, VOLUME_FLOW=1.5E-3, NODE_ID='Node','Node04', ROUGHNESS=1.0E-3, LENGTH=0.05/

&HVAC ID='Node05', TYPE_ID='NODE', DUCT_ID='Duct04', AMBIENT=.TRUE., XYZ=26.5,4.2,2.4/

&HVAC ID='Node06', TYPE_ID='NODE', DUCT_ID='Duct04', AMBIENT=.TRUE., XYZ=26.55,4.2,2.4/

&HVAC ID='Duct04', TYPE_ID='DUCT', DIAMETER=0.3048, VOLUME_FLOW=1.5E-3, NODE_ID='Node05','Node06', ROUGHNESS=1.0E-3, LENGTH=0.05/

&HVAC ID='Node07', TYPE_ID='NODE', DUCT_ID='Duct', VENT_ID='Vent23'/

&HVAC ID='Node08', TYPE_ID='NODE', DUCT_ID='Duct', VENT_ID='Vent22'/

&HVAC ID='Duct', TYPE_ID='DUCT', DIAMETER=5.0, DAMPER=.TRUE., NODE_ID='Node08','Node07', ROUGHNESS=0.0, LENGTH=10.0/

&SLCF QUANTITY='MASS FRACTION', SPEC_ID='CARBON MONOXIDE', VECTOR=.TRUE., PBX=1.0/

&SLCF QUANTITY='MASS FRACTION', SPEC_ID='CARBON MONOXIDE', VECTOR=.TRUE., PBY=1.0/

&SLCF QUANTITY='MASS FRACTION', SPEC_ID='CARBON MONOXIDE', VECTOR=.TRUE., PBZ=1.0/

&SLCF QUANTITY='MASS FRACTION', SPEC_ID='CARBON MONOXIDE', VECTOR=.TRUE., PBX=2.0/

&SLCF QUANTITY='VELOCITY', VECTOR=.TRUE., PBX=7.2/

&SLCF QUANTITY='VELOCITY', VECTOR=.TRUE., PBY=4.2/

&SLCF QUANTITY='VELOCITY', VECTOR=.TRUE., PBZ=2.45/

&SLCF QUANTITY='VELOCITY', VECTOR=.TRUE., PBX=27.2/

&TAIL /

purtroppo non mi fa copiare nessu n file al di fuori delle immagini, il che non ha neanche tanto senso, è questo quello che le serve ?

[luca.mi...@gmail.com]

Si è questo. Ma a te questo file da un errore?

[Francesco Roga]

no allora questo va ma perchè ho messo valori dell'area e del perimetro minori (ritornando al primo messaggio che ho scritto). Mettendo quelli che dovrebbero essere i valori giusti nel dumper : area 25m^2 e perimetro 20m, allora mi da errore

[luca.mi...@gmail.com]

postami il caso che ti da errore allora, perchè io ho provato a fare questa modifica e non mi da errori. Per i prossimi invii di listati, ti chiedo di postare file ripuliti, con solo il minimo indispensabile per cogliere il problema che manifesti.

[Francesco Roga]

&HVAC ID='Node', TYPE_ID='NODE', DUCT_ID='Duct03', AMBIENT=.TRUE., XYZ=6.5,4.2,2.4/

&HVAC ID='Node04', TYPE_ID='NODE', DUCT_ID='Duct03', AMBIENT=.TRUE., XYZ=6.55,4.2,2.4/

&HVAC ID='Duct03', TYPE_ID='DUCT', DIAMETER=0.3048, VOLUME_FLOW=5.0E-3, NODE_ID='Node','Node04', ROUGHNESS=1.0E-3, LENGTH=0.05/

&HVAC ID='Node05', TYPE_ID='NODE', DUCT_ID='Duct04', AMBIENT=.TRUE., XYZ=26.5,4.2,2.4/

&HVAC ID='Node06', TYPE_ID='NODE', DUCT_ID='Duct04', AMBIENT=.TRUE., XYZ=26.55,4.2,2.4/

&HVAC ID='Duct04', TYPE_ID='DUCT', DIAMETER=0.3048, VOLUME_FLOW=5.0E-3, NODE_ID='Node05','Node06', ROUGHNESS=1.0E-3, LENGTH=0.05/

&HVAC ID='Node07', TYPE_ID='NODE', DUCT_ID='Duct', VENT_ID='Vent23'/

&HVAC ID='Node08', TYPE_ID='NODE', DUCT_ID='Duct', VENT_ID='Vent22'/

&HVAC ID='Duct', TYPE_ID='DUCT', AREA=25.0, PERIMETER=20.0, DAMPER=.TRUE., NODE_ID='Node08','Node07', LENGTH=10.0/

l'errore dovrebbe stare qui dentro perchè qui tratta il duct in cui imposto le dimensioni. Giusto?

[luca.mi...@gmail.com]

A me non da problemi. Direi che è un problema di versione di FDS. Devi usare l'ultima versione, mentre ho visto stai usando pyrosim 2019. Avvia il file fds direttamente da CMDfds sul desktop e vedi se ti da errore.   Scusa se te lo chiedo, è un progetto reale o una tesi?

[Francesco Roga]

è un progetto di tesi

[luca.mi...@gmail.com]

Ok. Fammi sapere se facendo come ti ho detto il problema scompare. Se ti da ancora errore fammi uno screenshot al finestra di CMDfds.

[Francesco Roga]

buon pomeriggio, sto vedendo un pò le cose che lei ha detto, unica problematica è che dovrei continuare a lavorare col pyrosim 2019 per effettuare un confronto di valori al termine delle valutazioni. Avrei un'altra domanda. E' possibile sostituire un modello 3-D con uno 1-D mediante un duct? perchè ho notato che il passaggio del fumo dal primo modulo (prima mesh) al secondo ( come mostrato nella figura di sopra) è praticamente istantaneo e quindi non veritiero. Sa a cosa potrebbe essere dovuto?

[luca.mi...@gmail.com]

pyrosim è solo una interfaccia grafica, puoi lavorare con pyrosim 2019 (o anche blender che è anche opensource e te lo consiglio dato che stai lavorando in ambito universitario) e allo stesso tempo inizializzare le simulazioni sulla versione di fds più aggiornata.  Una non preclude l'altra.   Ricorda che il software CFD è fds e non Pyrosim.  Ho già visto in alcune pubblicazioni sostituire il modello 3D con 1D, ma in quel caso il condotto 1D non comunicava a sua volta con  secondo dominio  3D,  ma smaltiva direttamente i fumi (questo fa passare in secondo piano il tempo di trasporto nel condotto). Inoltre queste pubblicazioni non usavano proprio FDS, ma credo che non sia da escludere che anche fds lo possa fare. Credo che l'utilizzo di un condotto 1D alla fine del tunnel potrebbe essere una buona idea per avere un risparmio in termini computazionali, senza trascurare le condizioni al contorno.

Ti consiglio di leggere sulla User Guide il capitolo 12.12. Dato che comunque stai forzando il modello HVAC a svolgere una applicazione non tipica devi in qualche modo trovare una soluzione su misura per il tuo caso. Dato che sei agli inizi, ti consiglio di partire leggendo il capitolo dell'HVAC sulla user Guide, poi proverei a guardare sulla verification guide i casi già fatti dell'HVAC, che puoi far rigirare anche tu, li trovi qui https://github.com/firemodels/fds/tree/master/Verification/HVAC.  Inoltre questo modello, per come è stato fatto, dal mio punto di vista non ti aiuta a capire se le tempistiche di trasporto dei fumi vengono conteggiate. Infatti, la crescita del focolare è praticamente istantanea e 10 m sono pochissimi per valutare le tempistiche.

[Francesco Roga]

i file, per aprirli, li ho copiati in un documento, poi ho messo "apri con " pyrosim, ma mi dice che i file potrebbero essere corrotti. C'è un altro modo per aprirli?

[Francesco Roga]

no ok, ho risolto importando il file in pyrosim

[Francesco Roga]

comunque da pyrosim ho aperto il manuale e gli HVAC li tratta nel capitolo 17, nel 12 sono riportate le species. 

[luca.mi...@gmail.com]

I manuali si trovano qui https://pages.nist.gov/fds-smv/manuals.html, questo è il gruppo di discussione di FDS, non di Pyrosim.  Ripeto, Pyrosim è l'interfaccia grafica di FDS, ti conviene capire come funziona FDS prima di utilizzare Pyrosim.

[Francesco Roga]

Buonasera, avrei un'atra domanda, quando lancio il programma mi dice che i VENT 2 e VENT 4 non sono collegati ad una superficie. Ma non capisco quali siano questi vent(poichè nella ricerca non mi escono)  , inoltre mi risulta che tutti i vent siano collegati ad una superficie. 

tunnell 130m 3-D con JF.fds

Generated by PyroSim - Version 2019.1.0515

13-mar-2021 19.35.43

&HEAD CHID='tunnell_130m_3-D_con_JF'/

&TIME T_END=3600.0/

&DUMP RENDER_FILE='tunnell_130m_3-D_con_JF.ge1', COLUMN_DUMP_LIMIT=.TRUE., DT_RESTART=300.0, DT_SL3D=0.25/

&MESH ID='Mesh01', IJK=100,20,20, XB=0.0,50.0,0.0,10.0,0.0,10.0/

&MESH ID='Mesh02', IJK=100,20,20, XB=80.0,130.0,0.0,10.0,0.0,10.0/

&MESH ID='Mesh03', IJK=100,20,20, XB=50.0,80.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,

      SOOT_YIELD=0.1,

      RADIATIVE_FRACTION=0.35/

&DEVC ID='fed', QUANTITY='FED', XYZ=5.0,5.0,2.5/

&DEVC ID='fed01', QUANTITY='FED', XYZ=10.0,5.0,2.5/

&DEVC ID='fed02', QUANTITY='FED', XYZ=15.0,5.0,2.5/

&DEVC ID='fed03', QUANTITY='FED', XYZ=20.0,5.0,2.5/

&DEVC ID='fed04', QUANTITY='FED', XYZ=25.0,5.0,2.5/

&DEVC ID='fed05', QUANTITY='FED', XYZ=30.0,5.0,2.5/

&DEVC ID='fed06', QUANTITY='FED', XYZ=35.0,5.0,2.5/

&DEVC ID='fed07', QUANTITY='FED', XYZ=40.0,5.0,2.5/

&DEVC ID='fed08', QUANTITY='FED', XYZ=45.0,5.0,2.5/

&DEVC ID='fed09', QUANTITY='FED', XYZ=50.0,5.0,2.5/

&DEVC ID='fed10', QUANTITY='FED', XYZ=5.0,5.0,5.0/

&DEVC ID='fed11', QUANTITY='FED', XYZ=5.0,5.0,7.5/

&DEVC ID='fed12', QUANTITY='FED', XYZ=10.0,5.0,5.0/

&DEVC ID='fed13', QUANTITY='FED', XYZ=15.0,5.0,5.0/

&DEVC ID='fed14', QUANTITY='FED', XYZ=20.0,5.0,5.0/

&DEVC ID='fed15', QUANTITY='FED', XYZ=25.0,5.0,5.0/

&DEVC ID='fed16', QUANTITY='FED', XYZ=30.0,5.0,5.0/

&DEVC ID='fed17', QUANTITY='FED', XYZ=35.0,5.0,5.0/

&DEVC ID='fed18', QUANTITY='FED', XYZ=40.0,5.0,5.0/

&DEVC ID='fed19', QUANTITY='FED', XYZ=45.0,5.0,5.0/

&DEVC ID='fed20', QUANTITY='FED', XYZ=50.0,5.0,5.0/

&DEVC ID='fed21', QUANTITY='FED', XYZ=10.0,5.0,7.5/

&DEVC ID='fed22', QUANTITY='FED', XYZ=15.0,5.0,7.5/

&DEVC ID='fed23', QUANTITY='FED', XYZ=20.0,5.0,7.5/

&DEVC ID='fed24', QUANTITY='FED', XYZ=25.0,5.0,7.5/

&DEVC ID='fed25', QUANTITY='FED', XYZ=30.0,5.0,7.5/

&DEVC ID='fed26', QUANTITY='FED', XYZ=35.0,5.0,7.5/

&DEVC ID='fed27', QUANTITY='FED', XYZ=40.0,5.0,7.5/

&DEVC ID='fed28', QUANTITY='FED', XYZ=45.0,5.0,7.5/

&DEVC ID='fed29', QUANTITY='FED', XYZ=50.0,5.0,7.5/

&DEVC ID='fed30', QUANTITY='FED', XYZ=85.0,5.0,2.5/

&DEVC ID='fed31', QUANTITY='FED', XYZ=85.0,5.0,5.0/

&DEVC ID='fed32', QUANTITY='FED', XYZ=85.0,5.0,7.5/

&DEVC ID='fed33', QUANTITY='FED', XYZ=90.0,5.0,2.5/

&DEVC ID='fed34', QUANTITY='FED', XYZ=90.0,5.0,5.0/

&DEVC ID='fed35', QUANTITY='FED', XYZ=90.0,5.0,7.5/

&DEVC ID='fed36', QUANTITY='FED', XYZ=95.0,5.0,2.5/

&DEVC ID='fed37', QUANTITY='FED', XYZ=95.0,5.0,5.0/

&DEVC ID='fed38', QUANTITY='FED', XYZ=95.0,5.0,7.5/

&DEVC ID='fed39', QUANTITY='FED', XYZ=100.0,5.0,2.5/

&DEVC ID='fed40', QUANTITY='FED', XYZ=100.0,5.0,5.0/

&DEVC ID='fed41', QUANTITY='FED', XYZ=100.0,5.0,7.5/

&DEVC ID='fed42', QUANTITY='FED', XYZ=105.0,5.0,2.5/

&DEVC ID='fed43', QUANTITY='FED', XYZ=105.0,5.0,5.0/

&DEVC ID='fed44', QUANTITY='FED', XYZ=105.0,5.0,7.5/

&DEVC ID='fed45', QUANTITY='FED', XYZ=110.0,5.0,2.5/

&DEVC ID='fed46', QUANTITY='FED', XYZ=110.0,5.0,5.0/

&DEVC ID='fed47', QUANTITY='FED', XYZ=110.0,5.0,7.5/

&DEVC ID='fed48', QUANTITY='FED', XYZ=115.0,5.0,2.5/

&DEVC ID='fed49', QUANTITY='FED', XYZ=115.0,5.0,5.0/

&DEVC ID='fed50', QUANTITY='FED', XYZ=115.0,5.0,7.5/

&DEVC ID='fed51', QUANTITY='FED', XYZ=120.0,5.0,2.5/

&DEVC ID='fed52', QUANTITY='FED', XYZ=120.0,5.0,5.0/

&DEVC ID='fed53', QUANTITY='FED', XYZ=120.0,5.0,7.5/

&DEVC ID='fed54', QUANTITY='FED', XYZ=125.0,5.0,2.5/

&DEVC ID='fed55', QUANTITY='FED', XYZ=125.0,5.0,5.0/

&DEVC ID='fed56', QUANTITY='FED', XYZ=125.0,5.0,7.5/

&DEVC ID='fed57', QUANTITY='FED', XYZ=130.0,5.0,2.5/

&DEVC ID='fed58', QUANTITY='FED', XYZ=130.0,5.0,5.0/

&DEVC ID='fed59', QUANTITY='FED', XYZ=130.0,5.0,7.5/

&DEVC ID='fed60', QUANTITY='FED', XYZ=55.0,5.0,7.5/

&DEVC ID='fed61', QUANTITY='FED', XYZ=55.0,5.0,5.0/

&DEVC ID='fed62', QUANTITY='FED', XYZ=55.0,5.0,2.5/

&DEVC ID='fed63', QUANTITY='FED', XYZ=60.0,5.0,7.5/

&DEVC ID='fed64', QUANTITY='FED', XYZ=60.0,5.0,5.0/

&DEVC ID='fed65', QUANTITY='FED', XYZ=60.0,5.0,2.5/

&DEVC ID='fed66', QUANTITY='FED', XYZ=65.0,5.0,7.5/

&DEVC ID='fed67', QUANTITY='FED', XYZ=65.0,5.0,5.0/

&DEVC ID='fed68', QUANTITY='FED', XYZ=65.0,5.0,2.5/

&DEVC ID='fed69', QUANTITY='FED', XYZ=70.0,5.0,7.5/

&DEVC ID='fed70', QUANTITY='FED', XYZ=70.0,5.0,5.0/

&DEVC ID='fed71', QUANTITY='FED', XYZ=70.0,5.0,2.5/

&DEVC ID='fed72', QUANTITY='FED', XYZ=75.0,5.0,7.5/

&DEVC ID='fed73', QUANTITY='FED', XYZ=75.0,5.0,5.0/

&DEVC ID='fed74', QUANTITY='FED', XYZ=75.0,5.0,2.5/

&DEVC ID='fed75', QUANTITY='FED', XYZ=80.0,5.0,7.5/

&DEVC ID='fed76', QUANTITY='FED', XYZ=80.0,5.0,5.0/

&DEVC ID='fed77', QUANTITY='FED', XYZ=80.0,5.0,2.5/

&SURF ID='burn',

      COLOR='RED',

      HRRPUA=654.54,

      TMP_FRONT=300.0/

&OBST ID='Obstruction', XB=7.0,7.1,8.0,8.65,3.55,4.05, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=6.0,7.0,8.0,8.6,4.05,4.1, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=6.0,7.0,8.0,8.05,3.55,4.05, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=6.0,7.0,8.0,8.6,3.55,3.6, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=6.0,7.0,8.6,8.65,3.55,4.05, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=7.1,8.1,8.0,8.6,4.05,4.1, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=7.1,8.1,8.0,8.05,3.55,4.05, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=7.1,8.1,8.0,8.6,3.55,3.6, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=7.1,8.1,8.6,8.65,3.55,4.05, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=87.0,87.1,8.0,8.65,3.55,4.05, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=86.0,87.0,8.0,8.6,4.05,4.1, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=86.0,87.0,8.0,8.05,3.55,4.05, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=86.0,87.0,8.0,8.6,3.55,3.6, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=86.0,87.0,8.6,8.65,3.55,4.05, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=87.1,88.1,8.0,8.6,4.05,4.1, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=87.1,88.1,8.0,8.05,3.55,4.05, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=87.1,88.1,8.0,8.6,3.55,3.6, SURF_ID='INERT'/ 

&OBST ID='Obstruction', XB=87.1,88.1,8.6,8.65,3.55,4.05, SURF_ID='INERT'/ 

&VENT ID='Vent', SURF_ID='OPEN', XB=0.0,0.0,0.0,10.0,0.0,10.0/ 

&VENT ID='Vent17', SURF_ID='HVAC', XB=7.0,7.0,8.0,8.65,3.55,4.05, TEXTURE_ORIGIN=14.0,16.0,7.1/ 

&VENT ID='Vent18', SURF_ID='HVAC', XB=7.1,7.1,8.0,8.65,3.55,4.05, TEXTURE_ORIGIN=14.2,16.65,7.6/ 

&VENT ID='Vent19', SURF_ID='HVAC', XB=87.0,87.0,8.0,8.65,3.55,4.05, TEXTURE_ORIGIN=94.0,8.0,3.55/ 

&VENT ID='Vent20', SURF_ID='HVAC', XB=87.1,87.1,8.0,8.65,3.55,4.05/ 

&VENT ID='Vent13', SURF_ID='burn', XB=14.5,19.5,0.0,0.0,4.5,8.5/ 

&VENT ID='Vent14', SURF_ID='OPEN', XB=130.0,130.0,0.0,10.0,0.0,10.0/ 

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

&HVAC ID='Node02', TYPE_ID='NODE', DUCT_ID='Duct', VENT_ID='Vent18'/

&HVAC ID='Duct', TYPE_ID='DUCT', DIAMETER=0.3048, VOLUME_FLOW=5.0E-3, NODE_ID='Node','Node02', ROUGHNESS=0.0, LENGTH=0.1/

&HVAC ID='Node05', TYPE_ID='NODE', VENT_ID='Vent19'/

&HVAC ID='Node06', TYPE_ID='NODE', VENT_ID='Vent20'/

&TAIL /

[luca.mi...@gmail.com]

Quando ti trovi in queste situazioni, in cui non sai quale sia l'oggetto che causa problemi, l'unica maniera è eliminare gradualmente oggetti finchè non vedi che  la simulazione si avvia senza errori. In generale ti conviene testare ogni oggetto che inserisci e verificare che non ti dia errori prima di inserirlo in un caso di studio più grande. 

Quando utilizzi fds inoltre devi tenere sempre conto che se inserisci un oggetto di dimensioni inferiori alla tua cella, questi andranno ad occupare tutta la cella. Io consiglio di far coincidere sempre input e output, in altre parole disegna sempre le cose con dimensioni multiple della tua cella. Ad esempio vedo che le superfici parietali del tuo ventilatore sono più sottili delle dimensioni della cella.

[Francesco Roga]

ciao, ho fatto girare un altro caso, in cui ho collegato due duct dirttamente all'esterno, aumentando le dimensioni del tunnel, ma ho visto che comunque il comportamento dei due modelli , ossia quello completamente 3D e quello con tratto finale composto dai duct è diverso. Nella pratica il comportamento dei fumi è diverso, così come le velocità. Le allego il secondo modello, con i duct 

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