Velocity Slice values Issues

[2000-lang]

Hello,

merged_u_velocity_200_300s.csv

This is my velocity slice extracted data for the below code. I have few questions,

1. Why I got velocity spike at mesh boundries?
2. There is a lot of soot irregularly produced , and less fire seen in smv
3. Also, I am continuosly getting mesh alignement errors. User manual also was not completely helpful.

Can you please tell me best mesh defining for my geometry and requirement? Also is there any references of inputs in these type of cases?

I need 0.25m resolution for fine mesh and 2m for coarse mesh.

&HEAD CHID='street_canyon_treefiretestHW1',
TITLE='Street canyon with three burning particle trees testHW1' /

!------------------ MESHES ------------------

! Fine mesh 1: x = 112–117, y = 35–65
&MESH ID='FINE_TREES_1',
XB = 112.0,118.0, 24.0,66.0, 0.0,20.0,
IJK = 24,168,40,
MPI_PROCESS=1,
CHECK_MESH_ALIGNMENT=.FALSE. /

! Fine mesh 2: x = 117–122, y = 35–65
&MESH ID='FINE_TREES_2',
XB = 118.0,122.0, 24.0,66.0, 0.0,20.0,
IJK = 16,168,40,
MPI_PROCESS=2,
CHECK_MESH_ALIGNMENT=.FALSE. /

! Fine mesh 3: x = 112–117, y = 65–95
&MESH ID='FINE_TREES_3',
XB = 112.0,118.0, 66.0,102.0, 0.0,20.0,
IJK = 24,144,40,
MPI_PROCESS=3,
CHECK_MESH_ALIGNMENT=.FALSE. /

! Fine mesh 4: x = 117–122, y = 65–95
&MESH ID='FINE_TREES_4',
XB = 118.0,122.0, 66.0,102.0, 0.0,20.0,
IJK = 16,144,40,
MPI_PROCESS=4,
CHECK_MESH_ALIGNMENT=.FALSE. /

! Coarse resolution ~2 m
! X is split at 126, Y at 63 (2×2 tiling)

! Lower-left coarse (x: 0–126, y: 0–63)
&MESH ID='COARSE_LL',
XB = 0.0,112.0, 0.0,24.0, 0.0,90.0,
IJK = 56,12,180,
MPI_PROCESS=5,
CHECK_MESH_ALIGNMENT=.FALSE. /

! Lower-right coarse (x: 126–252, y: 0–63)
&MESH ID='COARSE_LR',
XB = 122.0,252.0, 0.0,24.0, 0.0,90.0,
IJK = 65,12,180,
MPI_PROCESS=6,
CHECK_MESH_ALIGNMENT=.FALSE. /

! Upper-left coarse (x: 0–126, y: 63–126)
&MESH ID='COARSE_UL',
XB = 0.0,112.0, 102.0,128.0, 0.0,90.0,
IJK = 56,13,180,
MPI_PROCESS=7,
CHECK_MESH_ALIGNMENT=.FALSE. /

! Upper-right coarse (x: 126–252, y: 63–126)
&MESH ID='COARSE_UR',
XB = 122.0,252.0, 102.0,128.0, 0.0,90.0,
IJK = 65,13,180,
MPI_PROCESS=8,
CHECK_MESH_ALIGNMENT=.FALSE. /

&MISC TMPA=293.0, CFL_MAX=0.8, RESTART=.FALSE. /

&TIME T_END=300.0 /

&SURF ID='INERT' /

&VENT PBX=0.0, SURF_ID='OPEN' /
&VENT PBX=252.0, SURF_ID='OPEN' /
&VENT PBY=0.0, SURF_ID='OPEN' /
&VENT PBY=128.0, SURF_ID='OPEN' /
&VENT PBZ=90.0, SURF_ID='OPEN' /

&WIND SPEED=6.0, Z_0=0.5, Z_REF=5.0, DIRECTION=270.0 /

!------------ WOOD FUEL AND PRODUCTS (User Guide) ------------

