Skip to content

FDS Verification: Add duct_flow_uglmat_hypre.fds + dataplot figures. #13855

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
marcosvanella merged 2 commits into from
Dec 6, 2024
Merged

FDS Verification: Add duct_flow_uglmat_hypre.fds + dataplot figures. #13855

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
57b6397
FDS Verification: Add duct_flow_uglmat_hypre.fds + dataplot figures.
marcosvanella Dec 4, 2024
eb91793
Merge remote-tracking branch 'origin/master'
marcosvanella Dec 6, 2024
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
5 changes: 5 additions & 0 deletions Utilities/Matlab/FDS_verification_dataplot_inputs.csv
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -109,15 +109,18 @@ d,cyl_test,Flowfields/cyl_test_3_git.txt,Flowfields/cyl_test_3_devc.csv,2,3,Time
d,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,iter,Preconditioner,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_tight_no_precon_devc.csv,2,3,Time,iter,No Preconditioner,k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (dancing\_eddies),Time (s),Pressure Iterations,0,2,1,0,100,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,iter,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_uglmat_devc.csv,2,3,Time,iter,UGLMAT,r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (dancing\_eddies),Time (s),Pressure Iterations,0,2,1,0,100,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,iter,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_ulmat_devc.csv,2,3,Time,iter,ULMAT,b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (dancing\_eddies),Time (s),Pressure Iterations,0,2,1,0,100,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,iter,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_uglmat_hypre_devc.csv,2,3,Time,iter,UGLMAT HYPRE,m-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (dancing\_eddies),Time (s),Pressure Iterations,0,2,1,0,100,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,iter,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_ulmat_hypre_devc.csv,2,3,Time,iter,ULMAT HYPRE,g-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (dancing\_eddies),Time (s),Pressure Iterations,0,2,1,0,100,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
d,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,cpu,Preconditioner,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_tight_no_precon_devc.csv,2,3,Time,cpu,No Preconditioner,k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,CPU Time (dancing\_eddies),Simulated Time (s),CPU Time (s),0,2,1,0,400,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_tight_precon,N/A,end,0,Pressure Solver,kd,k,TeX
f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,cpu,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_uglmat_devc.csv,2,3,Time,cpu,UGLMAT,r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,CPU Time (dancing\_eddies),Simulated Time (s),CPU Time (s),0,2,1,0,400,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_tight_precon,N/A,end,0,Pressure Solver,kd,k,TeX
f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,cpu,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_ulmat_devc.csv,2,3,Time,cpu,ULMAT,b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,CPU Time (dancing\_eddies),Simulated Time (s),CPU Time (s),0,2,1,0,400,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_tight_precon,N/A,end,0,Pressure Solver,kd,k,TeX
f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,cpu,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_uglmat_hypre_devc.csv,2,3,Time,cpu,UGLMAT HYPRE,m-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,CPU Time (dancing\_eddies),Simulated Time (s),CPU Time (s),0,2,1,0,400,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_tight_precon,N/A,end,0,Pressure Solver,kd,k,TeX
f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,cpu,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_ulmat_hypre_devc.csv,2,3,Time,cpu,ULMAT HYPRE,g-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,CPU Time (dancing\_eddies),Simulated Time (s),CPU Time (s),0,2,1,0,400,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_tight_precon,N/A,end,0,Pressure Solver,kd,k,TeX
d,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_1mesh_devc.csv,2,3,Time,pres,"FDS, 1 mesh",k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,pres,"FDS, 4 mesh, tol=1E-5 m/s",r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure (dancing\_eddies\_tight),Time (s),Pressure (Pa),0,2,1,-0.3,0.2,1,no,0.05 0.90,SouthWest,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_tight,Absolute Error,end,0.3,Pressure Solver,kd,k,TeX
d,dancing_eddies,Pressure_Solver/dancing_eddies_default_git.txt,Pressure_Solver/dancing_eddies_1mesh_devc.csv,2,3,Time,pres,"FDS, 1 mesh",k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_default_devc.csv,2,3,Time,pres,"FDS, 4 mesh, tol=5E-4 m/s",r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure (dancing\_eddies\_default),Time (s),Pressure (Pa),0,2,1,-0.3,0.2,1,no,0.05 0.90,SouthWest,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_default,Absolute Error,end,0.01,Pressure Solver,kd,k,TeX
d,dancing_eddies,Pressure_Solver/dancing_eddies_uglmat_git.txt,Pressure_Solver/dancing_eddies_1mesh_devc.csv,2,3,Time,pres,"FDS, 1 mesh",k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_uglmat_devc.csv,2,3,Time,pres,"FDS, 4 mesh, uglmat",r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure (dancing\_eddies\_uglmat),Time (s),Pressure (Pa),0,2,1,-0.3,0.2,1,no,0.05 0.90,SouthWest,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_uglmat,Absolute Error,end,0.01,Pressure Solver,kd,k,TeX
f,dancing_eddies,Pressure_Solver/dancing_eddies_uglmat_git.txt,Pressure_Solver/dancing_eddies_1mesh_devc.csv,2,3,Time,pres,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_ulmat_devc.csv,2,3,Time,pres,"FDS, 4 mesh, ulmat",b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure (dancing\_eddies\_uglmat),Time (s),Pressure (Pa),0,2,1,-0.3,0.2,1,no,0.05 0.90,SouthWest,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_uglmat,Absolute Error,end,0.01,Pressure Solver,kd,k,TeX
f,dancing_eddies,Pressure_Solver/dancing_eddies_uglmat_git.txt,Pressure_Solver/dancing_eddies_1mesh_devc.csv,2,3,Time,pres,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_uglmat_hypre_devc.csv,2,3,Time,pres,"FDS, 4 mesh, uglmat hypre",m-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure (dancing\_eddies\_uglmat),Time (s),Pressure (Pa),0,2,1,-0.3,0.2,1,no,0.05 0.90,SouthWest,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_uglmat,Absolute Error,end,0.01,Pressure Solver,kd,k,TeX
f,dancing_eddies,Pressure_Solver/dancing_eddies_uglmat_git.txt,Pressure_Solver/dancing_eddies_1mesh_devc.csv,2,3,Time,pres,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_ulmat_hypre_devc.csv,2,3,Time,pres,"FDS, 4 mesh, ulmat hypre",g-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure (dancing\_eddies\_uglmat),Time (s),Pressure (Pa),0,2,1,-0.3,0.2,1,no,0.05 0.90,SouthWest,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_uglmat,Absolute Error,end,0.01,Pressure Solver,kd,k,TeX
d,device_restart,Restart/device_restart_base_case_git.txt,Restart/device_restart_base_case_devc.csv,2,3,Time,temp,Uninterrupted Case,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Restart/device_restart_a_devc.csv,2,3,Time,temp,Interrupted Case,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (device\_restart),Time (s),Temperature (°C),0,60,1,0,100,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/device_restart,Absolute Error,end,5,Restart,kd,k,TeX
d,devc_interpolation,Miscellaneous/devc_interpolation_git.txt,Miscellaneous/devc_interpolation.csv,1,2,z,u,{\it u} prescribed,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Miscellaneous/devc_interpolation_velocity_line.csv,2,3,Z,U-VEL-1|U-VEL-2|U-VEL-3,U-VEL-1|U-VEL-2|U-VEL-3,r--|b--|g--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Velocity (devc\_interpolation),{\it z} (m),Velocity (m/s),0,1,1,0,1,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/devc_interpolation_velocity,Absolute Error,end_1_3,0.01,Miscellaneous,kd,k,TeX
Expand All @@ -144,8 +147,10 @@ d,droplet_distributions,Sprinklers_and_Sprays/droplet_distributions_2_git.txt,Sp
d,droplet_distributions,Sprinklers_and_Sprays/droplet_distributions_2_git.txt,Sprinklers_and_Sprays/droplet_distributions_2.csv,1,2,Time,Mass,Total Mass (Prescibed),ko,0,40000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/droplet_distributions_2_devc.csv,2,3,Time,USER_M,PDPA Mass (Sampled),k-,0,4000,,0,100000,-1.00E+09,1.00E+09,0,User Defined Distribution Total Mass,Time (s),Mass (kg),0,11,1,0,0.1,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/droplet_distributions_2_4_mass,Relative Error,end,2.00E-02,Sprinklers and Sprays,ro,r,TeX
d,duct_flow,Pressure_Solver/duct_flow_git.txt,Pressure_Solver/duct_flow.csv,1,2,Time,Flow,Ideal (Flow),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/duct_flow_devc.csv,2,3,Time,flow_in|flow_out,FFT IBM (flow\_in)|FFT IBM (flow\_out),k-|r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Volume Flow (duct\_flow),Time (s),Volume Flow (m³/s),0,60,1,0,1.5,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/duct_flow,Relative Error,end_1_2,0.05,Pressure Solver,kd,k,TeX
f,duct_flow,Pressure_Solver/duct_flow_uglmat_git.txt,Pressure_Solver/duct_flow.csv,1,2,Time,Flow,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/duct_flow_uglmat_devc.csv,2,3,Time,flow_in|flow_out,UGLMAT (flow\_in)|UGLMAT (flow\_out),k-|b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Volume Flow (duct\_flow),Time (s),Volume Flow (m³/s),0,60,1,0,1.5,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/duct_flow,Relative Error,end_1_2,0.05,Pressure Solver,kd,k,TeX
f,duct_flow,Pressure_Solver/duct_flow_uglmat_hypre_git.txt,Pressure_Solver/duct_flow.csv,1,2,Time,Flow,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/duct_flow_uglmat_hypre_devc.csv,2,3,Time,flow_in|flow_out,UGLMAT HYPRE (flow\_in)|UGLMAT HYPRE (flow\_out),k-|m-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Volume Flow (duct\_flow),Time (s),Volume Flow (m³/s),0,60,1,0,1.5,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/duct_flow,Relative Error,end_1_2,0.05,Pressure Solver,kd,k,TeX
d,duct_flow,Pressure_Solver/duct_flow_git.txt,Pressure_Solver/duct_flow.csv,1,2,Time,Iterations,Ideal (Iterations),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/duct_flow_devc.csv,2,3,Time,iter,FFT IBM (iter),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (duct\_flow),Time (s),Iterations,0,60,1,0,100,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/duct_flow_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
f,duct_flow,Pressure_Solver/duct_flow_uglmat_git.txt,Pressure_Solver/duct_flow.csv,1,2,Time,Iterations,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/duct_flow_uglmat_devc.csv,2,3,Time,iter,UGLMAT (iter),b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (duct\_flow),Time (s),Iterations,0,60,1,0,100,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/duct_flow_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
f,duct_flow,Pressure_Solver/duct_flow_uglmat_hypre_git.txt,Pressure_Solver/duct_flow.csv,1,2,Time,Iterations,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/duct_flow_uglmat_hypre_devc.csv,2,3,Time,iter,UGLMAT HYPRE (iter),m-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (duct\_flow),Time (s),Iterations,0,60,1,0,100,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/duct_flow_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
d,e_coefficient,Sprinklers_and_Sprays/e_coefficient_git.txt,Sprinklers_and_Sprays/e_coefficient.csv,1,1,Time,HRR,Expected,ko,0,100000,,14.9,15.1,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/e_coefficient_devc.csv,2,3,Time,HRR_1,FDS,k-,0,100000,,14.9,15.1,-1.00E+09,1.00E+09,0,Heat Release Rate (OBST),Time (s),Heat Release Rate (kW),0,30,1,0,60,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/e_coefficient_obst,Absolute Error,max,3,Sprinklers and Sprays,bs,b,TeX
d,e_coefficient,Sprinklers_and_Sprays/e_coefficient_git.txt,Sprinklers_and_Sprays/e_coefficient.csv,1,1,Time,HRR,Expected,ko,0,100000,,14.9,15.1,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/e_coefficient_devc.csv,2,3,Time,HRR_2,FDS,k-,0,100000,,14.9,15.1,-1.00E+09,1.00E+09,0,Heat Release Rate (GEOM),Time (s),Heat Release Rate (kW),0,30,1,0,60,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/e_coefficient_geom,Absolute Error,max,3,Sprinklers and Sprays,bs,b,TeX
d,emissivity,Radiation/emissivity_git.txt,Radiation/emissivity.csv,1,2,Time,Emissivity,Analytical (Emissivity),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Radiation/emissivity_devc.csv,2,3,Time,EMISSIVITY,FDS (EMISSIVITY),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Surface Emissivity (emissivity),Time (s),Emissivity,0,10,1,0,1,1,no,0.30 0.90,SouthWest,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/emissivity,Absolute Error,end,0.1,Radiation,kd,k,TeX
Expand Down
1 change: 1 addition & 0 deletions Verification/FDS_Cases.sh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -402,6 +402,7 @@ $QFDS -p 8 -d Pressure_Solver duct_flow.fds
$QFDS -p 8 -d Pressure_Solver duct_flow_ulmat.fds
$QFDS -p 8 -d Pressure_Solver duct_flow_ulmat_hypre.fds
$QFDS -p 8 -d Pressure_Solver duct_flow_uglmat.fds
$QFDS -p 8 -d Pressure_Solver duct_flow_uglmat_hypre.fds
$QFDS -p 5 -d Pressure_Solver hallways.fds
$QFDS -p 8 -d Pressure_Solver tunnel_demo.fds
$QFDS -p 8 -d Pressure_Solver tunnel_demo_glmat.fds
Expand Down
86 changes: 86 additions & 0 deletions Verification/Pressure_Solver/duct_flow_uglmat_hypre.fds
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,86 @@
Test pressure iteration routine for flow through a duct across meshes

