3dfirelab
commented
1 month ago I have notices the same problem. with this simple python script
import fdsreader
import numpy as np
import xarray as xr
import matplotlib.pyplot as plt
sim = fdsreader.Simulation('./')
sim.slices # to show available slice
hrrSlc = sim.slices[2] # to select HRRPUV
hrr, grid = hrrSlc.to_global(masked=True, fill=np.nan, return_coordinates=True)
times = hrrSlc.times
data_array = xr.DataArray(hrr, coords=[times,grid['y'], grid['x']], dims=["time", "y","x"])
data_array.isel(time=5).plot()
plt.show()I can read the 2D HRR slice from a simple burner simulation (see below for config file) run with version FDS6.7.7-4170-g2ce1a34-master
with the same simulation run now with FDS-6.9.1-825-g48a6b66-master the line
hrr, grid =hrrSlc.to_global(masked=True, fill=np.nan, return_coordinates=True)form the python script above does not ouput the hrr with the right dimension.
FDS config file:
&HEAD CHID='fixed_burner', TITLE='Fixed Burner Simulation' /
&TIME T_END=60.0 / ! Simulation runs for 600 seconds
&MESH IJK=10,10,10, XB=0.0,1.5, 0.0,1.5, 0.0,1.5, MULT_ID='mesh'/
&MULT ID='mesh',DX=1.5, DY=1.5, DZ=1.5, I_UPPER=1, J_UPPER=1, K_UPPER=0 / ! Mesh definition
&DUMP DT_HRR=1., DT_SLCF=5., DT_DEVC=5., SIG_FIGS=4, SIG_FIGS_EXP=2, /
&REAC FUEL='CELLULOSE', C=6, H=10, O=5, SOOT_YIELD=0.001, /
&SURF ID='BURNER', COLOR='RED', HRRPUA=1000.0, SPEC_ID='CELLULOSE', RAMP_Q='ramp'/ ! Burner surface properties
&RAMP ID='ramp', T= 0., F=0. /
&RAMP ID='ramp', T= 1., F=1. /
&RAMP ID='ramp', T= 4600., F=1. /
&VENT XB=1.,2.,1.,2.,0.0,0.0, SURF_ID='BURNER' / ! Positioning the burner on the floor
&VENT ID='Mesh Vent: Mesh1 [XMAX]', SURF_ID='OPEN', XB=3.0,3.0,0.0,3.0,0.0,1.5/
&VENT ID='Mesh Vent: Mesh1 [XMIN]', SURF_ID='OPEN', XB=0.0,0.0,0.0,3.0,0.0,1.5/
&VENT ID='Mesh Vent: Mesh1 [YMAX]', SURF_ID='OPEN', XB=0.0,3.0,3.0,3.0,0.0,1.5/
&VENT ID='Mesh Vent: Mesh1 [YMIN]', SURF_ID='OPEN', XB=0.0,3.0,0.0,0.0,0.0,1.5/
&VENT ID='Mesh Vent: Mesh1 [ZMAX]', SURF_ID='OPEN', XB=0.0,3.0,0.0,3.0,1.5,1.5/
&SLCF QUANTITY='TEMPERATURE', PBY=1.5, VECTOR=.TRUE. /
&SLCF QUANTITY='VELOCITY', PBY=1.5, VECTOR=.TRUE. /
&SLCF QUANTITY='HRRPUV', PBZ=0, /
&SLCF QUANTITY='HRRPUV', PBY=1.5, /
&SLCF QUANTITY='TEMPERATURE', PBZ=0, /
&BNDF QUANTITY='RADIATIVE HEAT FLUX' /
&BNDF QUANTITY='WALL TEMPERATURE' /
&TAIL /
lu-kas
Hi,
I am new to FDS.
I am encountering an error while recreating the AEST contour map from the tutorial on the following website: https://firedynamics.github.io/LectureFireSimulation/content/examples/02_intermediate/05_data_analysis_02.html#../../tools/03_analysis/02_fdsreaders
The error message is:
WARNING:root:Module vents: could not convert string to float: '!' The error can be safely ignored if not requiring the vents module. However, please consider to submit an issue on Github including the error message, the stack trace and your FDS input-file so we can reproduce the error and fix it as soon as possible!
The Python code correctly captures the domain, but this is not true for the data. Please find attached screenshots for the expected and actual results.
Expected:
Real:
Could this be due to different versions of FDS? The tutorial was created using FDS 6.7, while I am using the latest 6.9 version. I believe python is not able to get the data for the correct time step. It has only data for time step = 0 and the problem is with the following line of the code:
it = slc.get_nearest_timestep(50)