&SPEC ID='WOOD', FORMULA='CH1.7O0.74N0.002' /
&SPEC ID='OXYGEN', LUMPED_COMPONENT_ONLY=T /
&SPEC ID='NITROGEN', LUMPED_COMPONENT_ONLY=T /
&SPEC ID='WATER VAPOR', LUMPED_COMPONENT_ONLY=T /
&SPEC ID='CARBON MONOXIDE', LUMPED_COMPONENT_ONLY=T /
&SPEC ID='CARBON DIOXIDE', LUMPED_COMPONENT_ONLY=T /
&SPEC ID='SOOT', LUMPED_COMPONENT_ONLY=T /

&SPEC ID='AIR',
SPEC_ID(1)='OXYGEN', VOLUME_FRACTION(1)=1.,
SPEC_ID(2)='NITROGEN', VOLUME_FRACTION(2)=3.76,
BACKGROUND=T /

&SPEC ID='PRODUCTS_2',
SPEC_ID(1)='CARBON DIOXIDE', VOLUME_FRACTION(1)=0.96438,
SPEC_ID(2)='CARBON MONOXIDE',VOLUME_FRACTION(2)=0.00365,
SPEC_ID(3)='WATER VAPOR', VOLUME_FRACTION(3)=0.85,
SPEC_ID(4)='SOOT', VOLUME_FRACTION(4)=0.03196,
SPEC_ID(5)='NITROGEN', VOLUME_FRACTION(5)=3.84076 /

&REAC ID='WOOD_FIRE',
FUEL='WOOD',
HEAT_OF_COMBUSTION=16400.,
SPEC_ID_NU='WOOD','AIR','PRODUCTS_2',
NU=-1.,-1.02121,1. /

!------------ MATERIALS ------------

&MATL ID='MOISTURE',
DENSITY=1000.,
CONDUCTIVITY=0.6,
SPECIFIC_HEAT=4.19 /

&MATL ID='DRY_PINE',
DENSITY=500.,
CONDUCTIVITY=0.1,
SPECIFIC_HEAT=1.0 /

!------------ TRUNK (OBST, WOOD, simple burning) ------------

&SURF ID='TRUNK_WOOD',
MATL_ID(1,1)='DRY_PINE',
MATL_MASS_FRACTION(1,1)=1.0,
THICKNESS(1)=0.20,
HRRPUA=150.0,
SPEC_ID='WOOD',
COLOR='BROWN' /

!----------------- BUILDINGS (unchanged) -----------------
&OBST XB=90.0,108.0, 18.0,36.0, 0.0,18.0, SURF_ID='INERT' /
&OBST XB=90.0,108.0, 54.0,72.0, 0.0,18.0, SURF_ID='INERT' /
&OBST XB=90.0,108.0, 90.0,108.0, 0.0,18.0, SURF_ID='INERT' /
&OBST XB=126.0,144.0, 18.0,36.0, 0.0,18.0, SURF_ID='INERT' /
&OBST XB=126.0,144.0, 54.0,72.0, 0.0,18.0, SURF_ID='INERT' /
&OBST XB=126.0,144.0, 90.0,108.0, 0.0,18.0, SURF_ID='INERT' /

! Three trunks, same geometry as single_tree_fire (0.2 m x 0.2 m, 2 m high)
&OBST ID='TREE1_TRUNK',
XB=116.5,117.5, 26.5,27.5, 0.0,4.0,
SURF_ID='TRUNK_WOOD' /

&OBST ID='TREE2_TRUNK',
XB=116.5,117.5, 62.5,63.5, 0.0,4.0,
SURF_ID='TRUNK_WOOD' /

&OBST ID='TREE3_TRUNK',
XB=116.5,117.5, 98.5,99.5, 0.0,4.0,
SURF_ID='TRUNK_WOOD' /

!------------ CROWN NEEDLE SURFACE (WOOD FIRE SOURCE) ------------

&SURF ID='needle surface',
MATL_ID(1,1)='DRY_PINE',
MATL_MASS_FRACTION(1,1)=1.0,
MOISTURE_FRACTION=0.25,
SURFACE_VOLUME_RATIO=8000.,
LENGTH=0.10,
GEOMETRY='CYLINDRICAL',
HRRPUA=150.0,
SPEC_ID='WOOD',
BURN_AWAY=.TRUE.,
COLOR='FOREST GREEN' /

!------------ PARTICLES: PINE NEEDLES (STATIC, BURNING) ------------