&HEAD CHID='duct_flow_uglmat_hypre', TITLE='Flow through a duct, test of UGLMAT with HYPRE' /

&MESH IJK=16,16,16, XB=-.2,3.0,-.2,3.0,-.2,3.0, MULT_ID='mesh' /
&MULT ID='mesh', DX=3.2, DY=3.2, DZ=3.2, I_UPPER=1, J_UPPER=1, K_UPPER=1 /

MESH IJK=32,32,32, XB=-.2,6.2,-.2,6.2,-.2,6.2 /

&RADI RADIATION=.FALSE. /
&TIME T_END=60.0 /

&MISC STRATIFICATION=.FALSE. /
&PRES CHECK_POISSON=T, SOLVER='UGLMAT HYPRE' /

&SURF ID='DUCT', COLOR='GRAY 30', TRANSPARENCY=0.5 /

&OBST XB=1.0,1.1,1.0,2.0,1.0,2.0, SURF_ID='DUCT', COLOR='RED',TRANSPARENCY=0.5 /
&VENT XB=1.0,1.0,1.0,2.0,1.0,2.0, SURF_ID='HVAC', ID='SUCK', COLOR='RED',TRANSPARENCY=0.5/
&VENT XB=1.1,1.1,1.0,2.0,1.0,2.0, SURF_ID='HVAC', ID='BLOW', COLOR='RED',TRANSPARENCY=0.5/
&HVAC TYPE_ID='NODE',ID='SUCK',VENT_ID='SUCK',DUCT_ID='FAN'/
&HVAC TYPE_ID='NODE',ID='BLOW',VENT_ID='BLOW',DUCT_ID='FAN'/
&HVAC TYPE_ID='DUCT',ID='FAN',NODE_ID='SUCK','BLOW',VOLUME_FLOW=1.,AREA=1./

