lautenberger / elmfire

Eulerian Level set Model of FIRE spread
https://elmfire.io
Eclipse Public License 2.0
23 stars 11 forks source link

Non uniform fireline intensity when activating spotting model #19

Closed yqin123 closed 9 months ago

yqin123 commented 1 year ago

Hi Chris,

Sorry for the late post. As discussed before, when I apply uniform fuel type (102) and constant wind speed (north to south, 15 mph), the fireline intensity is not uniform when I activate the spotting model. When I deactivate the spotting model, this issue of non-uniform occurs at the boundaries of the simulation region. I suspect that it could be a problem due to my enforcing the value of the level-set variable at spotting ignited pixels to be -1 at every time step. Because it happens at the boundary in the no-spotting cases where the level-set variable is also enforced to be -1. But I am not sure how exactly it could impact the fireline intensity calculation. I attached the input files I used to run the cases with/without spotting. I pulled the newest version of ELMFIRE and used it to test the cases.test_cases.zip

Thanks, Yiren

lautenberger commented 1 year ago

Thank you Yiren, I will take a look after the 4th of July holiday!

Your initial thinking is that this is caused by forcing the level set function to -1 in cells that are ignited by spotting?

lautenberger commented 12 months ago

@yqin123 I tried running the input decks in test_cases.zip but it seems that there is not a complete set of inputs? Could you either 1) create and post a link to an updated zip file with all necessary inputs, or 2) Update the scripts to generate the inputs on the fly and test it in an empty directory? I'll look into this after that. Thanks!

yqin123 commented 11 months ago

Hi Chris,

Sorry for the late feedback.

Yes, my guess is the issue may be caused by forcing the level set function to -1 in cells that are ignited by spotting. I remember in the original ELMFIRE algorithm, the values of phi at pixels ignited by spotting are not forced to -1 at every time step. But I am not sure about it, because I did not see the same issue in my 1-D matlab tests.

I got some trouble downloading results from the server at UMD recently. I re-implemented ELMFIRE locally these days and here are the input files. (raster inputs and elmfire.data inputs are all given in the directory 'inputs'). test_cases.zip

Thanks, Yiren

yqin123 commented 9 months ago

Hi Chris,

The problem is solved. The fireline intensity is calculated from the spreading velocity of the fire. The velocity is not the maximum velocity predicted by the Rothermel model but the projected velocity normal to the fire perimeter. Therefore, the fireline intensity is different in all directions when affected by wind.

When firebrand showering is activated in ELMFIRE, independent spot fires will create non-uniform fireline intensity around them, and for the linear fire spreading without embers, fire spreading towards the edges also has a different fireline intensity due to the same issue.

Thanks, Yiren