&PART ID='needles',
SURF_ID='needle surface',
DRAG_COEFFICIENT=2.8,
STATIC=.TRUE.,
COLOR='FOREST GREEN' /

&PROP ID='needle image',
SMOKEVIEW_ID='TUBE',
SMOKEVIEW_PARAMETERS='L=0.1,D=0.0005' /

!------------ TREE CROWN GEOMETRY (CONE) ------------

&INIT ID='TREE1_CROWN',
PART_ID='needles',
SHAPE='CONE',
XYZ=117.0,27.0,0.0,
CROWN_BASE_WIDTH=3.0,
CROWN_BASE_HEIGHT=4.0,
TREE_HEIGHT=8.0,
N_PARTICLES_PER_CELL=6,
MASS_PER_VOLUME=6, ! was 6.0 (more fuel, longer burn)
DRY=.TRUE. /

&INIT ID='TREE2_CROWN',
PART_ID='needles',
SHAPE='CONE',
XYZ=117.0,63.0,0.0,
CROWN_BASE_WIDTH=3.0,
CROWN_BASE_HEIGHT=4.0,
TREE_HEIGHT=8.0,
N_PARTICLES_PER_CELL=6,
MASS_PER_VOLUME=6,
DRY=.TRUE. /

&INIT ID='TREE3_CROWN',
PART_ID='needles',
SHAPE='CONE',
XYZ=117.0,99.0,0.0,
CROWN_BASE_WIDTH=3.0,
CROWN_BASE_HEIGHT=4.0,
TREE_HEIGHT=8.0,
N_PARTICLES_PER_CELL=6,
MASS_PER_VOLUME=6,
DRY=.TRUE. /

!------------ OUTPUT ------------

&DUMP DT_HRR=5.0, DT_SLCF=0.5, DT_RESTART=60.0 /

&SLCF QUANTITY='U-VELOCITY', PBY=27.0 /
&SLCF QUANTITY='U-VELOCITY', PBY=63.0 /
&SLCF QUANTITY='U-VELOCITY', PBY=99.0 /
&SLCF QUANTITY='U-VELOCITY', PBX=117.0 /
&SLCF QUANTITY='U-VELOCITY', PBZ=6.0 /
&SLCF QUANTITY='U-VELOCITY', PBZ=9.0 /

&SLCF QUANTITY='V-VELOCITY', PBY=27.0 /
&SLCF QUANTITY='V-VELOCITY', PBY=63.0 /
&SLCF QUANTITY='V-VELOCITY', PBY=99.0 /
&SLCF QUANTITY='V-VELOCITY', PBX=117.0 /
&SLCF QUANTITY='V-VELOCITY', PBZ=6.0 /
&SLCF QUANTITY='V-VELOCITY', PBZ=9.0 /

&SLCF QUANTITY='W-VELOCITY', PBY=27.0 /
&SLCF QUANTITY='W-VELOCITY', PBY=63.0 /
&SLCF QUANTITY='W-VELOCITY', PBY=99.0 /
&SLCF QUANTITY='W-VELOCITY', PBX=117.0 /
&SLCF QUANTITY='W-VELOCITY', PBZ=6.0 /
&SLCF QUANTITY='W-VELOCITY', PBZ=9.0 /

&SLCF QUANTITY='TEMPERATURE', PBY=27.0 /
&SLCF QUANTITY='TEMPERATURE', PBY=63.0 /
&SLCF QUANTITY='TEMPERATURE', PBY=99.0 /
&SLCF QUANTITY='TEMPERATURE', PBX=117.0 /
&SLCF QUANTITY='TEMPERATURE', PBZ=9.0 /
&SLCF QUANTITY='TEMPERATURE', PBZ=6.0 /

&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT', PBY=27.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT', PBY=63.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT', PBY=99.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT', PBX=117.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT', PBZ=6.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT', PBZ=9.0 /

&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', PBY=27.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', PBY=63.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', PBY=99.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', PBX=117.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', PBZ=6.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', PBZ=9.0 /

&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', PBY=27.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', PBY=63.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', PBY=99.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', PBX=117.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', PBZ=6.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', PBZ=9.0 /

&TAIL /

[drjfloyd]

This is what I see when I run your input file. You have disconnected meshes with solid faces and are specifying a wind. You should be meshing the entire extents of your geometry. Fix that and try running again.