&OBST XB=1.1,5.0,1.0,1.0,1.0,2.0, SURF_ID='DUCT' /
&OBST XB=1.1,4.0,2.0,2.0,1.0,2.0, SURF_ID='DUCT' /
&OBST XB=1.1,5.0,1.0,2.0,1.0,1.0, SURF_ID='DUCT' /
&OBST XB=1.1,5.0,1.0,2.0,2.0,2.0, SURF_ID='DUCT' /

&OBST XB=5.0,5.0,1.0,5.0,1.0,2.0, SURF_ID='DUCT' /
&OBST XB=4.0,4.0,2.0,4.0,1.0,2.0, SURF_ID='DUCT' /
&OBST XB=4.0,5.0,2.0,5.0,1.0,1.0, SURF_ID='DUCT' /
&OBST XB=4.0,5.0,2.0,5.0,2.0,2.0, SURF_ID='DUCT' /

&OBST XB=1.0,4.0,4.0,4.0,1.0,2.0, SURF_ID='DUCT' /
&OBST XB=1.0,5.0,5.0,5.0,1.0,2.0, SURF_ID='DUCT' /
&OBST XB=1.0,4.0,4.0,5.0,1.0,1.0, SURF_ID='DUCT' /
&OBST XB=2.0,4.0,4.0,5.0,2.0,2.0, SURF_ID='DUCT' /

