diff --git a/Source/wall.f90 b/Source/wall.f90 index 8991b810ff7..c9dd3711b4f 100644 --- a/Source/wall.f90 +++ b/Source/wall.f90 @@ -164,10 +164,9 @@ SUBROUTINE WALL_BC(T,DT,NM) IF (WC%BOUNDARY_TYPE/=OPEN_BOUNDARY .AND. WC%BOUNDARY_TYPE/=INTERPOLATED_BOUNDARY) THEN CALL CALCULATE_ZZ_F(T,DT,WALL_INDEX=IW) + CALL CALCULATE_RHO_F(BC,B1,WALL_INDEX=IW) ENDIF - IF (WC%BOUNDARY_TYPE/=INTERPOLATED_BOUNDARY) CALL CALCULATE_RHO_F(BC,B1,WALL_INDEX=IW) - ENDDO WALL_CELL_LOOP !$OMP END DO @@ -220,7 +219,6 @@ SUBROUTINE WALL_BC(T,DT,NM) ENDIF CALL CALCULATE_ZZ_F(T,DT,CFACE_INDEX=ICF) - CALL CALCULATE_RHO_F(BC,B1,CFACE_INDEX=ICF) ENDDO CFACE_LOOP @@ -1328,9 +1326,9 @@ SUBROUTINE CALCULATE_ZZ_F(T,DT,WALL_INDEX,CFACE_INDEX,PARTICLE_INDEX) ! Only set species mass fraction in the ghost cell if it is not solid IF (PRESENT(WALL_INDEX)) THEN - IF (WALL_INDEX<=N_EXTERNAL_WALL_CELLS .AND. .NOT.CELL(IC)%SOLID .AND. .NOT.CELL(ICG)%SOLID) THEN - ZZP(BC%II,BC%JJ,BC%KK,1:N_TRACKED_SPECIES) = 2._EB*B1%ZZ_F(1:N_TRACKED_SPECIES) - B1%ZZ_G(1:N_TRACKED_SPECIES) - ZZP(BC%II2,BC%JJ2,BC%KK2,1:N_TRACKED_SPECIES) = 2._EB*B1%ZZ_F(1:N_TRACKED_SPECIES) - B1%ZZ_G(1:N_TRACKED_SPECIES) + IF (WALL_INDEX<=N_EXTERNAL_WALL_CELLS) THEN + ZZP(BC%II,BC%JJ,BC%KK,1:N_TRACKED_SPECIES) = B1%ZZ_F(1:N_TRACKED_SPECIES) + ZZP(BC%II2,BC%JJ2,BC%KK2,1:N_TRACKED_SPECIES) = B1%ZZ_F(1:N_TRACKED_SPECIES) ENDIF ENDIF @@ -1496,10 +1494,14 @@ SUBROUTINE CALCULATE_RHO_F(BC,B1,WALL_INDEX,CFACE_INDEX) ENDIF ENDIF +! Fill exterior ghost cells with surface density and temperature + IF (PRESENT(WALL_INDEX)) THEN - IF (WALL_INDEX<=N_EXTERNAL_WALL_CELLS .AND. BOUNDARY_TYPE/=OPEN_BOUNDARY) THEN - RHOP(BC%II ,BC%JJ, BC%KK ) = 2._EB*B1%RHO_F - B1%RHO_G - RHOP(BC%II2,BC%JJ2,BC%KK2) = 2._EB*B1%RHO_F - B1%RHO_G + IF (WALL_INDEX<=N_EXTERNAL_WALL_CELLS) THEN + RHOP(BC%II ,BC%JJ, BC%KK ) = B1%RHO_F + RHOP(BC%II2,BC%JJ2,BC%KK2) = B1%RHO_F + TMP(BC%II ,BC%JJ, BC%KK ) = B1%TMP_F + TMP(BC%II2,BC%JJ2,BC%KK2) = B1%TMP_F ENDIF ENDIF diff --git a/Utilities/Matlab/FDS_verification_dataplot_inputs.csv b/Utilities/Matlab/FDS_verification_dataplot_inputs.csv index cbe2fef8de3..72cc5640855 100644 --- a/Utilities/Matlab/FDS_verification_dataplot_inputs.csv +++ b/Utilities/Matlab/FDS_verification_dataplot_inputs.csv @@ -108,10 +108,10 @@ d,condensation_evaporation,Species/wall_cond_git.txt,Species/wall_cond.csv,1,2,T d,condensation_evaporation,Species/wall_cond_git.txt,Species/wall_cond.csv,1,2,Time,dP,Analytical,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Species/wall_cond_devc.csv,2,3,Time,pres,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure (wall\_cond),Time (s),Pressure (Pa),0,20000,1,-25000,0,1,no,0.05 0.80,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/wall_cond_P.pdf,Relative Error,end,0.05,Condensation,r^,r,TeX d,control_test_2,Controls/control_test_2_git.txt,Controls/control_test_2.csv,1,2,Time,Add|Multiply|Subtract|Divide|Power|PID,Add|Multiply|Subtract|Divide|Power|PID,ko|ro|bo|go|mo|co,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Controls/control_test_2_devc.csv,2,3,Time,CF Add|CF Multiply|CF Subtract|CF Divide|CF Power|CF PID,FDS Add|FDS Multiply|FDS Subtract|FDS Divide|FDS Power|FDS PID,k-|r-|b-|g-|m-|c-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Control Function Outputs (control\_test\_2),Time (s),Control Function Output,0,10,1,-15,40,1,no,0.05 0.90,EastOutside,,1.35,linear,FDS_User_Guide/SCRIPT_FIGURES/control_test_2,Relative Error,end,0.01,Controls,kd,k,TeX d,convective_cooling,Heat_Transfer/convective_cooling_git.txt,Heat_Transfer/convective_cooling.csv,1,2,Time,Temp,Exact (Temp),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/convective_cooling_devc.csv,2,3,Time,inner temp,FDS (inner temp),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Inner Temperature (convective\_cooling),Time (min),Temperature (°C),0,30,60,0,1200,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/convective_cooling,Relative Error,end,0.01,Heat Transfer,r^,r,TeX -d,convective_cooling,Heat_Transfer/convective_cooling_p1_git.txt,Heat_Transfer/convective_cooling_error.csv,1,2,dx,dx|dx^2,O(?x)|O(?x²),k--|k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/convective_cooling_error.csv,1,2,dx,relative error,FDS (Relative Error),ko-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Convergence (convective\_cooling),Grid Spacing (m),Relative Error,1.00E-05,1.00E-03,1,1.00E-08,1.00E-02,1,no,0.03 0.90,SouthEast,,1,loglog,FDS_Verification_Guide/SCRIPT_FIGURES/convective_cooling_error,N/A,end,0,Heat Transfer,kd,k,TeX +d,convective_cooling,Heat_Transfer/convective_cooling_p1_git.txt,Heat_Transfer/convective_cooling_error.csv,1,2,dx,dx|dx^2,O(\deltax)|O(\deltax²),k--|k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/convective_cooling_error.csv,1,2,dx,relative error,FDS (Relative Error),ko-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Convergence (convective\_cooling),Grid Spacing (m),Relative Error,1.00E-05,1.00E-03,1,1.00E-08,1.00E-02,1,no,0.03 0.90,SouthEast,,1,loglog,FDS_Verification_Guide/SCRIPT_FIGURES/convective_cooling_error,N/A,end,0,Heat Transfer,kd,k,TeX d,couch,Fires/couch_git.txt,Fires/couch.csv,1,2,Time,HRR,Expected Peak,k--,0,100000,,100,300,-1.00E+09,1.00E+09,0,Fires/couch_hrr.csv,2,3,Time,HRR,FDS (HRR),k-,0,100000,,100,300,-1.00E+09,1.00E+09,0,Heat Release Rate (couch),Time (min),Heat Release Rate (kW),0,10,60,0,4000,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/couch,Relative Error,max,0.2,Fires,r*,r,TeX -d,cyl_test,Flowfields/cyl_test_1_git.txt,Flowfields/cyl_test_1_devc.csv,2,3,Time,temp1,?y=2 mm,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Flowfields/cyl_test_2_devc.csv,2,3,Time,temp2,?y=2 cm,r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,"Cylindrical Test, LES (cyl\_test\_1-2)",Time (s),Temperature (°C),0,20,1,0,100,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/cyl_test,Relative Error,end,0.01,Flowfields,r^,r,TeX -d,cyl_test,Flowfields/cyl_test_3_git.txt,Flowfields/cyl_test_3_devc.csv,2,3,Time,temp1,?y=2 mm,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Flowfields/cyl_test_4_devc.csv,2,3,Time,temp1,?y=4 cm,r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,"Cylindrical Test, DNS (cyl\_test\_3-4)",Time (s),Temperature (°C),0,3,1,0,2000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/cyl_test_3-4,Relative Error,end,0.01,Flowfields,r^,r,TeX +d,cyl_test,Flowfields/cyl_test_1_git.txt,Flowfields/cyl_test_1_devc.csv,2,3,Time,temp1,\deltay=2 mm,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Flowfields/cyl_test_2_devc.csv,2,3,Time,temp2,\deltay=2 cm,r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,"Cylindrical Test, LES (cyl\_test\_1-2)",Time (s),Temperature (°C),0,20,1,0,100,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/cyl_test,Relative Error,end,0.01,Flowfields,r^,r,TeX +d,cyl_test,Flowfields/cyl_test_3_git.txt,Flowfields/cyl_test_3_devc.csv,2,3,Time,temp1,\deltay=2 mm,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Flowfields/cyl_test_4_devc.csv,2,3,Time,temp1,\deltay=4 cm,r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,"Cylindrical Test, DNS (cyl\_test\_3-4)",Time (s),Temperature (°C),0,3,1,0,2000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/cyl_test_3-4,Relative Error,end,0.01,Flowfields,r^,r,TeX 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 @@ -510,8 +510,8 @@ d,Needle_TGA,WUI/Needle_TGA_E_cat_git.txt,WUI/USDA_FPL_Longleaf_Pine_Needle_TGA_ d,Needle_TGA,WUI/Needle_TGA_F_cat_git.txt,WUI/USDA_FPL_Longleaf_Pine_Needle_TGA_Data.csv,2,3,F-Time,F-Temperature,Exp (Temperature),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,WUI/Needle_TGA_F_cat_devc.csv,2,3,Time,T,FDS (T),k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Case F; 60 °C/min; M=0.064; Inert,Time (min),Temperature (°C),0,20,60,0,800,1,no,0.05 0.90,East,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/Needle_TGA_F_Time,N/A,end,0,Needle TGA,kd,k,TeX d,Needle_TGA,WUI/Needle_TGA_G_cat_git.txt,WUI/USDA_FPL_Longleaf_Pine_Needle_TGA_Data.csv,2,3,G-Time,G-Temperature,Exp (Temperature),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,WUI/Needle_TGA_G_cat_devc.csv,2,3,Time,T,FDS (T),k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Case G; 60 °C/min; M=1.56; Air,Time (min),Temperature (°C),0,40,60,0,800,1,no,0.05 0.90,East,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/Needle_TGA_G_Time,N/A,end,0,Needle TGA,kd,k,TeX d,Needle_TGA,WUI/Needle_TGA_H_cat_git.txt,WUI/USDA_FPL_Longleaf_Pine_Needle_TGA_Data.csv,2,3,H-Time,H-Temperature,Exp (Temperature),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,WUI/Needle_TGA_H_cat_devc.csv,2,3,Time,T,FDS (T),k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Case H; 60 °C/min; M=1.33; Air,Time (min),Temperature (°C),0,30,60,0,800,1,no,0.05 0.90,East,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/Needle_TGA_H_Time,N/A,end,0,Needle TGA,kd,k,TeX -d,ns2d,NS_Analytical_Solution/ns2d_8_git.txt,NS_Analytical_Solution/ns2d_error.csv,1,2,dx,dx|dx^2,O(?x)|O(?x²),k--|k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,NS_Analytical_Solution/ns2d_error.csv,1,2,dx,rms error,FDS (RMS Error),k*-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,"Convergence, Inviscid Case (ns2d)",Grid Spacing (m),RMS Error (m/s),0.09,0.9,1,0.001,1,1,no,0.03 0.90,SouthEast,,1,loglog,FDS_Verification_Guide/SCRIPT_FIGURES/ns2d_error,N/A,end,0,NS Analytical Solution,kd,k,TeX -d,ns2d_nupt1,NS_Analytical_Solution/ns2d_8_nupt1_git.txt,NS_Analytical_Solution/ns2d_nupt1_error.csv,1,2,dx,dx|dx^2,O(?x)|O(?x²),k--|k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,NS_Analytical_Solution/ns2d_nupt1_error.csv,1,2,dx,rms error,FDS (RMS Error),k*-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,"Convergence, Viscous Case (ns2d)",Grid Spacing (m),RMS Error (m/s),0.09,0.9,1,0.001,1,1,no,0.03 0.90,SouthEast,,1,loglog,FDS_Verification_Guide/SCRIPT_FIGURES/ns2d_nupt1_error,N/A,end,0,NS Analytical Solution,kd,k,TeX +d,ns2d,NS_Analytical_Solution/ns2d_8_git.txt,NS_Analytical_Solution/ns2d_error.csv,1,2,dx,dx|dx^2,O(\deltax)|O(\deltax²),k--|k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,NS_Analytical_Solution/ns2d_error.csv,1,2,dx,rms error,FDS (RMS Error),k*-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,"Convergence, Inviscid Case (ns2d)",Grid Spacing (m),RMS Error (m/s),0.09,0.9,1,0.001,1,1,no,0.03 0.90,SouthEast,,1,loglog,FDS_Verification_Guide/SCRIPT_FIGURES/ns2d_error,N/A,end,0,NS Analytical Solution,kd,k,TeX +d,ns2d_nupt1,NS_Analytical_Solution/ns2d_8_nupt1_git.txt,NS_Analytical_Solution/ns2d_nupt1_error.csv,1,2,dx,dx|dx^2,O(\deltax)|O(\deltax²),k--|k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,NS_Analytical_Solution/ns2d_nupt1_error.csv,1,2,dx,rms error,FDS (RMS Error),k*-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,"Convergence, Viscous Case (ns2d)",Grid Spacing (m),RMS Error (m/s),0.09,0.9,1,0.001,1,1,no,0.03 0.90,SouthEast,,1,loglog,FDS_Verification_Guide/SCRIPT_FIGURES/ns2d_nupt1_error,N/A,end,0,NS Analytical Solution,kd,k,TeX d,ns2d_8_nupt1,NS_Analytical_Solution/ns2d_8_nupt1_git.txt,NS_Analytical_Solution/ns2d_8_nupt1_exact.csv,1,2,Time,u-vel,Analytical (u-vel),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,NS_Analytical_Solution/ns2d_8_nupt1_devc.csv,2,3,Time,UVEL,FDS (UVEL),k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Velocity (ns2d\_8\_nupt1),Time (s),Velocity (m/s),0,7,1,0.3,2.4,1,no,0.05 0.90,NorthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ns2d_8_nupt1,Convergent Series,end,0,NS Analytical Solution,kd,k,TeX d,ns2d_16_nupt1,NS_Analytical_Solution/ns2d_16_nupt1_git.txt,NS_Analytical_Solution/ns2d_16_nupt1_exact.csv,1,2,Time,u-vel,Analytical (u-vel),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,NS_Analytical_Solution/ns2d_16_nupt1_devc.csv,2,3,Time,UVEL,FDS (UVEL),k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Velocity (ns2d\_16\_nupt1),Time (s),Velocity (m/s),0,7,1,0.3,2.4,1,no,0.05 0.90,NorthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ns2d_16_nupt1,Convergent Series,end,0,NS Analytical Solution,kd,k,TeX d,ns2d_32_nupt1,NS_Analytical_Solution/ns2d_32_nupt1_git.txt,NS_Analytical_Solution/ns2d_32_nupt1_exact.csv,1,2,Time,u-vel,Analytical (u-vel),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,NS_Analytical_Solution/ns2d_32_nupt1_devc.csv,2,3,Time,UVEL,FDS (UVEL),k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Velocity (ns2d\_32\_nupt1),Time (s),Velocity (m/s),0,7,1,0.3,2.4,1,no,0.05 0.90,NorthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ns2d_32_nupt1,Convergent Series,end,0,NS Analytical Solution,kd,k,TeX @@ -728,7 +728,7 @@ d,velocity_bc_test,Flowfields/velocity_bc_test_git.txt,Flowfields/velocity_bc_te d,velocity_bc_test,Flowfields/velocity_bc_test_git.txt,Flowfields/velocity_bc_test_line.csv,2,3,x1p,Pres1,Thin Plate (Pres1),k-,-6,21,,-6,21,-1.00E+09,1.00E+09,0,Flowfields/velocity_bc_test_line.csv,2,3,x1p,Pres2|Pres3,Solid Block (Pres2)|Thick Plate (Pres3),r-|g-,-6,21,,-6,21,-1.00E+09,1.00E+09,0,Centerline Pressure (velocity\_bc\_test),Distance (m),Pressure (Pa),-5,20,1,-0.5,1,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/velocity_bc_test_pres,Absolute Error,mean_1_2,0.05,Flowfields,kd,k,TeX d,volume_flow_1,Flowfields/volume_flow_1_git.txt,Flowfields/volume_flow.csv,1,2,Time,vflow,Expected,ko,0,100000,,50,100,-1.00E+09,1.00E+09,0,Flowfields/volume_flow_1_devc.csv,2,3,Time,vflow,FDS,k-,0,100000,,50,100,-1.00E+09,1.00E+09,0,Flow Test (volume\_flow\_1),Time (s),Volume Flow (m³/s),0,100,1,0,0.015,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/volume_flow_1,Relative Error,mean,0.01,Flowfields,kd,k,TeX d,volume_flow_2,Flowfields/volume_flow_2_git.txt,Flowfields/volume_flow.csv,1,2,Time,vflow,Expected,ko,0,100000,,50,100,-1.00E+09,1.00E+09,0,Flowfields/volume_flow_2_devc.csv,2,3,Time,vflow,FDS,k-,0,100000,,50,100,-1.00E+09,1.00E+09,0,Flow Test (volume\_flow\_2),Time (s),Volume Flow (m³/s),0,100,1,0,0.015,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/volume_flow_2,Relative Error,mean,0.01,Flowfields,kd,k,TeX -d,vort2d,NS_Analytical_Solution/vort2d_40_git.txt,NS_Analytical_Solution/vort2d_error.csv,1,2,dx,dxmod|dx^2,O(?x)|O(?x²),k--|k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,NS_Analytical_Solution/vort2d_error.csv,1,2,dx,rms1|rms2|rms3,FDS (RMS error at t=0.0089 s)|FDS (RMS error at t=0.0178 s)|FDS (RMS error at t=0.0267 s),g*-|co-|r^-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Convergence (vort2d),Grid Spacing (m),RMS Error (m/s),0.0009,0.009,1,0.001,1,1,no,0.03 0.90,SouthEast,,1,loglog,FDS_Verification_Guide/SCRIPT_FIGURES/vort2d_error,N/A,end,0,NS Analytical Solution,kd,k,TeX +d,vort2d,NS_Analytical_Solution/vort2d_40_git.txt,NS_Analytical_Solution/vort2d_error.csv,1,2,dx,dxmod|dx^2,O(\deltax)|O(\deltax²),k--|k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,NS_Analytical_Solution/vort2d_error.csv,1,2,dx,rms1|rms2|rms3,FDS (RMS error at t=0.0089 s)|FDS (RMS error at t=0.0178 s)|FDS (RMS error at t=0.0267 s),g*-|co-|r^-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Convergence (vort2d),Grid Spacing (m),RMS Error (m/s),0.0009,0.009,1,0.001,1,1,no,0.03 0.90,SouthEast,,1,loglog,FDS_Verification_Guide/SCRIPT_FIGURES/vort2d_error,N/A,end,0,NS Analytical Solution,kd,k,TeX d,water_evap_1_const_gamma,Sprinklers_and_Sprays/water_evap_1_const_gamma_git.txt,Sprinklers_and_Sprays/water_evap_1_const_gamma.csv,1,2,Time,Rel. Hum,Exact (Rel. Hum),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/water_evap_1_const_gamma_devc.csv,2,3,Time,humid,FDS (humid),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Relative Humidity (water\_evap\_1\_const\_gamma),Time (s),Humidity (%),0,10,1,0,3,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/water_evap_1_const_gamma_humidity,Relative Error,end,0.01,Sprinklers and Sprays,kd,k,TeX d,water_evap_1_const_gamma,Sprinklers_and_Sprays/water_evap_1_const_gamma_git.txt,Sprinklers_and_Sprays/water_evap_1_const_gamma.csv,1,2,Time,h_gas|h_water,Expected (h\_gas)|Exact (h\_water),go|ro,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/water_evap_1_const_gamma_devc.csv,2,3,Time,h_gas|h_water,FDS (h\_gas)|FDS (h\_water),g-|r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Enthalpy Change (water\_evap\_1\_const\_gamma),Time (s),Enthalpy (kJ),0,10,1,-50,50,1,no,0.05 0.90,East,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/water_evap_1_const_gamma_enthalpy,Relative Error,end,0.01,Sprinklers and Sprays,kd,k,TeX d,water_evap_1_const_gamma,Sprinklers_and_Sprays/water_evap_1_const_gamma_git.txt,Sprinklers_and_Sprays/water_evap_1_const_gamma.csv,1,2,Time,dens,Exact (dens),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/water_evap_1_const_gamma_devc.csv,2,3,Time,dens,FDS (dens),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Density Change (water\_evap\_1\_const\_gamma),Time (s),Density (kg/m³),0,10,1,0,0.015,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/water_evap_1_const_gamma_density,Relative Error,end,0.01,Sprinklers and Sprays,kd,k,TeX diff --git a/Verification/Flowfields/helium_air_jet_floor.fds b/Verification/Flowfields/helium_air_jet_floor.fds index 14878d37e65..f9f84561a74 100644 --- a/Verification/Flowfields/helium_air_jet_floor.fds +++ b/Verification/Flowfields/helium_air_jet_floor.fds @@ -2,16 +2,24 @@ &MESH IJK=60,60,30, XB=0.0,1.50,0.0,1.50,0.0,0.75 / -&TIME T_BEGIN=-.1,T_END=2. / +&TIME T_BEGIN=0.0,T_END=2. / + +&MISC STRATIFICATION=F / + +&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='ZMAX', SURF_ID='OPEN' / &SPEC ID='HELIUM' / &SURF ID='helium', SPEC_ID(1)='HELIUM', MASS_FLUX(1)=0.040, RAMP_MF(1)='LEAKRAMP' / -&RAMP ID='LEAKRAMP',T=0.0,F=0.0 / +&RAMP ID='LEAKRAMP',T=0.0,F=1.0 / &RAMP ID='LEAKRAMP',T=1.0,F=1.0 / -&VENT XB=0.734,0.766,0.734,0.766,0.0,0.0, SURF_ID='helium', COLOR='RED' / +&VENT XB=0.725,0.775,0.725,0.775,0.0,0.0, SURF_ID='helium', COLOR='RED' / &SLCF PBY=0.75, QUANTITY='VELOCITY', VECTOR=.TRUE./ &SLCF PBY=0.75, QUANTITY='VOLUME FRACTION', SPEC_ID='HELIUM', CELL_CENTERED=.TRUE./