[2000-lang]

Thanks for the response. Actually, I am struck at the post processing step of the simulation. My requirement is to extract the instantaneous slice data. But the user manual shows how to extract the average not the instantaneous ones. So, I am literally confused how to do it for a single mesh and multi mesh. Any leads would be really helpful.

[mcgratta]

You could try averaging the data over a small segment in time. That would essentially be the "instantaneous".

[2000-lang]

Yes that is an option. But in large simulations and multi mesh cases like ex: For 600 seconds and my input has 8 meshes , if I want to extract slice data for a plane for every second, it is really difficult right ? And also each slice has 8 different divisions and merging all the data is not practically accurate right ?

[mcgratta]

Sounds to me that you probably should write your own parser. The format of the slice files is given in the back of the User's Guide. I wrote the fds2ascii routine a long time ago and it has not kept up with the changes in meshing. Because everyone seems to have a different need for reading and parsing those files, we generally leave it to the users. This is actually a good problem for an LLM. Ask one to write a Python (for example) routine to read the file and do whatever it is that you want to do. We could never write a single routine that would do everything that everyone would want.

[drjfloyd]

There are a couple of existing python tools that other users have developed near the end of:

https://github.com/firemodels/fds/wiki/Third-Party-Tools

[2000-lang]

Is this good now ? Can you please check with this input

! =====================================================
! FINE REGION (0.25 m)
! X=112–126, Y=34–94, Z=0–30
! Split in Y only (2 meshes)
! =====================================================

&MESH ID='F1',
XB=112.0,126.0, 34.0,64.0, 0.0,30.0,
IJK=56,120,15 /

&MESH ID='F2',
XB=112.0,126.0, 64.0,94.0, 0.0,30.0,
IJK=56,120,15 /

! =====================================================
! COARSE REGION (2 m)
! Entire domain split at:
! X = 112, 126
! Y = 34, 94
! Z = 30
! =====================================================

! ---------------- LEFT SIDE (0–112) ----------------

! Bottom layer Z=0–30
&MESH XB=0.0,112.0, 0.0,34.0, 0.0,30.0, IJK=56,17,15 /
&MESH XB=0.0,112.0, 34.0,94.0, 0.0,30.0, IJK=56,30,15 /
&MESH XB=0.0,112.0, 94.0,128.0,0.0,30.0, IJK=56,17,15 /

! Top layer Z=30–90
&MESH XB=0.0,112.0, 0.0,34.0, 30.0,90.0, IJK=56,17,30 /
&MESH XB=0.0,112.0, 34.0,94.0,30.0,90.0, IJK=56,30,30 /
&MESH XB=0.0,112.0, 94.0,128.0,30.0,90.0, IJK=56,17,30 /

! ---------------- RIGHT SIDE (126–252) ----------------

! Bottom layer Z=0–30
&MESH XB=126.0,252.0, 0.0,34.0, 0.0,30.0, IJK=63,17,15 /
&MESH XB=126.0,252.0, 34.0,94.0, 0.0,30.0, IJK=63,30,15 /
&MESH XB=126.0,252.0, 94.0,128.0,0.0,30.0, IJK=63,17,15 /

! Top layer Z=30–90
&MESH XB=126.0,252.0, 0.0,34.0, 30.0,90.0, IJK=63,17,30 /
&MESH XB=126.0,252.0, 34.0,94.0,30.0,90.0, IJK=63,30,30 /
&MESH XB=126.0,252.0, 94.0,128.0,30.0,90.0, IJK=63,17,30 /

! ---------------- CENTER FRONT/BACK (112–126 outside fine) ----------------

! Bottom strips Z=0–30
&MESH XB=112.0,126.0, 0.0,34.0, 0.0,30.0, IJK=7,17,15 /
&MESH XB=112.0,126.0, 94.0,128.0,0.0,30.0, IJK=7,17,15 /

! Top full column Z=30–90
&MESH XB=112.0,126.0, 0.0,128.0,30.0,90.0, IJK=7,64,30 /

&TAIL /

[Er9y714]

Note that TMPA in MISC line is in Celsius.
&MISC TMPA=293.0,

