Skip to content

Commit d1d0afd

Browse files
authored
Merge pull request #13676 from ericvmueller/master
FDS Validation: add initial NIST USFS Camp Swift case: burn block 1 a…
2 parents a8fce54 + 3f0b3d9 commit d1d0afd

File tree

4 files changed

+530
-0
lines changed

4 files changed

+530
-0
lines changed
Lines changed: 202 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,202 @@
1+
&HEAD CHID='CS_BB1_0p25dx', TITLE='Camp Swift burn block 1, 25cm resolution, constant U magnitude and dir' /
2+
3+
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
4+
NOTE:
5+
Because FDS/SMV can only handle single precision in some cases we have adjusted the coordinates
6+
relative to a new origin in EPSG:26914, based on the position of anemometer 42
7+
E: 666,295.73 meters
8+
N: 3,347,184.94 meters
9+
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
10+
11+
400 meshes
12+
&MULT ID='xblock',DX=10,DY=10,I_LOWER=0,I_UPPER=19,J_LOWER=0,J_UPPER=19/
13+
&MESH IJK=40,40,40, XB=-65.73,-55.73,-34.94,-24.94,0,30,MULT_ID='xblock',TRNZ_ID='TRNZ' /
14+
15+
&TIME T_BEGIN=0,T_END=180./
16+
17+
&MISC TMPA=32./
18+
19+
-- Constant inflow
20+
&WIND SPEED=3.4, DIRECTION=312.5, L=-500., Z_0=0.03, Z_REF=3.3/
21+
22+
-- Boundary conditions
23+
&VENT MB='XMIN', SURF_ID='OPEN' /
24+
&VENT MB='XMAX', SURF_ID='OPEN' /
25+
&VENT MB='YMIN', SURF_ID='OPEN' /
26+
&VENT MB='YMAX', SURF_ID='OPEN' /
27+
&VENT MB='ZMAX', SURF_ID='OPEN' /
28+
29+
&SURF ID='fuel_break'
30+
RGB=115,76,0
31+
MATL_ID='DIRT'
32+
THICKNESS=0.2 /
33+
34+
&MATL ID='DIRT'
35+
CONDUCTIVITY = 0.25
36+
SPECIFIC_HEAT = 2.
37+
DENSITY = 1300. /
38+
39+
&SURF ID='GROUND'
40+
DEFAULT=T
41+
MATL_ID='DIRT'
42+
THICKNESS=0.2 /
43+
44+
-- vegetation REAC, MATL, SURF, PART properties, plus INIT lines using 25cm fuel model hosted in firemodels/cad
45+
&CATF OTHER_FILES='./vegetation_model.fds' /
46+
47+
-- Outputs
48+
&DUMP DT_HRR=0.1, DT_DEVC=0.1, DT_PART=2, DT_SLCF=1, RESULTS_DIR='smv_out'/
49+
50+
&SLCF PBZ=1.0, QUANTITY='VELOCITY', VECTOR=T /
51+
&SLCF PBZ=1.0, QUANTITY='HRRPUV'/
52+
&SLCF PBZ=1.0, QUANTITY='TEMPERATURE'/
53+
54+
-- Measurement devices in field
55+
&DEVC ID='9NFBP', XYZ=22.94, 105.36, 0.9, QUANTITY='RADIATIVE HEAT FLUX GAS',ORIENTATION=-0.292,0.956,0. /
56+
&DEVC ID='15WFBP', XYZ=4.18, 67.87, 0.9, QUANTITY='RADIATIVE HEAT FLUX GAS',ORIENTATION=-0.309,0.951,0. /
57+
&DEVC ID='6CFBP', XYZ=32.62, 75.89, 0.9, QUANTITY='RADIATIVE HEAT FLUX GAS',ORIENTATION=-0.375,0.927,0. /
58+
&DEVC ID='13EFBP', XYZ=65.84, 80.80, 0.9, QUANTITY='RADIATIVE HEAT FLUX GAS',ORIENTATION=-0.588,0.809,0. /
59+
&DEVC ID='17SFBP', XYZ=49.83, 52.58, 0.9, QUANTITY='RADIATIVE HEAT FLUX GAS',ORIENTATION=-0.208,0.978,0. /
60+
61+
-- Other outputs
62+
cDEVC ID='x', QUANTITY='TEMPERATURE', SPATIAL_STATISTIC='MAXLOC X', XB=-10,110,-1,0,0,1 /
63+
cDEVC ID='D', QUANTITY='FIRE DEPTH', XB=-10,110,0,0,0.1,0.1 /
64+
65+
&DEVC ID='mass_BS_lwrbld',XB=-66,135,-35,165,0,1, QUANTITY='MPUV', PART_ID='bluestem_lower_blade', STATISTICS='VOLUME INTEGRAL' /
66+
&DEVC ID='mass_BS_lwrcyl',XB=-66,135,-35,165,0,1, QUANTITY='MPUV', PART_ID='bluestem_lower_cyl', STATISTICS='VOLUME INTEGRAL' /
67+
&DEVC ID='mass_BS_uprlgcyl',XB=-66,135,-35,165,0,1, QUANTITY='MPUV', PART_ID='bluestem_upper_large_cyl', STATISTICS='VOLUME INTEGRAL' /
68+
&DEVC ID='mass_BS_uprsmcyl',XB=-66,135,-35,165,0,1, QUANTITY='MPUV', PART_ID='bluestem_upper_small_cyl', STATISTICS='VOLUME INTEGRAL' /
69+
&DEVC ID='mass_BS_uprbld',XB=-66,135,-35,165,0,1, QUANTITY='MPUV', PART_ID='bluestem_upper_blade', STATISTICS='VOLUME INTEGRAL' /
70+
71+
&DEVC ID='mass_CW_lgcyl',XB=-66,135,-35,165,0,1, QUANTITY='MPUV', PART_ID='camphorweed_large_cyl', STATISTICS='VOLUME INTEGRAL' /
72+
&DEVC ID='mass_CW_smcyl',XB=-66,135,-35,165,0,1, QUANTITY='MPUV', PART_ID='camphorweed_small_cyl', STATISTICS='VOLUME INTEGRAL' /
73+
&DEVC ID='mass_CW_leaf',XB=-66,135,-35,165,0,1, QUANTITY='MPUV', PART_ID='camphorweed_leaf', STATISTICS='VOLUME INTEGRAL' /
74+
75+
&DEVC ID='mass_TA_bld',XB=-66,135,-35,165,0,1, QUANTITY='MPUV', PART_ID='threeawn_blade', STATISTICS='VOLUME INTEGRAL' /
76+
&DEVC ID='mass_TA_cyl',XB=-66,135,-35,165,0,1, QUANTITY='MPUV', PART_ID='threeawn_cyl', STATISTICS='VOLUME INTEGRAL' /
77+
78+
-- Ignitor
79+
80+
&SURF ID='LINEFIRE_W1',HRRPUA=1000,RAMP_Q='RAMPIGN_W1',RGB=250,0,0, XYZ=0.73,116.56,0, SPREAD_RATE=1.314 /
81+
&RAMP ID='RAMPIGN_W1',T= 9,F=0 /
82+
&RAMP ID='RAMPIGN_W1',T=10,F=1 /
83+
&RAMP ID='RAMPIGN_W1',T=15,F=1 /
84+
&RAMP ID='RAMPIGN_W1',T=16,F=0 /
85+
86+
&VENT XB = -25, -24, 105, 106, 0, 0, SURF_ID = 'LINEFIRE_W1' /
87+
&VENT XB = -24, -23, 105, 106, 0, 0, SURF_ID = 'LINEFIRE_W1' /
88+
&VENT XB = -24, -23, 106, 107, 0, 0, SURF_ID = 'LINEFIRE_W1' /
89+
&VENT XB = -23, -22, 106, 107, 0, 0, SURF_ID = 'LINEFIRE_W1' /
90+
&VENT XB = -22, -21, 106, 107, 0, 0, SURF_ID = 'LINEFIRE_W1' /
91+
&VENT XB = -22, -21, 107, 108, 0, 0, SURF_ID = 'LINEFIRE_W1' /
92+
&VENT XB = -21, -20, 107, 108, 0, 0, SURF_ID = 'LINEFIRE_W1' /
93+
&VENT XB = -20, -19, 107, 108, 0, 0, SURF_ID = 'LINEFIRE_W1' /
94+
&VENT XB = -20, -19, 108, 109, 0, 0, SURF_ID = 'LINEFIRE_W1' /
95+
&VENT XB = -19, -18, 108, 109, 0, 0, SURF_ID = 'LINEFIRE_W1' /
96+
&VENT XB = -18, -17, 108, 109, 0, 0, SURF_ID = 'LINEFIRE_W1' /
97+
&VENT XB = -17, -16, 108, 109, 0, 0, SURF_ID = 'LINEFIRE_W1' /
98+
&VENT XB = -17, -16, 109, 110, 0, 0, SURF_ID = 'LINEFIRE_W1' /
99+
&VENT XB = -16, -15, 109, 110, 0, 0, SURF_ID = 'LINEFIRE_W1' /
100+
&VENT XB = -15, -14, 109, 110, 0, 0, SURF_ID = 'LINEFIRE_W1' /
101+
&VENT XB = -15, -14, 110, 111, 0, 0, SURF_ID = 'LINEFIRE_W1' /
102+
&VENT XB = -14, -13, 110, 111, 0, 0, SURF_ID = 'LINEFIRE_W1' /
103+
&VENT XB = -13, -12, 110, 111, 0, 0, SURF_ID = 'LINEFIRE_W1' /
104+
&VENT XB = -13, -12, 111, 112, 0, 0, SURF_ID = 'LINEFIRE_W1' /
105+
&VENT XB = -12, -11, 111, 112, 0, 0, SURF_ID = 'LINEFIRE_W1' /
106+
&VENT XB = -11, -10, 111, 112, 0, 0, SURF_ID = 'LINEFIRE_W1' /
107+
&VENT XB = -11, -10, 112, 113, 0, 0, SURF_ID = 'LINEFIRE_W1' /
108+
&VENT XB = -10, -9, 112, 113, 0, 0, SURF_ID = 'LINEFIRE_W1' /
109+
&VENT XB = -9, -8, 112, 113, 0, 0, SURF_ID = 'LINEFIRE_W1' /
110+
&VENT XB = -9, -8, 113, 114, 0, 0, SURF_ID = 'LINEFIRE_W1' /
111+
&VENT XB = -8, -7, 113, 114, 0, 0, SURF_ID = 'LINEFIRE_W1' /
112+
&VENT XB = -7, -6, 113, 114, 0, 0, SURF_ID = 'LINEFIRE_W1' /
113+
&VENT XB = -7, -6, 114, 115, 0, 0, SURF_ID = 'LINEFIRE_W1' /
114+
&VENT XB = -6, -5, 114, 115, 0, 0, SURF_ID = 'LINEFIRE_W1' /
115+
&VENT XB = -5, -4, 114, 115, 0, 0, SURF_ID = 'LINEFIRE_W1' /
116+
&VENT XB = -5, -4, 115, 116, 0, 0, SURF_ID = 'LINEFIRE_W1' /
117+
&VENT XB = -4, -3, 115, 116, 0, 0, SURF_ID = 'LINEFIRE_W1' /
118+
&VENT XB = -3, -2, 115, 116, 0, 0, SURF_ID = 'LINEFIRE_W1' /
119+
&VENT XB = -3, -2, 116, 117, 0, 0, SURF_ID = 'LINEFIRE_W1' /
120+
&VENT XB = -2, -1, 116, 117, 0, 0, SURF_ID = 'LINEFIRE_W1' /
121+
&VENT XB = -1, 0, 116, 117, 0, 0, SURF_ID = 'LINEFIRE_W1' /
122+
&VENT XB = -1, 0, 117, 118, 0, 0, SURF_ID = 'LINEFIRE_W1' /
123+
&VENT XB = 0, 1, 117, 118, 0, 0, SURF_ID = 'LINEFIRE_W1' /
124+
&VENT XB = 1, 2, 117, 118, 0, 0, SURF_ID = 'LINEFIRE_W1' /
125+
&VENT XB = 1, 2, 118, 119, 0, 0, SURF_ID = 'LINEFIRE_W1' /
126+
&VENT XB = 2, 3, 118, 119, 0, 0, SURF_ID = 'LINEFIRE_W1' /
127+
&VENT XB = 3, 4, 118, 119, 0, 0, SURF_ID = 'LINEFIRE_W1' /
128+
&VENT XB = 3, 4, 119, 120, 0, 0, SURF_ID = 'LINEFIRE_W1' /
129+
&VENT XB = 4, 5, 119, 120, 0, 0, SURF_ID = 'LINEFIRE_W1' /
130+
&VENT XB = 5, 6, 119, 120, 0, 0, SURF_ID = 'LINEFIRE_W1' /
131+
&VENT XB = 5, 6, 120, 121, 0, 0, SURF_ID = 'LINEFIRE_W1' /
132+
&VENT XB = 6, 7, 120, 121, 0, 0, SURF_ID = 'LINEFIRE_W1' /
133+
&VENT XB = 7, 8, 120, 121, 0, 0, SURF_ID = 'LINEFIRE_W1' /
134+
&VENT XB = 7, 8, 121, 122, 0, 0, SURF_ID = 'LINEFIRE_W1' /
135+
&VENT XB = 8, 9, 121, 122, 0, 0, SURF_ID = 'LINEFIRE_W1' /
136+
&VENT XB = 9, 10, 121, 122, 0, 0, SURF_ID = 'LINEFIRE_W1' /
137+
&VENT XB = 9, 10, 122, 123, 0, 0, SURF_ID = 'LINEFIRE_W1' /
138+
&VENT XB = 10, 11, 122, 123, 0, 0, SURF_ID = 'LINEFIRE_W1' /
139+
&VENT XB = 11, 12, 122, 123, 0, 0, SURF_ID = 'LINEFIRE_W1' /
140+
&VENT XB = 11, 12, 123, 124, 0, 0, SURF_ID = 'LINEFIRE_W1' /
141+
&VENT XB = 12, 13, 123, 124, 0, 0, SURF_ID = 'LINEFIRE_W1' /
142+
&VENT XB = 13, 14, 123, 124, 0, 0, SURF_ID = 'LINEFIRE_W1' /
143+
&VENT XB = 13, 14, 124, 125, 0, 0, SURF_ID = 'LINEFIRE_W1' /
144+
&VENT XB = 14, 15, 124, 125, 0, 0, SURF_ID = 'LINEFIRE_W1' /
145+
&VENT XB = 15, 16, 124, 125, 0, 0, SURF_ID = 'LINEFIRE_W1' /
146+
&VENT XB = 15, 16, 125, 126, 0, 0, SURF_ID = 'LINEFIRE_W1' /
147+
&VENT XB = 16, 17, 125, 126, 0, 0, SURF_ID = 'LINEFIRE_W1' /
148+
&VENT XB = 17, 18, 125, 126, 0, 0, SURF_ID = 'LINEFIRE_W1' /
149+
&VENT XB = 17, 18, 126, 127, 0, 0, SURF_ID = 'LINEFIRE_W1' /
150+
&VENT XB = 18, 19, 126, 127, 0, 0, SURF_ID = 'LINEFIRE_W1' /
151+
&VENT XB = 19, 20, 126, 127, 0, 0, SURF_ID = 'LINEFIRE_W1' /
152+
&VENT XB = 19, 20, 127, 128, 0, 0, SURF_ID = 'LINEFIRE_W1' /
153+
&VENT XB = 20, 21, 127, 128, 0, 0, SURF_ID = 'LINEFIRE_W1' /
154+
&VENT XB = 21, 22, 127, 128, 0, 0, SURF_ID = 'LINEFIRE_W1' /
155+
&VENT XB = 21, 22, 128, 129, 0, 0, SURF_ID = 'LINEFIRE_W1' /
156+
&VENT XB = 22, 23, 128, 129, 0, 0, SURF_ID = 'LINEFIRE_W1' /
157+
&VENT XB = 23, 24, 128, 129, 0, 0, SURF_ID = 'LINEFIRE_W1' /
158+
&VENT XB = 23, 24, 129, 130, 0, 0, SURF_ID = 'LINEFIRE_W1' /
159+
160+
-- grid stretching in z-direction
161+
162+
&TRNZ ID='TRNZ', CC= 0.750000000000000, PC= 0.250000000000000 /
163+
&TRNZ ID='TRNZ', CC= 1.500000000000000, PC= 0.500000000000000 /
164+
&TRNZ ID='TRNZ', CC= 2.250000000000000, PC= 0.750000000000000 /
165+
&TRNZ ID='TRNZ', CC= 3.000000000000000, PC= 1.000000000000000 /
166+
&TRNZ ID='TRNZ', CC= 3.750000000000000, PC= 1.250000000000000 /
167+
&TRNZ ID='TRNZ', CC= 4.500000000000000, PC= 1.500000000000000 /
168+
&TRNZ ID='TRNZ', CC= 5.250000000000000, PC= 1.750000000000000 /
169+
&TRNZ ID='TRNZ', CC= 6.000000000000000, PC= 2.000000000000000 /
170+
&TRNZ ID='TRNZ', CC= 6.750000000000000, PC= 2.266866753288884 /
171+
&TRNZ ID='TRNZ', CC= 7.500000000000000, PC= 2.551738209332682 /
172+
&TRNZ ID='TRNZ', CC= 8.250000000000000, PC= 2.855829091649024 /
173+
&TRNZ ID='TRNZ', CC= 9.000000000000000, PC= 3.180436077523081 /
174+
&TRNZ ID='TRNZ', CC= 9.750000000000000, PC= 3.526943327183480 /
175+
&TRNZ ID='TRNZ', CC= 10.500000000000000, PC= 3.896828386015205 /
176+
&TRNZ ID='TRNZ', CC= 11.250000000000000, PC= 4.291668484977166 /
177+
&TRNZ ID='TRNZ', CC= 12.000000000000000, PC= 4.713147266090125 /
178+
&TRNZ ID='TRNZ', CC= 12.750000000000000, PC= 5.163061961673211 /
179+
&TRNZ ID='TRNZ', CC= 13.500000000000000, PC= 5.643331057942069 /
180+
&TRNZ ID='TRNZ', CC= 14.250000000000000, PC= 6.156002475647093 /
181+
&TRNZ ID='TRNZ', CC= 15.000000000000000, PC= 6.703262302634889 /
182+
&TRNZ ID='TRNZ', CC= 15.750000000000000, PC= 7.287444115569565 /
183+
&TRNZ ID='TRNZ', CC= 16.500000000000000, PC= 7.911038930562729 /
184+
&TRNZ ID='TRNZ', CC= 17.250000000000000, PC= 8.576705825142758 /
185+
&TRNZ ID='TRNZ', CC= 18.000000000000000, PC= 9.287283276856622 /
186+
&TRNZ ID='TRNZ', CC= 18.750000000000000, PC= 10.045801266853291 /
187+
&TRNZ ID='TRNZ', CC= 19.500000000000000, PC= 10.855494200059773 /
188+
&TRNZ ID='TRNZ', CC= 20.250000000000000, PC= 11.719814697042841 /
189+
&TRNZ ID='TRNZ', CC= 21.000000000000000, PC= 12.642448316366462 /
190+
&TRNZ ID='TRNZ', CC= 21.750000000000000, PC= 13.627329270222726 /
191+
&TRNZ ID='TRNZ', CC= 22.500000000000000, PC= 14.678657200349445 /
192+
&TRNZ ID='TRNZ', CC= 23.250000000000000, PC= 15.800915085768803 /
193+
&TRNZ ID='TRNZ', CC= 24.000000000000000, PC= 16.998888358707653 /
194+
&TRNZ ID='TRNZ', CC= 24.750000000000000, PC= 18.277685310211840 /
195+
&TRNZ ID='TRNZ', CC= 25.500000000000000, PC= 19.642758872466420 /
196+
&TRNZ ID='TRNZ', CC= 26.250000000000000, PC= 21.099929870703900 /
197+
&TRNZ ID='TRNZ', CC= 27.000000000000000, PC= 22.655411843849329 /
198+
&TRNZ ID='TRNZ', CC= 27.750000000000000, PC= 24.315837539740158 /
199+
&TRNZ ID='TRNZ', CC= 28.500000000000000, PC= 26.088287197899440 /
200+
&TRNZ ID='TRNZ', CC= 29.250000000000000, PC= 27.980318740463275 /
201+
202+
&TAIL /

0 commit comments

Comments
 (0)