&OBST XB=1.0,1.0,4.0,5.0,1.0,5.0, SURF_ID='DUCT' /
&OBST XB=2.0,2.0,4.0,5.0,2.0,4.0, SURF_ID='DUCT' /
&OBST XB=1.0,2.0,4.0,4.0,2.0,4.0, SURF_ID='DUCT' /
&OBST XB=1.0,2.0,5.0,5.0,2.0,4.0, SURF_ID='DUCT' /

&OBST XB=1.0,5.0,1.0,1.0,4.0,5.0, SURF_ID='DUCT' /
&OBST XB=1.0,4.0,2.0,2.0,4.0,5.0, SURF_ID='DUCT' /
&OBST XB=1.0,5.0,1.0,2.0,4.0,4.0, SURF_ID='DUCT' /
&OBST XB=1.0,5.0,1.0,2.0,5.0,5.0, SURF_ID='DUCT' /

&OBST XB=5.0,5.0,1.0,5.0,4.0,5.0, SURF_ID='DUCT' /
&OBST XB=4.0,4.0,2.0,4.0,4.0,5.0, SURF_ID='DUCT' /
&OBST XB=4.0,5.0,2.0,5.0,4.0,4.0, SURF_ID='DUCT' /
&OBST XB=4.0,5.0,2.0,5.0,5.0,5.0, SURF_ID='DUCT' /