Meshes are aligning but note that the first two meshes have cell size ratio (aspect ratio) of 8. It is recommended not to have higher than 2 to avoid numerical instability according to User guide 4.2 Common error statements.

The meshes adjacent to the fine mesh has a transition of a factor of 8. This is recommended not to be larger than 2.

[2000-lang]

Can you please comment on this. I tried real hard to define the mesh with proper aspect ratio. Thanks for the guidance again.

&HEAD CHID='mesh_test', TITLE='Mesh only test' /

! =====================================================
! 0.25 m FIRE REGION | Z = 0–20
! =====================================================
&MESH ID='FIRE1', IJK=72,140,40, XB=108,126,24,59, 0,10 /
&MESH ID='FIRE2', IJK=72,140,40, XB=108,126,59,94, 0,10 /
&MESH ID='FIRE3', IJK=72,140,40, XB=108,126,24,59,10,20 /
&MESH ID='FIRE4', IJK=72,140,40, XB=108,126,59,94,10,20 /

! =====================================================
! 0.5 m BUFFER RING | Z = 0–20
! =====================================================
&MESH ID='BUF_L', IJK=8,156,40, XB=104,108,20,98, 0,20 /
&MESH ID='BUF_R', IJK=8,156,40, XB=126,130,20,98, 0,20 /
&MESH ID='BUF_F', IJK=36,8,40, XB=108,126,20,24, 0,20 /
&MESH ID='BUF_B', IJK=36,8,40, XB=108,126,94,98, 0,20 /

! =====================================================
! 1 m TRANSITION RING | Z = 0–20
! =====================================================
&MESH ID='TR_L_Z0', IJK=8,94,20, XB= 96,104,12,106, 0,20 /
&MESH ID='TR_R_Z0', IJK=8,94,20, XB=130,138,12,106, 0,20 /
&MESH ID='TR_F_Z0', IJK=26,8,20, XB=104,130,12, 20, 0,20 /
&MESH ID='TR_B_Z0', IJK=26,8,20, XB=104,130,98,106, 0,20 /

! =====================================================
! 0.5 m PLUME LAYER | Z = 20–30
! =====================================================
&MESH ID='PLM05_1', IJK=52,78,20, XB=104,130,20,59,20,30 /
&MESH ID='PLM05_2', IJK=52,78,20, XB=104,130,59,98,20,30 /

! =====================================================
! 1 m TRANSITION RING UPPER | Z = 20–40
! =====================================================
&MESH ID='TR_L_Z20', IJK=8,94,20, XB= 96,104,12,106,20,40 /
&MESH ID='TR_R_Z20', IJK=8,94,20, XB=130,138,12,106,20,40 /
&MESH ID='TR_F_Z20', IJK=26,8,20, XB=104,130,12, 20,20,40 /
&MESH ID='TR_B_Z20', IJK=26,8,20, XB=104,130,98,106,20,40 /

! =====================================================
! 1 m PLUME LAYER | Z = 30–40
! =====================================================
&MESH ID='PLM10', IJK=26,78,10, XB=104,130,20,98,30,40 /

! =====================================================
! 2 m OUTER DOMAIN | Z = 0–90 / 40–90
! =====================================================
&MESH ID='OUT_L1', IJK=48,32,45, XB= 0, 96, 0, 64,0,90 /
&MESH ID='OUT_L2', IJK=48,32,45, XB= 0, 96, 64,128,0,90 /
&MESH ID='OUT_R1', IJK=57,32,45, XB=138,252, 0, 64,0,90 /
&MESH ID='OUT_R2', IJK=57,32,45, XB=138,252, 64,128,0,90 /
&MESH ID='OUT_F', IJK=21,6,45, XB= 96,138, 0, 12,0,90 /
&MESH ID='OUT_B', IJK=21,11,45, XB= 96,138,106,128,0,90 /
&MESH ID='OUT_TC', IJK=21,47,25, XB= 96,138, 12,106,40,90 /

&MISC TMPA=20.0 /

&TIME T_END=10.0 /

&VENT PBX=0.0, SURF_ID='OPEN' /
&VENT PBX=252.0, SURF_ID='OPEN' /
&VENT PBY=0.0, SURF_ID='OPEN' /
&VENT PBY=128.0, SURF_ID='OPEN' /
&VENT PBZ=90.0, SURF_ID='OPEN' /

&TAIL /

[Er9y714]

In the latest mesh the cells are cubic and my previous remarks were resolved. It looks better.

[2000-lang]

Thanks for the feedback. But when I run this its crashing. I am not able to figure out what is wrong below:

&HEAD CHID='mesh_test', TITLE='Mesh only test' /

! =====================================================
! 0.25 m FIRE REGION | Z = 0–20
! =====================================================
&MESH ID='FIRE1', IJK=72,140,40, XB=108,126,24,59, 0,10 /
&MESH ID='FIRE2', IJK=72,140,40, XB=108,126,59,94, 0,10 /
&MESH ID='FIRE3', IJK=72,140,40, XB=108,126,24,59,10,20 /
&MESH ID='FIRE4', IJK=72,140,40, XB=108,126,59,94,10,20 /

! =====================================================
! 0.5 m BUFFER RING | Z = 0–20
! =====================================================
&MESH ID='BUF_L', IJK=8,156,40, XB=104,108,20,98, 0,20 /
&MESH ID='BUF_R', IJK=8,156,40, XB=126,130,20,98, 0,20 /
&MESH ID='BUF_F', IJK=36,8,40, XB=108,126,20,24, 0,20 /
&MESH ID='BUF_B', IJK=36,8,40, XB=108,126,94,98, 0,20 /

! =====================================================
! 1 m TRANSITION RING | Z = 0–20
! =====================================================
&MESH ID='TR_L_Z0', IJK=8,94,20, XB= 96,104,12,106, 0,20 /
&MESH ID='TR_R_Z0', IJK=8,94,20, XB=130,138,12,106, 0,20 /
&MESH ID='TR_F_Z0', IJK=26,8,20, XB=104,130,12, 20, 0,20 /
&MESH ID='TR_B_Z0', IJK=26,8,20, XB=104,130,98,106, 0,20 /

! =====================================================
! 0.5 m PLUME LAYER | Z = 20–30
! =====================================================
&MESH ID='PLM05_1', IJK=52,78,20, XB=104,130,20,59,20,30 /
&MESH ID='PLM05_2', IJK=52,78,20, XB=104,130,59,98,20,30 /

! =====================================================
! 1 m TRANSITION RING UPPER | Z = 20–40
! =====================================================
&MESH ID='TR_L_Z20', IJK=8,94,20, XB= 96,104,12,106,20,40 /
&MESH ID='TR_R_Z20', IJK=8,94,20, XB=130,138,12,106,20,40 /
&MESH ID='TR_F_Z20', IJK=26,8,20, XB=104,130,12, 20,20,40 /
&MESH ID='TR_B_Z20', IJK=26,8,20, XB=104,130,98,106,20,40 /

! =====================================================
! 1 m PLUME LAYER | Z = 30–40
! =====================================================
&MESH ID='PLM10', IJK=26,78,10, XB=104,130,20,98,30,40 /

! =====================================================
! 2 m OUTER DOMAIN | Z = 0–90 / 40–90
! =====================================================
&MESH ID='OUT_L1', IJK=48,32,45, XB= 0, 96, 0, 64,0,90 /
&MESH ID='OUT_L2', IJK=48,32,45, XB= 0, 96, 64,128,0,90 /
&MESH ID='OUT_R1', IJK=57,32,45, XB=138,252, 0, 64,0,90 /
&MESH ID='OUT_R2', IJK=57,32,45, XB=138,252, 64,128,0,90 /
&MESH ID='OUT_F', IJK=21,6,45, XB= 96,138, 0, 12,0,90 /
&MESH ID='OUT_B', IJK=21,11,45, XB= 96,138,106,128,0,90 /
&MESH ID='OUT_TC', IJK=21,47,25, XB= 96,138, 12,106,40,90 /

! =====================================================
! STANDARD WOOD FUEL (pro duction-tested, proper soot)
! =====================================================

&SPEC ID='WOOD', FORMULA='CH1.7O0.74N0.002', MW=24.5 /
&SPEC ID='CARBON MONOXIDE', FORMULA='CO', MW=28.01 /
&SPEC ID='CARBON DIOXIDE', FORMULA='CO2', MW=44.01 /
&SPEC ID='SOOT', FORMULA='C', MW=12.0 /