&OBST XB=1.0,4.0,4.0,4.0,4.0,5.0, SURF_ID='DUCT' /
&OBST XB=1.0,5.0,5.0,5.0,4.0,5.0, SURF_ID='DUCT' /
&OBST XB=2.0,4.0,4.0,5.0,4.0,4.0, SURF_ID='DUCT' /
&OBST XB=1.0,4.0,4.0,5.0,5.0,5.0, SURF_ID='DUCT' /

&VENT MB='XMIN', SURF_ID='OPEN' /
&VENT MB='XMAX', SURF_ID='OPEN' /
&VENT MB='YMIN', SURF_ID='OPEN' /
&VENT MB='YMAX', SURF_ID='OPEN' /
&VENT MB='ZMIN', SURF_ID='OPEN' /
&VENT MB='ZMAX', SURF_ID='OPEN' /

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

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

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

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

&BNDF QUANTITY='VELOCITY ERROR' /

&DEVC XB=1.0,1.0,1.0,2.0,1.0,2.0, QUANTITY='U-VELOCITY', SPATIAL_STATISTIC='AREA INTEGRAL', QUANTITY_RANGE(1)=0, ID='flow_in' /
&DEVC XB=1.0,1.0,1.0,2.0,4.0,5.0, QUANTITY='U-VELOCITY', SPATIAL_STATISTIC='AREA INTEGRAL', QUANTITY_RANGE(2)=0, ID='flow_out', CONVERSION_FACTOR=-1. /
&DEVC XYZ=0.,0.,0., QUANTITY='MAXIMUM VELOCITY ERROR', ID='error' /
&DEVC XYZ=0.,0.,0., QUANTITY='PRESSURE ITERATIONS', ID='iter' /

&TAIL /
Loading