&REAC ID='WOOD_FIRE', FUEL='WOOD', C=1.0, H=1.7, N=0.002, O=0.74, SOOT_YIELD=0.004, CO_YIELD=0.025 /

&SURF ID='INERT', RGB=128,128,128 /
&SURF ID='BUILDING_WALL', RGB=255,255,255 /

&SURF ID='TRUNK_WOOD',
MATL_ID(1,1)='DRY_PINE',
MATL_MASS_FRACTION(1,1)=1.0,
THICKNESS(1)=0.20,
HRRPUA=600.0,
SPEC_ID='WOOD',
COLOR='BROWN' /

&MATL ID='MOISTURE',
DENSITY=1000.,
CONDUCTIVITY=0.6,
SPECIFIC_HEAT=4.19 /

&MATL ID='DRY_PINE',
DENSITY=500.,
CONDUCTIVITY=0.1,
SPECIFIC_HEAT=1.0 /

! =====================================================
! WIND & VENTS
! =====================================================

&MISC TMPA=20.0 /

&TIME T_END=10. /

&WIND SPEED=6.0, Z_0=0.5, Z_REF=5.0, DIRECTION=270.0 /

&VENT PBX=0.0, SURF_ID='OPEN' /
&VENT PBX=252.0, SURF_ID='OPEN' /
&VENT PBY=0.0, SURF_ID='OPEN' /
&VENT PBY=128.0, SURF_ID='OPEN' /
&VENT PBZ=90.0, SURF_ID='OPEN' /

! =====================================================
! BUILDINGS
! =====================================================

&OBST XB=90.0,108.0, 18.0,108.0, 0.0,18.0, SURF_ID='BUILDING_WALL' /
&OBST XB=126.0,144.0, 18.0,108.0, 0.0,18.0, SURF_ID='BUILDING_WALL' /

! =====================================================
! TREE TRUNKS & CROWNS (Proper woo d fire)
! =====================================================

! Trunks
&OBST ID='TREE1_TRUNK', XB=116.9,117.9, 33.9,34.9, 0.0,4.0, SURF_ID='TRUNK_WOOD' /
&OBST ID='TREE2_TRUNK', XB=116.9,117.9, 62.9,63.9, 0.0,4.0, SURF_ID='TRUNK_WOOD' /
&OBST ID='TREE3_TRUNK', XB=116.9,117.9, 91.9,92.9, 0.0,4.0, SURF_ID='TRUNK_WOOD' /

&SURF ID='needle surface',
MATL_ID(1,1)='DRY_PINE',
MATL_MASS_FRACTION(1,1)=1.0,
MOISTURE_FRACTION=0.12,
SURFACE_VOLUME_RATIO=8000.,
LENGTH=0.10,
GEOMETRY='CYLINDRICAL',
HRRPUA=150.0,
SPEC_ID='WOOD',
BURN_AWAY=.TRUE.,
COLOR='FOREST GREEN' /

&PART ID='needles',
SURF_ID='needle surface',
DRAG_COEFFICIENT=2.8,
STATIC=.TRUE.,
COLOR='FOREST GREEN' /

&PROP ID='needle image',
SMOKEVIEW_ID='TUBE',
SMOKEVIEW_PARAMETERS='L=0.1,D=0.0005' /

! Three particle crowns, identical geometry/physics to single_tree_fire,
! just shifted into the canyon
&INIT ID='TREE1_CROWN',
PART_ID='needles',
SHAPE='CONE',
XYZ=117.0,34.0,0.0,
CROWN_BASE_WIDTH=3.0,
CROWN_BASE_HEIGHT=4.0,
TREE_HEIGHT=8.0,
N_PARTICLES_PER_CELL=1,
MASS_PER_VOLUME=6.0,
DRY=.TRUE. /

&INIT ID='TREE2_CROWN',
PART_ID='needles',
SHAPE='CONE',
XYZ=117.0,63.0,0.0,
CROWN_BASE_WIDTH=3.0,
CROWN_BASE_HEIGHT=4.0,
TREE_HEIGHT=8.0,
N_PARTICLES_PER_CELL=1,
MASS_PER_VOLUME=6.0,
DRY=.TRUE. /

&INIT ID='TREE3_CROWN',
PART_ID='needles',
SHAPE='CONE',
XYZ=117.0,92.0,0.0,
CROWN_BASE_WIDTH=3.0,
CROWN_BASE_HEIGHT=4.0,
TREE_HEIGHT=8.0,
N_PARTICLES_PER_CELL=1,
MASS_PER_VOLUME=6.0,
DRY=.TRUE. /

!------------ OUTPUT ------------

&DUMP DT_HRR=5.0, DT_SLCF=1.0 /

&SLCF QUANTITY='U-VELOCITY', PBY=27.0 /
&SLCF QUANTITY='U-VELOCITY', PBY=63.0 /
&SLCF QUANTITY='U-VELOCITY', PBY=99.0 /
&SLCF QUANTITY='U-VELOCITY', PBX=117.0 /
&SLCF QUANTITY='U-VELOCITY', PBZ=6.0 /
&SLCF QUANTITY='U-VELOCITY', PBZ=9.0 /

&SLCF QUANTITY='V-VELOCITY', PBY=27.0 /
&SLCF QUANTITY='V-VELOCITY', PBY=63.0 /
&SLCF QUANTITY='V-VELOCITY', PBY=99.0 /
&SLCF QUANTITY='V-VELOCITY', PBX=117.0 /
&SLCF QUANTITY='V-VELOCITY', PBZ=6.0 /
&SLCF QUANTITY='V-VELOCITY', PBZ=9.0 /

&SLCF QUANTITY='W-VELOCITY', PBY=27.0 /
&SLCF QUANTITY='W-VELOCITY', PBY=63.0 /
&SLCF QUANTITY='W-VELOCITY', PBY=99.0 /
&SLCF QUANTITY='W-VELOCITY', PBX=117.0 /
&SLCF QUANTITY='W-VELOCITY', PBZ=6.0 /
&SLCF QUANTITY='W-VELOCITY', PBZ=9.0 /

&SLCF QUANTITY='TEMPERATURE', PBY=27.0 /
&SLCF QUANTITY='TEMPERATURE', PBY=63.0 /
&SLCF QUANTITY='TEMPERATURE', PBY=99.0 /
&SLCF QUANTITY='TEMPERATURE', PBX=117.0 /
&SLCF QUANTITY='TEMPERATURE', PBZ=9.0 /
&SLCF QUANTITY='TEMPERATURE', PBZ=6.0 /

&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT', PBY=27.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT', PBY=63.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT', PBY=99.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT', PBX=117.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT', PBZ=6.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT', PBZ=9.0 /

&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', PBY=27.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', PBY=63.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', PBY=99.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', PBX=117.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', PBZ=6.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', PBZ=9.0 /

&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', PBY=27.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', PBY=63.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', PBY=99.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', PBX=117.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', PBZ=6.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', PBZ=9.0 /

&TAIL /

SLURM SCRIPT:

#!/bin/bash
#SBATCH --job-name=simple_test
#SBATCH --output=simple_test.out
#SBATCH --error=simple_test.err
#SBATCH --ntasks=26
#SBATCH --cpus-per-task=1
#SBATCH --time=90:00:00
#SBATCH --partition=long
#SBATCH --mem=64G

module purge

module load tbb

module load compiler-rt/2024.2.1

module load compiler-intel-llvm/2024.2.1

module load mpi/2021.13

export I_MPI_FABRICS=shm:
export FI_PROVIDER=tcp

mpirun -n 26 ./fds_impi_intel_linux finaltest.fds > fds_output.log

[drjfloyd]

You input ran with the current source. You can try running with the current test bundle at https://github.com/firemodels/test_bundles

A few comments on your input file:

Delete your &SPEC lines. Once you define C,H,O, and N on REAC all those species get defined by FDS. I suggest you review the User's Guide in particular the Chapters on defining reactions and defining species and look at the various example.

Before running Delete your &DUMP input so you get more output over the 10 s you are running the case. Look at your hrr in the hrr.csv. Is that the hrr you expect?

Make sure you pay attention to units for inputs. The User's Guide gives the units for all FDS inputs. For example here you have defined the initial temperature to be 293 C.

&MISC TMPA=293.0, CFL_MAX=0.8, RESTART=.FALSE. /